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In september nam prof. dr. J. de Boer, hoogleraar in de Theoretische Natuurkunde te 
Amsterdam afscheid, wegens het bereiken van de pensioengerechtigde leeftijd. 
Prof. De Boer kon als een coryfee, in het wereldje van Nederlandse natuurkundigen 


beschouwd worden. 


Bij velen zal De Boer bekend zijn vanwege de talrijke organisatorische functies die hij 
heeft bekleed. Vooral tijdens zijn jarenlang voorzitterschap van het bestuur van de 
Stichting FOM heeft hij metterdaad wetenschappelijk onderzoek aan universiteiten, 
hogescholen en instituten gestimuleerd. Het NTvN-A is echter niet de plaats om 


daarop verder in te gaan. 


Ten einde zijn wetenschappelijke verdiensten nog eens te benadrukken werd op 17 

en 18 september jl. in Amsterdam een symposium georganiseerd waar leerlingen en 
vrienden van De Boer een overzicht geven van de nieuwste ontwikkelingen op hun vak- 
gebied. Ontwikkelingen die veelal zijn voortgesproten of in ieder geval nauw samen- 
hangen met het werk van De Boer en zijn medewerkers. De sprekers waren bereid hun 
tekst voor publikatie te bewerken en het voor u liggende nummer is het resultaat van 
hun inspanning. Het wijkt in zoverre af van wat de redactie in het algemeen nastreeft, 
dat het meer als een symposiumwverslag moet worden gezien dan als de beschrijving 
van een onderzoekgebied in Nederland. In het laatste geral is begrijpelijkheid voor 

de niet-specialist een vereiste; bij een symposiumverslag is dat gewenst maar uiter- 
aard niet altijd mogelijk. U zult dan ook voorzover u geen theoretisch fysicus bent, 
waarschijnlijk een aantal artikelen in dit nummer erg moeilijk vrinden. Laat dit u 
echter niet weerhouden met lezen te beginnen: de artikelen zijn naar opklimmende 
moeilijkheid gerangschikt en ook de niet-specialist kan toch kernmismaken met een 
aantal boeiende ontwikkelingen op dit gebied van onderzoek. 


A. van der Woude 


NEEN EEN 





Theoretici in de fysica 


Toen mij gevraagd werd in dit nummer, 
dat zo’n fraaie presentatie van een aan- 
tal theoretische aspecten van de statis- 
tische fysica geeft, een gastcommentaar te 
schrijven, kwam als vanzelfsprekend de 
rol van de theoretici in de fysica bij mij 
naar voren — een probleem waarmee ik 
in de afgelopen vijftig jaar trouwens ook 
herhaaldelijk geconfronteerd werd. Ik 
zou daarom hierover bij deze gelegenheid 
iets willen zeggen. 


Het is thans alweer ruim honderd jaar ge- 
leden dat de eerste expliciet voor de be- 
oefening van de theoretische fysica be- 
stemde leerstoelen in ons land ontston- 
den. De eerste leerstoel in Nederland, 
expliciet voor de mathematische fysica, 
werd in Leiden in 1878 met de benoe- 
ming van Lorentz (1878-1912) ingesteld. 
Deze leerstoel was ontstaan door op- 
splitsing van de leerstoel van Rijke, die 
zelf in 1882 door Kamerlingh Onnes 
werd opgevolgd. In Amsterdam werd in 
1877 Van der Waals (1877-1908) de 
eerste hoogleraar in de natuurkunde, 
maar dit was geen specifieke leerstoel 
voor de mathematische fysica: hij do- 
ceerde alle vakken en pas in 1697 ont- 
stond een opsplitsing: Van der Waals 
kreeg de mathematische fysica in 

engere zin tot opdracht. Naast hem wer- 
den toen de experimentatoren Sissing en 
Zeeman benoemd. In Groningen was al 
in 1893 Windt benoemd tot lector in de 
mathematische fysica, naast de experi- 
mentator Haga. In 1903 werd het lec- 
toraat omgezet in een protessoraat voor 


Van der Waals jr. Bij diens vertrek naar 
Amsterdam in 1908 werd het wederom 
een lectoraat, achtereenvolgens bezet 
door Ornstein en Zernike. Pas in 1920 
werd met de benoeming van Zernike tot 
hoogleraar de onafhankelijke leerstoel 
in de mathematische fysica hersteld. In 
Utrecht werd in 1915 naast Julius een 
onafhankelijke leerstoel in de mathema- 
tische fysica ingesteld, waarop Ornstein 
werd benoemd. Ornstein nam echter de 
gehele directie van het Physisch Labo- 
ratorium — ook met het experimentele 
werk — op zich, zodat in 1926 opnieuw 
een afsplitsing plaatsvond met de be- 
noeming van Kramers (1926-1935). 

In de jaren twintig hadden dus alle uni- 
versiteiten een aparte leerstoel in de ma- 
thematische fysica verkregen. De TH in 
Delft volgde in 1939 met een aparte 
leerstoel die door Kronig werd bezet. 
Allèên in Leiden bestond een apart in- 
stituut voor de theoretische fysica, met 
de beroemde leeskamer Bosscha. 

Thans, een vijftigtal jaren later, heeft 
bijna iedere instelling een apart theore- 
tisch fysisch instituut, waarin gemiddeld 
een vijftal hoogleraren — inclusief de 
vroegere lectoren — op verschillende ge- 
bieden van de theoretische fysica werk- 
zaam is in samenwerking met een vaste 
groep van universitaire wetenschap- 
pelijke medewerkers van ongeveer ge- 





Figuur voorpagina: Fase diagram in de 
p-T-x-ruimte. 
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lijke omvang. Daarbij zijn aan ieder van 
de universiteiten nog eens gemiddeld 
een kleine tien universitaire en FOM- 
medewerkers werkzaam aan hun pro- 
motie, waardoor de totale omvang van 
de wetenschappelijke staf — hoog- 
leraren, vaste medewerkers en promo- 
vendi — vaak al tegen de twintig loopt? 
Ik wil hiermee nadrukkelijk niet zeg- 
gen, dat de relatieve omvang van de 
theoretische fysica binnen het geheel 
van de fysica is toegenomen: preciese 
gegevens om dit na te gaan heb ik niet 
tot mijn beschikking. Het aantal leden 
van de NNV nam in de laatste 50 jaren 
globaal ook met ruim een factor 10 toe, 
van ruim 200 in de jaren twintig tot ruim 
2500 thans. Overigens steeg de Neder- 
landse bevolking in die vijftig jaren met 
nog geen factor 2. Het aantal natuur- 
kundestudenten steeg in deze vijftig ja- 
ren ook met ongeveer een factor 10, 
zodat het er op lijkt dat, althans in ons 
land, de relatieve omvang van de theore- 
tische fysica in de laatste vijftig jaren niet 
zo veel veranderd is. 

Wel is door deze ontwikkelingen wat ver- 
loren gegaan: de wisselwerking tussen 
theorie en experiment, die je als theore- 
ticus op een laboratorium ervaart, is mij 
altijd als een goede ervaring bijgebleven: 
de stimulans voor de theorie van de 
stroom nieuwe experimentele gegevens, 
waarmee je op een goed lab continu ge- 
confronteerd wordt en omgekeerd: de 
bevruchtende werking, die een goede 
theoretische begeleiding op het experi- 
menteel onderzoek kan uitoefenen. Het 
is deze directe wisselwerking van theorie 
en experiment, die eigenlijk door de 
meer gescheiden institutionalisering van 
de theoretische fysica in de laatste vijf- 
tig jaren voor een deel verloren is ge- 
gaan. Daar staat anderzijds tegenover 
dat de bundeling van krachten van theo- 
retici, zoals die weerspiegeld wordt in 
goed bemande theoretische instituten, 
met zijn aanzienlijke versterking van 
onderlinge contacten, een niet te keren 
ontwikkeling is dat, als men afsluiting 
naar buiten voorkomt, de theoretische 
fysica ook zeker ten goede komt. 

Het is duidelijk dat de onafhankelijke en 
expliciete vormgeving van de theoretische 
fysica in afzonderlijke theoretische insti- 
tuten, het resultaat is van een ontwikke- 
ling van de fysica in breedte en diepte 
met steeds verdere specialisatie. Deze 
ontwikkeling leidt tot de vaak gehoorde 
klacht van de experimentatoren dat 
theoretici ‘hun eigen gang maar gaan en 
zich weinig bemoeien met dicht bij huis 
verkregen resultaten’. Als remedie hier- 
tegen zijn dan ook op verschillende labo- 
ratoria ‘huistheoretici’ werkzaam, die 
zich juist die interpreterende en stimu- 


lerende medewerking met het experi- 


mentele onderzoek op het lab tot taak 
gesteld hebben. Op de FOM-instituten 
zijn zelfs zeer omvangrijke theoretische 
groepen gehuisvest, die dezelfde doel- 
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stelling hebben en het onderzoek ook 
zeer doeltreffend steunen. Specialisatie 
is niet alleen in de samenleving, maar 


_ook in de wetenschap een harde reali- 


teit geworden, die zijn tol eist. Ook 
binnen de fysica wordt het meer en meer 
zo, dat onderzoekers uit verschillende 
gebieden elkaar vaak niet meer begrij- 
pen of verstaan. Het vindt zijn weerslag 
in de wetenschappelijke tijdschriften, 

die steeds minder ‘algemeen’ worden 


en zich geleidelijk steeds meer opdelen 


in secties van een meer gespecialiseerd 


karakter. Hetzelfde geldt ook voor de 


wetenschappelijke conferenties, waar- 
door men in feite nog alleen echte con- 
tacten onderhoudt met onderzoekers 
van ongeveer dezelfde interesse. Ik wil 
hier echter op het verschijnsel van de 
toenemende differentiatie naar vakge- 
bieden — hoe belangrijk die ook is — 
niet nader ingaan. 


Ia 
\ al 
__WISIECUS 


De opsplitsing tussen experimentele en 
theoretische fysica is namelijk eigenlijk 
toch van een andere aard. Het gaat hier- 
bij immers niet om een verdere differen- 
tiatie tussen verschillende vakgebieden 
van de fysica, zoals hoge-energiefysica, 
statistische fysica of vaste-stoffysica, maar 
om een zich geleidelijk steeds verder toe- 
spitsende taakverdeling binnen ieder van 
de vakgebieden. Hierbij groeien theoreti- 
ci en experimentatoren niet uit elkaar: 
zij publiceren vaak in dezelfde tijdschrif- 
ten en bezoeken dezelfde internationale 
conferenties, waar confrontatie van ex- 
perimenteel onderzoek en theoretische 
ontwikkelingen aan de orde is. Het is in 
dit verband ook van belang er op te wij- 
zen dat in FOM de theoretische werk- 
groepen zich niet tot een aparte werk- 
gemeenschap hebben aaneengesloten, 
maar een deel zijn gebleven van de werk- 
gemeenschap waar ze qua vakgebied 
thuishoren. 

Ik zou de relatie tussen theorie en expe- 
riment nog wat willen verduidelijken aan 
de hand van enkele beschouwingen van 
de fysicus en wetenschapshistoricus Tho- 







mas Kuhn? die teruggaan naar de tijd dat 


de natuurkunde zelf als erkende discipli- 


ne tussen de wetenschappen ontstond. 


Feitelijk lijkt de fysica als erkende disci- 
pline van onderzoek tussen de verschil- 
lende wetenschappen nog maar een twee 


eeuwen oud te zijn. Nog in Newtons tijd 
was daar enerzijds de ‘natuurfilosofie’, 
die voornamelijk de wiskunde, de sterren- 
kunde, de statica en bewegingsleer en de 
geometrische optica omvatte — een 
voortzetting van antieke tradities, maar 
in een ingrijpend methodisch vernieuw- 
ingsproces betrokken — en anderzijds de 
zich ontwikkelende nieuwe empirische 
‘baconiaanse’ wetenschap, die zich met 
scheikunde, warmte, elektriciteit en 
magnetisme bezig hield, en die vaak — 
uit angst voor metafysische en aprio- 
ristische elementen — de wiskunde op 
enige afstand hielden. Pas bij de over- 
gang van de 18° naar de 19° eeuw vond 
— hoofdzakelijk in Frankrijk — een in- 
grijpend ‘mathematiseringsproces’ 
plaats* waarbij in betrekkelijk korte tijd 
empirische wetenschapsgebieden zich 
door gebruik van methoden van de wis- 
kunde, ontwikkelden tot volwassen ge- 
bieden van wetenschappelijk onderzoek. 
Men denke daarbij aan Fourier, Laplace en 
Carnot voor de warmtetheorie, Poisson, 
Biot en Ampère voor elektriciteit en mag- 
netisme en Fresnel voor de optica. Hier- 
door ontstond, tezamen met de alreeds 
in een eerder stadium tot ontwikkeling 
gekomen mechanica (Fermat, Newton, 
Lagrange), elasticiteitstheorie en hydro- 
dynamica (de Bernouilli’s Euler, enz), 
een volwaardige nieuwe discipline — de 
fysica — tussen de wiskunde enerzijds en 
de scheikundige wetenschappen ander- 
zijds. In deze nieuwe discipline werden 
belangrijke wiskundige tradities uit de 
natuurfilosofie en experimentele tradi- 
ties uit de empirische (“baconiaanse’) 
wetenschappen geïntegreerd. Het leek 
dat hierdoor een langdurige antithese 
tussen de meer wiskundige en de empi- 
rische wetenschappen werd opgelost, 
maar — zoals Kuhn terecht opmerkt — 
het zou beter zijn te zeggen dat deze 
‘tegenstelling’ in de kiem bleef bestaan, 
maar nu binnen de fysica zelf, waar het 
een bron van zowel individuele als | 
professionele spanning bleef, die ook 
thans nog wel eens doorklinkt. 

Na de geweldige bijdrage, die de Franse 
wetenschap in het begin van de 19° eeuw 
tot dit mathematiseringsproces en daar- 
bij tot het ontstaan van de fysica als vak- 
gebied tussen wiskunde en scheikunde 
leverde, werd het vaandel van de fysica 
in belangrijke mate overgenomen door 
Engeland en Duitsland. Ik wil op enkele 
recente studies® over de ontwikkeling 
in het Duitsland van begin 19° eeuw 
nog wat nauwkeuriger ingaan. 

Door Cormmach® wordt Franz. E. Neu- 
mann in Königsberg in de jaren dertig 
en veertig gezien als de meest prominen- 
te voorstander van de integratie van de 
kwantitatieve methoden, die aan de 
wiskunde ontleend zijn, met de be- 
staande empirische methoden, die vaak 
nog een zeer kwalitatief karakter had- 
den. Neumann sprak hier van de ‘mathe- 


matische fysica’ als synthese tussen het 
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rationele denken van de wiskunde en 

het kwalitatieve, “anschauliche’, em- 
pirische onderzoek. 

Ideaal daarbij is de vereniging van 
wiskunde en experimentele kwaliteiten 
in één persoon, zoals ook bij Helmholtz 
en Kirchhof, beiden zeer sterk door 
Neumann beïnvloed, het geval was. Toch 
was dit een ideaal dat allêèn door de heel 
groten kon worden waargemaakt. Al 
spoedig zou er opnieuw een splitsing, 
maar nu in de vorm van een taakverde- 
ling tussen experimentele en theoretische 
fysica ontstaan: al in 1849 splitste 
Weber het fysisch instituut in Göttingen 
in een theoretische en een experimentele 
afdeling. Kohlrausch werd hier in 1866 
de grote pur sang experimentator, wiens 
beroemde Lehrbuch der Praktischen 
Physik (1870) voor vele generaties een 
eerste handleiding in nauwkeurig meten 
betekende en een solide basis voor de 
experimentele fysica gaf. Anderzijds 
kwamen ook aparte leerstoelen voor 

de theoretische fysica van de grond: van 
Boltzmann in 1869 in Graz en voor 


Kirchhof in 1875 in Berlijn. In 1889 werd 
Max Planck in Berlijn als opvolger van 
Kirchhof de leider van het nieuwe 
“Institut für Theoretische Physik’. 

Zo was korte tijd na de synthese 

tussen ratio en empirie in het begin van 
de 19° eeuw de verdeling tussen theorie 
en experiment weer een feit geworden, 
maar nu toch op geheel andere wijze: als 
specialisatie in twee aspecten van dezelf- 
de wetenschap, die elkaar nodig hebben 
en wederzijds aanvullen. En zo is de 
fysica ook in de laatste honderd jaren 
gebleven: een uniek voorbeeld van een 
kwantitatieve wetenschap met een spe- 
cialisatie in de breedte: de verschillende 
vakgebieden die tot de fysica gerekend 
werden en in de diepte: waarbij de theo- 
retische en experimentele research el- 
kaar aanvullen. Ik heb in dit commen- 
taar trachten te betogen, dat theoreti- 

ci en experimentatoren bij elkaar horen 
— als het kan in één laboratorium. De 
communicatie is thans echter goed 
genoeg om deze eenheid van plaats niet 
meer als een voorwaarde voor goed con- 


tact tussen theorie en experiment te 
vereisen: tijdschriften, letters en confe- 
renties zijn de media waar men elkaar 
ontmoeten kan. Ik stel me evenwel 
niet voor dat hiermee de discussie, die 
het uitgangspunt van dit commentaar 
vormde, de wereld uit is. 

J. de Boer 


Referenties 

1. The development of science in the Nether- 
lands during the last half century. Leiden, Van 
Doesburgh (1930). 

2. FOM-rapport: De beoefening van de theo- 
retische natuurkunde in Nederland. FOM 
32709/2 (1973). 

3. Thomas Kuhn: Zhe function of measure- 
ment in modern science ISIS 52 (1961) p. 
161/193; The Essential Tension (1977) p. 178/ 
224. 

4. Zie verschillende bijdragen in: R. 
MeCormmach (Ed.) Hist. Stud. in the Phys. 
Sc., deel 4 (1974) en 8 (1977). 

5. R. McCormmach: Editors Forword in Hist. 
Stud. Phys. Sc. 3 (1971). Zie ook: K.L. Cane- 
va: From Galvanism to Electrodynamics: The 
Transformation of German Physics and Its 
Social Context in R. McCormmach (Ed) Hist. 
Stud. Phys. Sc. 9 (1978) p. 63/160. 


RN RED TL TEL NDB RE EN AR nnee 





Enige persoonlijke reminiscenties 


aan Jan de Boer 


E.G.D. Cohen 


Jan de Boers werk is veel te uitgebreid en gevarieerd om in één artikel te bespreken. Ik 
zal mij daarom beperken tot die aspecten waar ik zelf min of meer:bij betrokken ben 
geweest. Twee van de belangrijkste gebieden waar De Boer op gewerkt heeft zijn: 

1. De invloed van de quantummechanica op de evenwichts- en niet-evenwichtseigen- 
schappen van verdunde gassen. Hiermee is zijn werk over het quantummechanische 
principe van overeenstemmende toestanden direct verbonden. 2. De dichtheids- of 
viriaalontwikkeling van de hogere verdelingsfuncties, speciaal de paar-verdelingsfunc- 
tie, in een gas in evenwicht. Ik zal zijn publikaties op deze twee gebieden in enig detail 
bespreken! , de rest zal meer summier behandeld moeten worden. 


Ik zal beginnen met De Boers werk in de 
quantummechanica. Laat ik eerst iets zeg- 
gen over hoe de toestand was toen De 
Boer zijn wetenschappelijke carriëre be- 
gon in het midden van de jaren dertig. 
Merkwaardigerwijze werden de quantum- 
mechanica en de quantumstatistiek van 
de jaren twintig eerder toegepast op de 
transportcoëfficiënten van een verdund 
gas — dit is een gas van dermate lage dicht- 
heid dat alleen binaire botsingen optre- 
den — dan op de toestandsvergelijking. 


In 1933 en 1934 losten Vehling en Uhlen- 
beck de quantummechanische transport- 
vergelijking op voor een verdund gas, een 
vergelijking die al was voorgesteld in 1928 
door Nordheim!®. Zij verkregen uitdruk- 
kingen voor de viscositeit n en het warmte- 
geleidingsvermogen À, die afhingen van de 
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intermoleculaire potentiaal. Deze quan- 
tummechanische transportvergelijking 
verschilde van de klassieke Boltzmann- 
vergelijking in de Stoszzahl Ansatz, dat is 
de veronderstelling over het gemiddelde 
aantal binaire botsingen in een gas. Deze 
veronderstelling was op twee punten ge- 
wijzigd om met de quantummechanica in 
overeenstemming te zijn. Ten eerste was 
de klassieke botsingsdoorsnede voor twee 
moleculen vervangen door de quantum- 
mechanische botsingsdoorsnede, die het 
golfkarakter van de botsende moleculen 
in rekening brengt. Dit zogenaamde dif- 
fractie-effect wordt bepaald door 


\/o waar \= hlOnmkaT)® 


evenredig is aan de De Broglie-golflengte 
van een molecule dat zich beweegt met 
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de gemiddelde impuls niv « (mk T)” en 
waar 0 de afmeting (diameter) van een 
molecule is. Hierin is 7m de massa en v de 
gemiddelde snelheid van een molecule, Kg 
is de constante van Boltzmann en 7’ de 
absolute temperatuur. Aangezien A a 7 
zullen diffractie-effecten belangrijk wor- 
den bij lage temperaturen. Verder is er 
een tweede quantummechanisch effect op 
de botsingsdoorsnede vanwege het feit 
dat gesymmetriseerde golffuncties moe- 
ten worden gebruikt voor de twee botsen- 
de moleculen: symmetrisch of anti-sym- 
metrisch in de coördinaten van de mole- 
culen naar gelang het gas Bose-Einstein 
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(B.E.) of Fermi-Dirac (F.D) statistiek 
volgt. | 

Afgezien van dit twee-deeltjes statistiek- 
effect, is er ook een veel-deeltjes statistiek- 
effect. Dit effect houdt rekening met de 
bezetting door andere deeltjes van de toe- 
stand waarin de twee botsende moleculen 
zullen komen na de botsing. Indien name- 
lijk de eindtoestand al bezet is, zal de 
waarschijnlijkheid dat de twee botsende 
deeltjes in deze toestand terecht komen 
toe(af)nemen als het gas B.E. (F.D) sta- 
tistiek volgt. Dit statistiekeffect wordt be- 
paald door de parameter n\° waar n 

de deeltjesdichtheid is. Deze quantum- 
statistiekeffecten zullen dus belangrijk 
worden bij lage temperaturen en relatief 
hoge dichtheden, als 


n\? = (N/a)® > 1, waaraan”! 


de gemiddelde afstand tussen twee gas- 
moleculen is. Als \ > a is zullen de golf- 
pakketjes, die naburige moleculen be- 
schrijven, elkaar overlappen en aan de 
quantumstatistiek zal streng de hand ge- 
houden worden. 


De uitdrukking die Vehling en Uhlenbeck 
voor de viscositeit vonden was, afgezien 
van een voor ons onbelangrijke factor: 


(1) 


Het accent bij nijgass betekent dat bij de 
berekeningen van nyjass de klassieke uit- 
drukking voor n met een quantummecha- 
nische botsingsdoorsnede gebruikt moet 
worden. Deze quantummechanische bot- 
singsdoorsnede kan worden uitgedrukt 

in de quantummechanische fasen n,(k), 
die de botsing karakteriseren (voor een 
gegeven intermoleculair potentiaal) van 
twee deeltjes die een draai-impuls 

[LU + 1) en een kinetische energie van 
relatieve beweging n° k*/m hebben. Hier 
is k het golfgetal dat met de relatieve be- 
weging geassocieerd is. Deze fasen kunnen 
worden gevonden door oplossing van het 
quantummechanische twee-deeltjes pro- 
bleem. Het +(—) teken in (1) is van toe- 
passing als het gas B.E. (F.D) statistiek 
volgt. In het midden van de jaren dertig 
waren de fasen en nyiass(7') alleen nog 
maar berekend voor het harde bollenmo- 
del van moleculen en was niets bekend 
over meer realistische moleculaire model- 
len. 


nn, T) = Nidass (CT) [1 + nX' Ön (TD). 


Voor ik doorga met mijn bespreking van 
de berekening van de transportcoëfficiën- 
ten bij lage temperaturen, waar de boven- 
genoemde quantumeffecten een rol zullen 
spelen, wil ik eerst iets zeggen over de in- 
vloed van de quantummechanica op de 
toestandsvergelijking van een gas!2. De 
theorie hiervoor is ontwikkeld in 1936- 
1937 door Uhlenbeck en Beth. Zij be- 
schouwden de tweede viriaalcoëfficiënt, 
dat is de coëfficiënt die de afwijkingen 
van de ideale gaswet bepaalt, tengevolge 
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van de wisselwerkingen tussen twee deel- 
tjes. In Kamerlingh Onnes’ viriaalontwik- 
keling voor de druk: 


pV B(T) 
nn | 


C(T) 
RT y + 


vi + .…., (2) 


is de tweede viriaalcoëfficiënt B(T) ge- 
geven in de klassieke statistische mecha- 
nica door de uitdrukking: 


B(T)= 
nm 2aN f [e-VOUKBT _ 1]r?dr. (3) 


Hier is N het getal van Avogadro; V het 
volume van het gas en R de gasconstante. 
Het intermoleculaire potentiaalveld V(r) 
is hier bolvormig verondersteld. Uhlen- 
beck en Beth generaliseerden deze klas- 
sieke uitdrukking tot de quantummecha- 
nica en vonden dat B(T) uitgedrukt kon 
worden in de verstrooiingsfasen n‚(k) en 
in de discrete energieniveaus van het po- 
tentiaalveld V(r): 


B(T) = Bia(T) -NV2 N° z (2! +1) x 


5 Jos # ke JmkepT} (amyll dk + 


+ XY er Fdisr/kB7 | (4) 
discr 
Hier is Bij (T') = + N\°/25/? de tweede 
viriaalcoëfficiënt voor een ideaal BE. (—) 
of F.D. (+) gas; de som over / gaat voor 
B.E. (F.D.) statistiek over alleen even 
(oneven) waarden van /; de integratie over 
k is over alle waarden van de energie van 
de relatieve beweging n° k?/m en de som 
B over alle discrete energieniveaus van 
ISCT 
V(r), indien aanwezig. Beth en Uhlenbeck 
in 1937 en later Kahn en Uhlenbeck in 
1938 gaven wel een uitvoerige discussie 
van (4), in het bijzonder over de afhanke- 
lijkheid van het gedrag van B(T) bij zeer 
lage temperaturen van de fase no(k), maar 
zij konden B(T) niet uitrekenen. Uhlen- 
beck dacht in feite dat dit onmogelijk was 
voor een enigszins realistische intermole- 
culaire potentiaal. 


Nu hadden De Boer en Michels in 1938 


zo’n potentiaal afgeleid door de experi- 


mentele temperatuur afhankelijkheid van 
de tweede viriaalcoëfficiënt van helium 
boven 50 K met die van (3) te vergelijken. 
Zij namen aan dat V(r) gekarakteriseerd 
kon worden door slechts twee parameters: 
een wisselwerkingsenergie en een lengte. 
Hoewel dit al eerder door Lennard-Jones 
was voorgesteld, waren het De Boeren 
Michels die de intermoleculaire wisselwer- 
king tussen twee moleculen voor het eerst 
in de nu standaardvorm schreven (verge- 
lijk figuur la): 


V (rx) 











Fig. 1. (a) Een aantal intermoleculaire 
potentialen (S, A, C,‚ BM) analoog aan (6) 
met dezelfde o maar verschillende e; BM 
is de De Boer-Michels-potentiaal; r* = r/o. 

(b) B*T*/2 (B*=B/No’, T*= 
kT/e) als functie van T* voor de poten- 
tialen S, A, C, BM, x is experiment. 


Vr) = 4e [(o/r)*? — (o/)®1. (5) 
In (5) is e de diepte van de potentiaal, dat 
ise= V(r in) waar Fmin volgt uit (dV(r)/ 
dr) oi = Cen o is de ‘diameter’ van het 
molecule dit is V(o) = 0. Voor helium 
vonden De Boer en Michels e = 14,03 x 
107 ergen o= 2,56 Â en voor deze waar- 
den van e en o wordt (5) de De Boer- 
Michels(BM)potentiaal genoemd. 
Zoals al eerder gezegd, leek het onmoge- 
lijk om het quantummechanische ver- 
strooiingsprobleem, dat is de Schrödinger- 
vergelijking van twee deeltjes, met een 
wisselwerkingspotentiaal (5) op te lossen 
en daaruit de fasen n;(k) te berekenen. 
Toch is dit net wat De Boer deed, onge- 
twijfeld gestimuleerd door zijn regelmatige 
bezoeken aan Utrecht, waar Uhlenbeck 
hoogleraar was en Kahn zijn student. Met 
wat bij hedendaagse standaarden uiterst 
primitieve rekenfaciliteiten genoemd 
moeten worden, loste De Boer Schrödin- 
gers vergelijking voor de BM-potentiaal 
numeriek op en bepaalde negen en twintig 
fasen n,(k): van /= O tot 6 en voor ko van 
0,25 tot 4,00. Dit maakte een berekening 
van B(T) alleen beneden temperaturen 
van 5 K mogelijk. Er is geen discreet 
niveau voor de BM-potentiaal, maar er is 
een virtueel niveau vlakbij nul, dat het 
gedrag van B(T) bij zeer lage temperatu- 
ren bepaalt via het gedrag van no (Kk) bij 
kleine k. De vergelijking met het experi- 
ment is gegeven in figuur 1b (de BM-krom- 
me). Om de invloed van de quantum- 
mechanica te appreciëren moet men be- 
denken dat de klassieke waarde van 
B*T*3/2 bij T*= 0,3 bijvoorbeeld, gelijk 
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aan —9 is! De quantummechanische theo- 
rie geeft de juiste orde van grootte en het 
juiste globale gedrag. Toch zullen de grote 
afhankelijkheid van V(r) en daarom de 
essentiële beperkingen van de theorie, ge- 
zien de onbekendheid met de ware poten- 
tiaal, duidelijk zijn. Het is een bewijs van 
De Boers gezond oordeel dat hij hier altijd 
de nadruk op heeft gelegd en nooit ge- 
probeerd heeft ‘perfecte’ overeenstem- 
ming met het experiment te krijgen. Toch 
liet dit werk duidelijk zien dat voor het 
juiste gas — helium in dit geval — quan- 
tumeffecten een belangrijke rol spelen 

bij lage temperaturen, zelfs voor zulke 
_fklassieke’ macroscopische eigenschappen 
als de toestandsvergelijking. 


Een van de meest belangrijke ontdekkin- 
gen van De Boer was een preciese formu- 
lering van wat het juiste gas’ betekent. Ik 
zinspeel hier op de invoering van de quan- 
tummechanische parameter — of zoals hij 
nu algemeen wordt genoemd, de De Boer- 
parameter — en de quantummechanische 
wet van overeenstemmende toestanden. 
De klassieke wet van overeenstemmende 
toestanden van Van der Waals leerde dat 
als men de druk p, het volume V en de 
temperatuur 7 van een gas dimensieloos 
maakte door ze te delen door hun waar- 
den p‚., Ve en Te bij het kritieke punt, 
een universele toestandsvergelijking, gel- 
dig voor alle gassen, uit de Van der Waal- 
se toestandsvergelijking volgt. 


De Boer merkte nu allereerst op dat voor 
een intermoleculaire potentiaal van de 
vorm (5), die van slechts twee parameters 
e en o afhangt, óók een universele toe- 
| standsvergelijking geldt: 

p*=pr(V* TE), (6) 
waar p*= po’ Je, V*= V/No’ en T*= 
kT/e. Dit volgt uit het feit dat in het alge- 
meen de toestandsvergelijking van de 
vorm: p =P(kKT,N, V‚,m,e,0) moet zijn. 
Construeert men nu uit e en o een groot- 
heid van de dimensie van een druk: e/o° 
dan kan de dimensieloze druk p* alleen 
nog maar een functie zijn van de dimensie- 
loze grootheden V* en 7*. In de quan- 
tummechanica echter treedt de Planckse 
constante / op. Gezien de dimensie van 

h (energie x tijd), kan men nu een nieuwe 
dimensieloze constante construeren, die 
dus ook in de toestandsvergelijking kan 
voorkomen: A* = h/o(me)”. De toe- 
standsvergelijking kan daarom in de quan- 
tummechanica de vorm hebben: 
pr pF TAS), (7) 
en is dus niet meer universeel aangezien 
hij (parametrisch) afhangt van A*. Hoe 
groter A*,hoe lichter en kleiner het mole- 
cule en hoe zwakker de intermoleculaire 
wisselwerking, des te groter de afwijkin- 
gen van het klassieke gedrag zullen zijn. 
(7) is een voorbeeld van het quantum- 
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mechanische principe van de correspon- 
derende toestanden. 

In 1948 pasten De Boer en Lunbeck een 
gegeneraliseerde vorm van dit quantum- 
mechanische principe van corresponderen- 
de toestanden toe op een groot aantal 
eigenschappen van een aantal stoffen. 
Anneke Sengers zal dit werk in meer detail 
bespreken, dus ik zal mij beperken tot 
één heel frappant voorbeeld. Door dit 
principe toe te passen op een aantal stof- 
fen met toenemende waarden van A* 
(Xe: A* = 0,064; Kr: A* = 0,102; Ar: 
A= 0187: Ne: A*= 0,501: Hy: A" = 
1,73; He“: A*= 2,64; He®: A* = 3,05) 
konden De Boer en Lunbeck door extra- 
polatie de kritieke temperatuur en druk, 
de inwendige energie en dichtheid bij 0 K 
en de temperatuur afhankelijkheid van de 
dampdruk van He® voorspellen, een stof 
waarvan de thermodynamische eigen- 
schappen toen nog niet bekend waren. 
De overeenstemming met later gedane ex- 
perimenten was opvallend en de mensen 
in Amerika vroegen ‘hoe krijgt hij het 
voor elkaar? 


Ik kom nu terug op de berekening van de 
transport coëfficiënten van een gas met 
de intermoleculaire potentiaal (5) bij lage 
temperaturen. In 1943 berekende De Boer 
de viscositeit n en het warmtegeleidings- 
vermogen À, van gasvormig helium bene- 
den 5 K met behulp van dezelfde fasen 
die hij had gebruikt in 1938 voor de 
tweede viriaal coëfficiënt. 

De overeenstemming met het experiment 
was heel bevredigend. De berekeningen 
lieten zien dat de statistiekeffecten , en 
analoog ô, erg klein waren behalve voor 
temperaturen beneden 1 K. Men kan zich 
afvragen of bij de tegenwoordige stand 
van experimenteren deze interessante veel- 
deeltjes effecten in n en À niet gemeten 
zouden kunnen worden. De berekening 
van de tweede viriaal coëfficiënt en van 

n en \ van andere lichte gassen, zoals He® 
en para- en ortho-waterstof en hun meng- 
sels zijn uitgevoerd door veel studenten 
van De Boer in de jaren veertig en vijftig. 
Vele interessante en vaak subtiele quan- 
tumeffecten ten gevolge van diffractie en 
(twee-deeltjes) symmetrie-effecten werden 
gevonden, die essentieel waren voor het 
begrijpen van de evenwichts- en transport- 
eigenschappen van deze lichte gassen bij 


lage temperaturen. 


Ik kom nu tot het tweede gedeelte van 
mijn artikel, namelijk het werk van De 
Boer aan de paar-verdelingsfunctie van 
een gas in evenwicht!?. De paar-verde- 
lingsfunctie n9 (r, ‚72 ) geeft aan het ge- 
middelde aantal moleculeparen dat zich 
in een systeem op de plaatsen r‚ en rz 
bevindt. 


In evenwicht hangt n, alleen af van de 
afstand rj, = |r, — rs |, zodat dan 

nz (ri, F2) = 1 (712). Voor een ideaal gas 
geldt 7, (r) = n°, maar in een echt gas zal 


3 
gl), (a) 
Ë u 
05 2 3 
— r* 
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Fig. 2. (a) g(r) bepaald door fourier- 
inversie uit röntgenexperimenten van 
kwik (r*=r/o, o=2,98Á). 

(b) g(r) berekend uit vergelijking 
(11), met de De Boer-Michels-potentiaal 


(r*=r/o). 


nz (r) van de intermoleculaire potentiaal 
V(r) afhangen en in het algemeen zal 
daarom: 

n(r)=n°g(r), (8) 
zijn, waar g(r) de radiële verdelingsfunctie 
is. Voor grote waarden van 7, als de twee- 
deeltjes ver van elkaar zijn en hun wissel- 
werking verwaarloosd kan worden, nadert 
z(r) tot 1. De g(r)-functie is ingevoerd 
door Ornstein en Zernike in 1914 in een 
theorie die poogde de kritische opales- 
centie, dat is de abnormaal grote lichtver- 
strooiing bij het kritieke punt, te verkla- 
ren. Ornstein en Zernike lieten zien dat 
voor grote waarden van r geldt: 

gr) lae"/r, (9) 
waarin K omgekeerd evenredig aan de 
compressibiliteit van het gas is. Dus, bij 
nadering van het kritieke punt, als Kk — O0 
krijgt g(r) — 1 een lange dracht, aangezien 
dan g(r) — 1 > 1I/r. Later, in 1927, heb- 
ben Zernike en Prins deze theorie toege- 
past op de röntgenverstrooiing van vloei- 
stoffen. In 1932 lieten Debije en anderen 
zien hoe men door fourierinversie de g(r) 
kon bepalen uit röntgenspectra en g(r)- 
krommen, zoals die in figuur 2a voor 
vloeibaar kwik werden verkregen. Hoe- 
wel een aantal primitieve moleculaire be- 
rekeningen van g(r) gedaan waren, werd 
de eerste theoretische berekening van (7) 
op basis van de statistische mechanica uit- 
gevoerd door Kirkwood in 1935. Kirk- 
wood vond de volgende benaderde uit- 
drukking van de radiële verdelingsfunctie: 
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waar B = I/kBT en rj de afstand tussen de 
moleculen i en j is. Hogere benaderingen 
konden slechts met grote moeite uitgere- 
kend worden. Aan de andere kant hadden 
Uhlenbeck en Kahn in 1938 laten zien, 
hoe men op systematische wijze een dicht- 
heidsontwikkeling van de thermodynami- 
sche eigenschappen van een gas kon ver- 
krijgen op grond van een cluster ontwik- 
kelingsmethode van Ursell. De Boer gene- 
raliseerde deze methode in 1939 en ver- 
kreeg zodoende een dichtheidsontwikke- 
ling van de paar verdelingsfunctie. Hij 
vond: 


gr) = e-BV2) E + 


+ nfdr3 (ooren _ i (even en is 


+ n° mj 


In (11) heb ik alleen de eerste twee termen 
opgeschreven: De Boer gaf nog verschil- 
lende hogere-ordetermen. Voor lage dicht- 
heden reduceert (10) zich tot (11). Ver- 
volgens berekende De Boer met de BM- 
potentiaal numeriek g(r) en verkreeg zo 
een kromme (vergelijk figuur 2b) die een 
analoge structuur had als de door fourier- 
inversie van röntgenspectra gevonden 
kromme van figuur 2a. De Boer gaf een 
hele eenvoudige fysische afleiding van (10) 
en een hele interessante verklaring van het 
oscillerend gedrag van g(r) die ik, helaas, 
vanwege plaatsgebrek hier niet kan repro- 
duceren. Men kan ze echter in zijn disser- 
tatie en in het klassieke overzichtsartikel 
over moleculaire verdelingsfuncties in de 
“Reports and Progress in Physics” van 
1949 vinden. Dit artikel bevat veel mate- 
riaal dat oorspronkelijk in 1940 in zijn 
bijzonder fraai geschreven dissertatie was 
verschenen, maar dat, vanwege de oorlog, 
niet internationaal bekend was geworden. 
De Boers dichtheidsontwikkeling van de 
z(r) was gedaan vóór het analoge werk 
van Mayer en Montroll in 1941 en hij 
gebruikte deze ontwikkeling vóór hen en 
vóór het veel uitgebreidere latere werk 
van Bird, Spotz en Hirschfelder in 1950, 
om voor het eerst de derde viriaal coêffi- 
ciënt C(T) voor een realistische inter- 
moleculaire potentiaal te berekenen. 


(1) 


Ik wil nu iets zeggen over de buitengewone 
wetenschappelijke activiteit die De Boer 
gedurende de oorlogsjaren ontplooide. 
Bijzonder interessant is allereerst zijn 
werk over de transporteigenschappen van 
verdunde meeratomige gassen. Hier gene- 
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raliseerde hij de gewone Boltzmann-ver- 
gelijking voor één-atomige gassen tot 
meeratomige gassen, waarbij de inwendige 
vrijheidsgraden van de moleculen quan- 
tummechanisch behandeld moesten wor- 
den. Dit werk werd pas veel later, in 1964, 
gepubliceerd tesamen met het werk van 
Wang Chang en Uhlenbeck, die onafhan- 
kelijk van De Boer tot analoge resultaten 
waren gekomen. De publikatie geschiedde 
in volume 2 van de serie “Studies in Statis- 
tical Mechanics”, waarvan De Boer en 
Uhlenbeck de oprichters waren. 


Het is ook, geloof ik, de moeite waard te 
vermelden dat de Boer in deze tijd afge- 
leid heeft wat nu de B-B-G-K-Y-hiërarchie 
genoemd wordt naar Bogolubov, Born en 
Green, Kirkwood en Yvon; echter — 
typisch voor hem — heeft hij dit nooit 
gepubliceerd, omdat hij vond dat zolang 
men er niets mee deed, publiceren niet de 
moeite waard was. Verder rekende hij in 
een zeer redelijke benadering het niet- 
bolvormige potentiaalveld tussen twee 
waterstofmoleculen uit. Dit, of een ana- 
loog potentiaalveld, zou van nut kunnen 
zijn om de transportcoëfficiënten van niet- 
bolvormige moleculen te begrijpen, spe- 
ciaal in een uitwendig magnetisch of elek- 
trisch veld. Ik geloof dat als zulke poten- 
tialen werden gebruikt om de evenwichts- 
(B(T))- en transporteigenschappen (n, \) 
van verdunde gassen bestaande uit niet- 
bolvormige moleculen uit te rekenen, 
men veel zou kunnen leren over de niet- 
bolvormigheid van zulke moleculen, bij- 
voorbeeld door de theoretische resultaten 
te vergelijken met de experimenten van 
Beenakker en Knaap? . Dit zou dan uit- 
eindelijk kunnen leiden tot een begrijpen 
van de eigenschappen van niet-bolvormige 
moleculen, op basis van de intermolecu- 
laire potentiaal dat nu alleen nog bestaat 
voor bolvormige moleculen. 

Ik kom nu tot enige onderzoekingen die 
De Boer heeft gedaan met sommige van 
zijn studenten in de jaren vijftig. 


Een onderzoek, dat mij persoonlijk na 
aan het hart ligt, kwam voort uit De 
Boers idee om een clusterontwikkeling te 
gebruiken voor een theorie van vloeistof- 
fen. Hij paste dit idee toe op het celmodel 
van Lennard-Jones en Devonshire voor 
vloeistoffen. In dit model wordt veronder- 
steld dat elke molecule zich altijd dichtbij 
één roosterpunt van een virtueel rooster 
in de vloeistof bevindt. Elke molecule 
beweegt zich, onafhankelijk van de andere 
moleculen, om zijn roosterpunt in een cel 
ter grootte van het volume per molecule 
en in een potentiaalveld dat bepaald 
wordt door de (MN — 1) op hun rooster- 
plaatsen geplaatste andere moleculen. Het 
idee van de Boers celclustertheorie van de 
vloeistof is om dit ‘onafhankelijke mole- 
cule’ model als eerste benadering van een 
systematische theorie te beschouwen. Als 
eerste correctie wordt dan genomen de 


beweging van twee moleculen in een dub- 


belcel (of een cluster van twee cellen), in 
het potentiaalveld van de (VN — 2) op hun 
rooster geplaatste andere moleculen enz. 
Een motivering voor deze vloeistoftheorie 
was vloeibaar helium II, omdat de essen- 
tiële Bose-Einstein-statistiek in dit model 
ingevoerd kan worden door te eisen dat 
de golffunctie die de twee moleculen in 
een celcluster van twee cellen beschrijft, 
symmetrisch in de coördinaten van deze 
moleculen moet zijn. Merkwaardiger- 
wijze is deze celcluster-theorie van de 
vloeistof, die het onderwerp van mijn 
dissertatie was®, uiteindelijk noch op 
vloeibaar helium noch op gewone vloei- 
stoffen toegepast, maar op vaste stoffen. 
Verleden jaar nog hebben Barker, Hender- 
son en Abraham* een versie van deze 
theorie gebruikt om de vaste-stoflijn in 
het fasediagram van een twee-dimensio- 
nale vloeistof, met een intermoleculaire 
potentiaal van de vorm (5) te berekenen, 
ten einde aan te tonen dat, in tegenstelling 
tot recente theoretische voorspellingen?” , 
smelten in twee dimensies een gewone 
eerste orde fase-overgang is. 


Andereinteressante onderzoekingen waren 
die met Hijmans over een orde-wanorde 
theorie van fase-overgangen en met Van 
Leeuwen en Groeneveld?’ over een nieuwe 
methode om de radiële verdelingsfunctie 
uit te rekenen. In dit laatste werk werd, 
wat nu de “hyper-netted chain” benade- 
ring wordt genoemd, ontwikkeld op een 
bijzonder fraaie en overzichtelijke manier. 
Deze methode, in combinatie met andere 
analoge benaderingsmethoden, hebben 
met behulp van moderne rekenmachines. 
de berekening van de radiële verdelings- 
functies in verdichte vloeistoffen erg ver- 
beterd. 


Voor ik eindig met een discussie over een 
eigen onderzoek dat als een voortzetting 
of generalisatie beschouwd kan worden _ 
van De Boers werk over de g(r)-functie in 
evenwicht, wil ik eerst een paar persoon- 
lijke woorden zeggen over De Boer en zijn 
werk zoals dit naar voren komt uit een 
studie van zijn publikaties en uit de erva- 
ringen van een ex-student. 


Laat ik allereerst zeggen dat De Boer ‘een 
gelukkige hand’ heeft, wat geen toeval is 
aangezien hij “very good sense and very 
sound judgement” heeft. Hij verstaat de 
kunst om zinvolle vereenvoudigende be- 
naderingen te maken in moeilijke theo- 
retische berekeningen, die dan bijna altijd 
tot goede overeenstemming met het ex- 
periment leiden. Hij heeft een uitzonder- 
lijk ‘plastisch’ fysisch inzicht en — zoals 
duidelijk is uit zijn overzichtsartikelen — 
een heel duidelijke en gebalanceerde op- 
vatting over wat het verband tussen de 
theorie en het experiment moet zijn. In 
additie tot zijn dissertatie is er een aantal 
overzichtsartikelen die, volgens mij, nog 
steeds zeer de moeite waard zijn om te 
lezen, zoals het al eerder vermelde artikel 
127 


in “Reports on Progress in Physics” en 
twee hoofdstukken in het standaard werk 
“Molecular Theory of Gases and Liquids” 
van Hirschfelder, Curtiss en Bird: één 
over de quantumtheorie van de toestands- 
vergelijking (hoofdstuk 6) en één over de 
quantumtheorie van de transportver- 
schijnselen (hoofdstuk 10). 


Door van het werk van Uhlenbeck en 
Kirkwood te behouden wat waardevol 
was en te combineren met wat nuttig was, 
heeft De Boer belangrijke bijdragen ge- 
leverd tot de ontwikkeling van de statis- 
tische mechanica in de jaren dertig, veer- 
tig en vijftig. Een duidelijk teken dat hij 
alom beschouwd werd als een van de 
leidende figuren in de statistische me- 
chanica is dat Kirkwood, Mayer, Hirsch- 
felder en Onsager hem hun studenten na 
hun promotie toestuurden. Ook herinner 
ik mij nog goed hoe geïmponeerd ik was op 
het internationale symposium ““Transport 
Processes in Statistical Mechanics” dat 

in 1956 in Brussel gehouden werd. Alle 
grootmeesters van de statistische mecha- 
nica waren aanwezig: Onsager, De Mayers, 
Kirkwood, Uhlenbeck enz. Om een mij 
onbekende reden was De Boer niet aan- 
wezig. De eerste morgen presenteerden 
Kirkwood en anderen hun afleidingen van 
de Boltzmann-vergelijking uit de grond- 
slagen van de statistische mechanica. Ge- 
durende de discussie werd Kirkwood heel 
intensief door iemand ondervraagd over 
enige stappen in zijn afleiding. Kirkwood 
werd meer en meer in een hoek gedreven 
totdat hij eindelijk, min of meer ten einde 
raad, zei dat hij de gewraakte stappen in 
zijn afleiding gebaseerd had op resultaten 
van De Boer. Deze opmerking van Kirk- 
wood was duidelijk beslissend, want de 


ondervrager stopte onmiddellijk en de dis- 


cussie kwam tot een abrupt einde. 

Als student heeft mij misschien het meest 
beïnvloed, het voorbeeld dat De Boer 
stelde van wetenschappelijke eerlijkheid 
en integriteit. lets werkelijk begrijpen, 
jezelf niets wijs te maken en vooral om 
‘achter de formules’ te kijken om te zien 
wat zij fysisch betekenen en dit dan te 
kunnen zeggen ‘met je vingers’. Dat was 
de boodschap die hij zijn studenten gaf, 
in het bijzonder, tijdens de Donderdag- 
middag theoretische seminaria die, volgens 
mij, het belangrijkste onderdeel van de 
fysische opvoeding waren. Je niet laten 
imponeren door geleerde schijn maar, op 
misschien typisch nuchtere Hollandse 
wijze, te kijken wat de werkelijke inhoud 
is. Dit was voor een filosofisch georiën- 
teerde en humanistisch geïnteresseerde 
student een moeilijke maar ontzaglijk 
waardevolle opvoeding. 


Laat ik tot slot iets zeggen over enig recent 
werk dat ik met mijn student T.R. Kirk- 
patrick en met J.R. Dorfman® gedaan heb 
over de paar-verdelingsfunctie in een sta- 
tionaire toestand van een gas of vloeistof 
die niet in evenwicht zijn. Dit werk kan 
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gezien worden als een uitbreiding van De 
Boers werk over de paar-verdelingsfunctie 
in evenwicht tot het op één na eenvou- 
digste geval: de stationaire toestand. De 
stationaire toestand die wij in het bijzon- 
der bestudeerd hebben is degene, die zich 
voordoet in de zogenaamde Rayleigh- 
Bénard-cel. In dit geval wordt een vloei- 
stof (of gas) in een vat verwarmd, bijvoor- 
beeld van boven af, zodat er een constante 
warmtestroom door de vloeistof van boven 
naar beneden stroomt. Men kan zich nu 
afvragen wat de paar-verdelingsfunctie 

in een vloeistof in een dergelijke statio- 
naire niet-evenwichtstoestand is. Deze 
vraag is, zelfs voor een verdund gas, veel 
moeilijker te beantwoorden dan in even- 
wicht, aangezien er in niet-evenwicht geen 
nuttig algemeen beginpunt is, zoals de 
ensemble theorie van Gibbs in evenwicht. 
Zelfs de oude vraag van Ehrenfest? : ‘kan 
men bewijzen dat een irreversibel proces 
bij vaste uitwendige randvoorwaarden tot 
een unieke stationaire toestand zal naderen 
en kan men deze stationaire toestand op 
de een of andere wijze als de “relatief 
waarschijnlijkste”’ toestand karakteriseren 
en deze toestand dan kenmerken als het 
minimum van een of andere functie, die 
als generalisatie van de H-functie be- 
schouwd kan worden?’ is nog steeds niet 
beantwoord. Ik kan hier niet op de bere- 
keningen die wij gedaan hebben ingaan. 
Ik wil alleen zeggen dat wij de paar-verde- 
lingsfunctie voor een verdund gas op basis 
van kinetische theorie via een uitbreiding 
van de Boltzmann-vergelijking en voor 
een vloeistof via een uitbreiding van de 
hydrodynamische vergelijkingen hebben 
berekend. Soortgelijke berekeningen zijn 
uitgevoerd door Van der Zwan, Bedeaux 
en Mazur!® en anderen. De paar-verde- 
lingsfunctie nz (r, ‚rz } is nu niet meer, zo- 
als in evenwicht, bolsymmetrisch maar 
hangt af van r, en r‚. In het bijzonder 
hangt n, (ri ‚rz ) af van de oriëntatie van 
Fr, —rz ten opzichte van de richting van 
de warmtestroom. Deze asymmetrie van 
de paar-verdelingsfunctie zal te voorschijn 
komen in de lichtverstrooiing en heel re- 
cent zijn enige van de theoretische voor- 
spelde effecten, die afwijkingen van de 
verstrooiing in evenwicht betreffen, expe- 


rimenteel gezien!! . Een nog moeilijkere 


maar ook interessantere vraag is wat de 
paar-verdelingsfunctie — en de lichtver- 
strooing — is als de vloeistof van onderen 
af verwarmd wordt. In dat geval treedt 
vanwege de zwaartekracht, bij toenemend 
temperatuurverschil tussen beneden- en 
bovenwand van het vat, abrupt een insta- 
biliteit op, waar het karakter van de sta- 
tionaire toestand van de vloeistof veran- 
dert van een rustende tot een bewegende 
vloeistof, dat wil zeggen bij ‘te grote’ 
temperatuurverschillen treedt convectie 
in de vloeistof op. Dicht bij deze instabi- 
liteit vertoont de paar-verdelingsfunctie 
een analoge lange dracht als bij het kri- 
tieke punt in evenwicht (vergelijk (9). 
Deze lange dracht van de paar-verdelings- 


functie bij een instabiliteit zal tot een 
analogon van de kritieke opalescentie in 
evenwicht leiden, hoewel de fysische oor- 
zaak van de lange dracht hier een geheel 
andere js dan in evenwicht. Deze “insta- 
biliteits-opalescentie’ kan in principe aan 
de verstrooiing van de elektromagnetische 
golven gezien worden. De golflengte, no- 
dig om dit verschijnsel waar te nemen, ligt 
echter niet in het zichtbare gebied, zoals 
in evenwicht, maar in het microgolven 
gebied, zodat wellicht andere experimen- 
teertechnieken gebruikt zullen moeten 
worden. 


Dit is alles wat ik over dit recente werk 
wil zeggen. Het is een eerste poging om 
De Boers werk over de paar-verdelings- 
functie in evenwicht tot niet-evenwicht te 
generaliseren. De oplossing van dit pro- 
bleem bevindt zich nog in een rudimentair 
stadium en de resultaten, die tot nu toe 
verkregen zijn, zijn zeker niet zo klassiek 
als die De Boer verkreeg in ‘evenwicht. 
Denkende aan al het werk van De Boer 
dat ik hier genoemd heb, voel ik mij in 
eenzelfde situatie als Hajo Meijer vele jaren 
geleden was. Casimir had toen Hajo ge- 
vraagd om alle publikaties van De Boer te 
verzamelen en samen te vatten ter voor- 
bereiding, naar ik aanneem, op De Boers 
verkiezing tot lid van de Koninklijke 
Nederlandse Akademie van Wetenschap- 
pen. Toen Hajo zijn rapport aan Casimir 
gaf, vroeg deze: ‘ben je niet geïmpo- 
neerd?’, waarop Hajo — net als ik nu — 
antwoordde: “ja, heel erg’. 
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Recente ontwikkelingen in de theorie 
van de Boltzmann-vergelijking 


M.H. Ernst 


De Boltzmann-vergelijking, het fundament van de kinetische gastheorie, is ruim een 
eeuw oud, en veelvuldig door fysici en mathematen bestudeerd, toegepast en gegene- 
raliseerd. Zodra niet-lineaire effecten (schokgolven, relaxatieverschijnselen van hoog- 
energetische deeltjes, en dergelijke) een rol gaan spelen, zijn vele vragen niet te beant- 
woorden met de nu beschikbare kennis. Door de complexe structuur van de botsings- 
term zijn er praktisch geen exacte oplossingen van deze niet-lineaire vergelijking bekend. 
Het is daarom van het grootste belang vereenvoudigde modellen te bestuderen waar- 


mee de gerezen vragen beantwoord kunnen worden. Door puur toeval, naar het schijnt, 


isin 1971 door Bobylev, Krook en Wu een speciale exacte oplossing van de Boltzmann- 
vergelijking voor een bepaald geval gevonden. Deze oplossing heeft betrekking op een 
ruimtelijk homogeen gas van Maxwell-moleculen en heeft bovendien een heel eenvou- 
dige structuur. Krook en Wu uitten het vermoeden dat deze speciale oplossing een fun- 
damentele rol zou kunnen spelen in de beschrijving van de gang naar evenwicht. Ten- 
gevolge hiervan ontstond een enorme opleving in de belangstelling voor exact oplosbare 
modellen, die tot interessante resultaten geleid hebben. Over deze speciale oplossing, 
over het vermoeden en de daaruit voortvloeiende ontwikkelingen wordt hier verslag 


gedaan. 


In 1872 stelde Boltzmann! zijn welbeken- 
de niet-lineaire integro-differentiaal ver- 
gelijking op. Deze nog steeds meest fun- 
damentele vergelijking in de kinetische 
theorie van verdunde gassen? beschrijft 
de tijdsevolutie van de ééndeeltjes-verde- 
lingsfunctie f(w‚r, t) onder invloed van 
binaire botsingen tussen de deeltjes in het 
gas, en luidt in zijn meest gebruikelijke 
vorm: 


òf tu. Vf= 
=fdvi Sdûgl(g,. If fi —ffil- 


Hier isf du dr het aantal deeltjes in het 
volume-elementje du dr ten tijde f, 

dû = sin x dx dp de infinitesimale ruimte- 
hoek en f een éénheidsvector in de ver- 
strooiingsrichting. Verder is f, =f(v1, 
r,t)enz. De snelheden (u ‚v1) na de bot- 
sing: 


! 1 A 
v=G+i8 
/ Ì DN 
vi=G- 78 


gehoorzamen de behoudswetten van 
energie en impuls met G = 5 +U,). 


In de structuur van de botsingsterm uit de 
Boltzmann-vergelijking (B.V .) is het meest 
essentiële element de Stosszahlansatz. 
Deze zegt dat de kans op een binaire bot- 
sing evenredig is met het produkt van de 
één-deeltjes-verdelingsfuncties van de 
botsende deeltjes. De botsingsterm bevat 
de differentiële werkzame doorsnede 
1(g,x) waaring = |v —v, |de relatieve 
snelheid en x de verstrooiingshoek is. 

In geval het gas ruimtelijk homogeen is, 
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| ontbreekt de gradiënt term Vf in de B.V. 


Het totale aantal deeltjes N en de totale 
kinetische energie 5 NKT zijn dan behou- 
den grootheden. We meten nu snelheden 
in eenheden KT. De uiteindelijke even- 
wichtsverdeling f(u, oo) is de Maxwellse 
snelheidsverdeling — exp (-—v®/2). 
Kennen we eenmaal f als functie van snel- 
heid, plaats en tijd voor gegeven begin- en 
randvoorwaarden, dan hebben we een 
praktisch volledige beschrijving van de 
gang naar evenwicht zowel in ruimtelijk 
homogene gassen (relaxatieprocessen), als 
ook in ruimtelijk inhomogene gassen 
(hydrodynamica, transporteigenschappen, 
schokgolven). 


Helaas is de oplossing voor algemene inter- 
moleculaire krachten niet bekend. De 
moeilijkheid hierbij is vooral de ingewik- 
kelde structuur van de botsingsterm, die 
kwadratisch is in f. 

Wegens het ontbreken van exacte oplos- 
singen is het van groot belang vereenvou- 
digde modellen voor de niet-lineaire B.V. 
te bestuderen. Een vereiste is enerzijds 
dat de essentiële structuur met belang- 
rijke eigenschappen bewaard blijft, zoals 
behoudswetten en H-theorema. Ander- 
zijds moet het mogelijk zijn de algemene 
oplossing, of althans speciale oplossingen, 
te vinden. 


De recente ontwikkelingen® rond de niet- 
lineaire B.V. zullen we illustreren aan de 
hand van twee typische voorbeelden: 
Maxwell-moleculenen super-hardedeeltjes. 
Beide zijn modellen met een continue 
snelheidsvariabele. Ook beperken we ons 
tot ruimtelijk uniforme systemen geken- 
merkt door relaxatieprocessen. 





Prof. dr. M.H. Ernst (46 ) studeerde natuur- 
kunde aande Universiteit van Amsterdam. 
Zijn promotieonderzoek werd onder lei- 
ding van prof. dr. E.G.D. Cohen uitge- 
voerd, deels op het Instituut voor Theo- 
retische Fysica te Amsterdam, deels op 
Rockefeller University in New York. Hij 
promoveerde in 1965 bij prof. dr. J. de 
Boer. Vervolgens was hij verbonden aan 
het National Bureau of Standards te 
Washington, de University of Maryland en 
de Katholieke Universiteit van Nijmegen. 
In 1972 werd hij benoemd tot lector in 
de theoretische natuurkunde aan de Rijks- 
universiteit van Utrecht, waar hij nog 
steeds werkzaam is. Door de jaren heen 
richtte zijn belangstelling zich vooral op 
de niet-evenwicht statistische mechanica: 
transporteigenschappen in gassen en vloei- 
stoffen, fluctuatieverschijnselen, korte en 
lange tijdsgedrag van correlatiefuncties. 
Gedurende de laatste jaren heeft hij veel 
gewerkt aan exacte oplossingen van de 
Boltzmann-vergelijking. 


Maxwell moleculen 

De complexe structuur van de botsings- 
term in de B.V. wordt mede bepaald door 
de botsingsfrequentie g/(g, x). Indien de 
deeltjes elkaar afstoten met een kracht 
evenredig met 7=S (r is hun onderlinge af- 
stand) dan is g/(g, x) onafhankelijk van 
de snelheid als de exponent s = 5. Derge- 
lijke deeltjes worden Maxwell-moleculen 
genoemd. In feite is de x-afhankelijkheid 
van de botsingsfrequentie niet essentieël 
voor onze beschouwingen. Reden waarom 
we isotrope verstrooiing aannemen (dit is 
een mathematisch model; er bestaan im- 
mers geen deeltjes met een constante 
waarde voor g/(g, x)). Door nu tijden te 
meten in eenheden 7 — de gemiddelde 
vrije tijd tussen twee botsingen — neemt 
de B.V. voor het ruimtelijk homogene 
geval de vorm aan: 


dû ll 
of= fav | A AN 


terwijl voor (v'‚v, ) de standaardrelaties 
gelden. Aan de hand van deze vergelijking 
kunnen we enige recente ontwikkelingen 
toelichten. 


BKW-mode 
In 1976, ongeveer 100 jaar na Boltzmann 
en Maxwell, ontdekken Bobylev*, en 
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Fig. 1. De BK W-oplossing R(v, t) = 
fslv, t)/flv, eo) voor Maxwell-moleculen 
als functie van v voor verschillende tijden 
met T = NMt—to). 


Krook en Wu5 een speciale oplossing van 
| de B.V. voor Maxwell-moleculen. Deze 
oplossing — naar hun vinders BKW-mode 
genoemd — luidde een reveille in voor de 
belangstelling in exacte oplossingen van 
de B.V. Sinds 1978 zijn reeds een 60-tal 
publikaties over dit onderwerp in de 
literatuur verschenen, waarvoor we ver- 
wijzen naar een recent overzichtsartikel® 
De belangstelling werd nog extra geprik- 
keld doordat Krook en Wu een vermoe- 
den formuleerden over de mogelijk fun- 
damentele betekenis van deze oplossing, 
waarover straks meer. 

De BKW-mode f,(v, t) heeft een zeer 
eenvoudige structuur: 


—v? /2a 


RS 
(27a)2/2 


fw, t= 


Hierin zijn A en B positieve constanten, 
die met behulp van deeltjes- en energiebe- 
houd uit te drukken zijn in «, waarbij 


a(t)= 1 — exp ME — f0)]. 


Verder is A een eigenwaarde, die voor ons 
model gelijk is aan 1/6. Door substitutie 
in de bovenstaande B.V. is de oplossing 
rechtstreeks te verifiëren. Hij geldt ook 
voor echte Maxwell-moleculen, waarbij 

\ een gewogen hoekgemiddelde is van de 
werkzame doorsnede. 


De BKW-mode bevat slechts één willekeu- 
rige constante, de fase ft, die een klasse 
van begintoestanden karakteriseert. Meer 
details zijn te zien in figuur 1, waar de 
verhouding R‚(v, £) = fs(v, £) / f(w, co) is 
uitgezet als functie van v voor verschillen- 
de tijden. Voor grote tijden nadert R‚(v, £) 
naar de evenwichtsoplossing R (wv, co) = 1, 
waarbij vooral opvalt de extreem langzame 
(niet-uniforme) nadering tot evenwicht 
bij hoge snelheden. 
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Vermoeden van Krook en Wu 

We bekijken een beginverdeling R(v, 0), 
waar de hoge-energiestaart nog aan ont- 
breekt, dat wil zeggen dat R(v, 0) voor 
grote v sneller naar nul gaat dan een Max- 
well-verdeling. Dit is bijv. het geval voor 
de beginverdeling in de BKW-mode. De 
vraag is nu op welke manier de Maxwellse 
staartverdeling uiteindelijk tot stand 
komt. Dit is bijvoorbeeld van belang voor 
chemische en nucleaire reacties in gassen, 
die alleen boven een bepaalde (hoge) 
drempelwaarde plaatsvinden. Krook en 
Wu denken met name aan kernfusie, waar 
slechts een paar botsingen tussen hoog- 
energetische deeltjes nodig zijn om de 
reactie op gang te brengen. Extreem kleine 
concentraties kunnen dus van belang zijn. 
Zij speculeren nu dat een willekeurige be- 
gintoestand zonder hoge-energiestaart 
eerst relaxeert naar een toestand welke in 
essentie voorgesteld kan worden door de 
BKW-mode. Het daarop volgende stadium 
in de tijdsontwikkeling wordt in essentie 
beschreven door de BKW-mode met een 
geschikt gekozen fase ft. Een wat exactere 
formulering luidt: bij geschikt gekozen 
parameter f, bestaat een t, zodanig dat 
voor alle t > t,‚ geldt dat | R(v, ) — 
Rv, t)|<e, uniform in v. 





Is dit vermoeden waar? De eerste studies 
hierover zijn gemaakt door Tjonen Wué 7 
door het numeriek oplossen van de B.V. 
voor een model van deeltjes met twee 
translatie-vrijheidsgraden, dat nauw sa- 
menhangt met het hierboven besproken 
Maxwell-model met isotrope verstrooiing. 
Typische resultaten zijn te zien in figuur 
2, waar de verhouding R(x, f) = fw, #)/ 
f(v, oo) te zien is als functie van t bij ver- 
schillende waarden van de energiex =5v?. 
Energieën wòrden hier gemeten in een- 
heden KT. Twee begintoestanden, beide 
zonder hoge-energiestaart, vervallen bij 
hoge energieën op verschillende wijze. De 
ene vervalswijze (stippellijnen) is kwalita- 
tief in overeenstemming met de BK W- 
mode. In figuur 1 zien we immers dat de 
BKW-mode, bij voldoende hoge energie, 
de evenwichtsverdeling altijd van de on- 
derkant nadert. De vervalswijze van de 


Fig. 2. De numerieke oplossing’ R(x,t) 
= flv, t)/f{v, oo) van een model-B. V. als 
functie van t bij verschillende energieën 
x= zv? Iedere beginverdeling bestaat uit 


twee deltapieken: voor de stippellijnen 


staan deze op x = 0,4 en 2,0; voor de ge- 
trokken lijnen op x = 0,4 en 3,6. 
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andere begintoestand (getrokken lijnen) 
is kwalitatief verschillend van de BKW- 
mode: de aanvankelijk lege staart wordt 
zeer snel bevolkt, en de bezetting schiet 
over de uiteindelijke evenwichtswaarde 
heen. Er ontstaat een tijdelijke overbevol- 
king van hoge-energietoestanden, soms 
een factor 10 groter dan de uiteindelijke 
evenwichtswaarde. Dit kan weer van be- 
lang zijn in chemische en nucleaire reacties 
met een bepaalde activeringsenergie. In 
die tijdelijke overbezetting zullen reactie- 
snelheden aanzienlijk groter zijn dan 
schattingen gebaseerd op de Maxwellse 
snelheidsverdeling. 

Het laatste relaxatieproces is niet noodza- 
kelijkerwijs in tegenspraak met het ver- 
moeden van Krook en Wu. De mogelijk- 
heid bestaat dat voor dit model het sta- 
dium van nadering tot de BKW-mode lan- 
ger duurt dan 30 7. In dat geval laat 
figuur 2 zien dat het vermoeden triviaal 
waar is voor x-waarden tot ongeveer 100 
maal de thermische energie, omdat R en 
R, praktisch de evenwichtswaarde 1 be- 
reikt hebben binnen een tijd van 30 7. 
Hoe de relaxatie verloopt bij veel hogere 
x-waarden is nog steeds een open vraag. 
Echter, ook al zou het vermoeden juist 
zijn, toch zijn de hier gevonden concen- 
traties e-100 — 1075 te klein om van enig 
praktisch nut te zijn. Daarbij dient wel 
opgemerkt te worden dat het juist in niet- 
lineaire vergelijkingen moeilijk schatten is 
in hoeverre de keuze van het model (aan- 
tal translatie-vrijheidsgraden, botsings- 
frequentie) de grootte-orde van de rele- 
vante energiewaarden kwalitatief kan wij- 
zigen. Er bestaan verschillende analytische 
argumenten tegen de geldigheid van het 
vermoeden. Deze zijn gebaseerd op de 
eerste paar termen in de Sonine-reeksen 
voor de verdelingsfunctie (zie hierna), en 
lenen zich derhalve niet tot uitspraken 
over het gedrag van de verdelingsfunctie 
bij willekeurig hoge energiewaarden. 
Samenvattend kan gezegd worden dat de 
speciale BKW-oplossing voor Maxwell- 
moleculen van weinig praktisch belang 
lijkt voor de beschrijving van de algemene 
gang naar evenwicht zelfs bij relatief hoge 
energieën. Tot nu toe ontbreken echter 
een rigoreus bewijs vóór en onomstote- 
lijke argumenten tegen de geldigheid van 
het vermoeden. 


Ex pansies in Sonine-Polynomen 

Voor het gemak beperken we de discussie 
verder tot isotrope verdelingen, alleen af- 
hankelijk van v = lvl. Een verdeling kan, 
behalve door f(v, t) zelf, ook beschreven 
worden door zijn reeks momenten. Nu is 
het in het algemeen niet mogelijk uit de 
B.V. een stelsel (niet-lineaire) bewegings- 
vergelijkingen voor deze momenten af te 
leiden, die uitsluitend de momenten zelf 
bevatten. 

Voor Maxwell-moleculen is dit echter wel 
het geval, en het speciale er van is dat de 
vergelijking voor het n-de moment alleen 
de momenten 0 tot en met 7 bevat. De 
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momentenvergelijkingen zijn dus succes- 
sief oplosbaar, zoals in 1956 door Trues- 
dell gevonden is®. 

Soortgelijke resultaten, eveneens in het 
geval van Maxwell-moleculen, zijn in de 
laatste jaren ook verkregen® voor de zoge- 
naamde Sonine-momenten c,‚(t), gedefi- 
nieerd als de coëfficiënten in de reeks- 
ontwikkeling: 


fw, £) = f(v, 00) B en (D SW). 


De Sonine-polynomen S,(v) vormen een 
stelsel orthogonale polynomen, waarvan 
de precieze vorm hier niet ter zake doet. 
Wegens deeltjes- en energiebehoud is ver- 
dercg =l enc, =0. Ter illustratie geven 
we de bewegingsvergelijkingen voor de 
Sonine-momenten van bovenstaand 
Maxwell-model met isotrope verstrooiïng. 
Deze luiden: 


Ë +en=(n +1)! X Ck Cn kt 
n=0 


Alle c,‚’s zijn in principe uit te rekenen 
voor willekeurig grote n, en daarmee is 
ook de verdelingsfunctie zelf be paalbaar. 
Zijn we alleen geinteresseerd in de verde- 
ling van energieën tot bijvoorbeeld 15 KT, 
dan kan een hoge nauwkeurigheid voor 
flv, t) bereikt worden door de Sonine- 
reeks af te breken bij, zeg, n = 18. Het 
aantal termen dat bijdraagt tot cg is 88. 
Voor een redelijk nauwkeurige weergave 
van relaxatieprocessen in de hoge-energie- 
staart van de verdeling zou n echter zeer 
groot moeten zijn, bijvoorbeeld n = 10°. 
Dit is ver boven de capaciteit van de tegen- 
woordige rekenmachines. De reden is dat 
het aantal termen, dat bijdraagt tot c,, 
zeer sterk toeneemt met stijgende n ten 
gevolge van de niet-lineariteit in het rech- 
terlid van de c‚-vergelijking. Dit aantal is 
bijvoorbeeld voor n = 10° reeds van de 
orde 10%. | 


Inmiddels zijn er voor vele Maxwell-ach- 
tige modellen oplossingen van de niet- 
lineaire B.V. verkregen in de vorm van 
reeksontwikkelingen in orthogonale poly- 
nomen. Oplossingen in gesloten vorm zijn 
echter niet bekend voor dergelijke model- 
len, behalve de reeds besproken BKW- 
mode f‚(v, t) en de Maxwell-verdeling 
fu, co). 

Dit is wel het geval voor een speciaal 
model van super-harde deeltjes, dat we 
bespreken in het resterende deel van dit 
artikel. Analytische oplossingen zijn van 
het grootste gewicht omdat ze antwoor- 
den kunnen geven — al is het slechts voor 
een speciaal model — op de boven gerezen 
vragen en problemen: relaxatie processen 
bij alle mogelijke energieën; niet-unifor- 
miteiten in de nadering tot evenwicht; de 
energie-afhankelijkheid van de bijbeho- 
rende karakteristieke tijden; tijdelijke 
overbezettingen van hoge-energietoe- 
standen enz. 


super-harde 
deeltjes 





harde 
schijven 





TL 0 TL 


Fig. 3. Hoekafhankelijkheid van de werk- 
zame doorsnede Ig, x) voor twee twee- 
dimensionale modellen. 


Super-harde deeltjes | 
De essentie van het model van super-harde 
deeltjes, ingevoerd door Hendriks en 
schrijver dezes’, is dat de deeltjes be we- 
gen in een plat vlak en dat hun botsingen 
beschreven worden door de werkzame 
doorsnede I(g, x) —g | sin x |. 

In binaire botsingen zijn energie en impuls 
behouden. Ook laat figuur 3 zien dat de 
hoekafhankelijkheid van de werkzame 
doorsnede voor dit model analoog is aan 
die voor harde schijven. De kwalificatie 
super-hard is ontleed aan de eigenschap 
dat Z(g, x) lineair toeneemt met g. Bij _ 
harde bollen en schijven is (2, x) onaf- 
hankelijk van g. Bij zachtere interacties, 
zoals voor Maxwell-moleculen, neemt 


1(g, x) af met g. 


Een nadeel van het model is dat /(g, x) te 
snel toeneemt met g. Dit is geen realis- 
tische eigenschap van intermoleculaire 
wisselwerkingen. Omdat het verstrooiings- 
mechanisme te efficiënt is bij hoge ener- 
gieën, zijn de typische relaxatietijden bij 
die energieën in het algemeen te kort. 
Beperken we ons weer tot isotrope snel- 
heidsverdelingen, dan kan de niet-lineaire 
B.V. voor dit model in gesloten vorm 
worden opgelost voor algemene begin- 
voorwaarden. Om de algemene oplossing 
te vinden is het beter over te gaan op de 
energie-verdelingsfunctie F(x, t), gede- 


| finieerd door de relatie f(u, f) 2nvdv = 


F(«, t)dx. Hierin isx = jv° de energie, 
gemeten in eenheden kT. Als we verder 
het totale aantal deeltjes en de totale 
energie op 1 normeren, kan de B.V. voor 
dit model geschreven worden als: 
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(B, +x+1) Fx) = 


Ie du [* dy Fy) Fu — 9). 


Deze vergelijking kan herleid worden tot 
een niet-lineaire partiële differentiaal- 
vergelijking van eerste orde, die oplosbaar 
is, zodat alle gewenste informatie berekend 
kan worden. Allereerst een eenvoudig 
voorbeeld? van een beginverdeling, die 
een superpositie is van twee Maxwell-ver- 
delingen exp[—x/7;(0)] met verschillende 
temperatuur 7,(O) (i = 1,2) met 7, (0) + 
T,(O) > 1 en 7,(O) <1. De verdelings- 
functie behoudt deze vorm in de loop van 
de tijd: 


F«,‚t)= Za;(t) exp[—x/7;(6)]. 


De a;(t)’s hangen af van 7, (t) en 7, (t), 
die zelf bepaald worden als wortels uit 
een vierkantsvergelijking met bekende tijd- 
afhankelijke coëfficiënten. Voor grote f 
nadert F(x, t) naar de Maxwell-verdeling 
F(x, oo) =e”*, doordat 7, (t) naar 1 na- 
dert, en 7, (t) omgekeerd evenredig met 

t afneemt. 

In dit voorbeeld zijn twee kwalitatief ver- 
schillende gevallen te onderscheiden: Als 
T,(O) <1,danisop t= 0 de hoge-energie- 
staart dun bevolkt en gedraagt R(x, t) 
zich kwalitatief als in figuur 1. De tijd- 
schalen zijn echter korter. Als daarentegen 
T,‚(O) > 1, dan hebben we op f =O een 
overbevolkte staartverdeling en R(x, t) 
zal voor grote energieën van de bovenkant 
naar 1 naderen. 

Tenslotte bespreken we nog hoe een wil- 
lekeurige begintoestand tot evenwicht na- 
dert. Voor voldoend grote x en t geldt: 


R(x‚t)=e F(x, t) =zexcd, 


waarin c(f) samenhangt met het tweede 
moment van de energieverdelingsfunctie, 
volgens c(f)= 5 <x? >— 1. In dit model 
kan c(t) expliciet berekend worden, en 
wordt voor grote t van de orde ef. De 
coëfficiënt c(t) kan n nulpunten (n = 0, 
1, 2 ...) hebben, afhankelijk van de ge- 
geven beginverdeling. De verdelingsfunc- 
tie zelf zal dan n maal(n=0,1,2 ...) 
met afnemende amplitude om zijn even- 
wichtswaarde oscilleren, zoals geschetst 
in figuur4 voorn = 2. Het grote x-t-gedrag 
is hier dus geheel verschillend van dat in 
de BKW-mode. 


Niet-uniformiteiten 

De verdelingsfunctie nadert niet-uniform 
tot de Maxwell-verdeling. Dit manifesteert 
zich door een toename van karakteristieke 
tijden in hogere energieën. Voor lagere 
energieën, tot enige malen K7’, is de typi- 
sche relaxatietijd natuurlijk de vrije-weg- 
lengte 7. De toename van de karakter- 
istieke tijden is als volgt in te zien. 


Na een binaire botsing is de totale energie 
van de twee beschouwde deeltjes in het 
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Fig. 4. Oscillerend gedrag van de verde- 
lingsfunctie in het superharde-deeltjes- 
model, bij vaste (willekeurige grote) ener- 
gie X. 


gunstigste geval geconcentreerd op één. 
Stel, de staart van de verdeling is aanvan- 
kelijk onbezet boven een zekere grens- 
waarde, dan kan de grenswaarde hoogstens 
verdubbelen na iedere gunstige botsing. 
Er zijn dus vele botsingstijden 7 nodig 
om energiewaarden te bereiken, die ver 
boven de initiële grenswaarde liggen. 
Voor Maxwell-moleculen kunnen de 
karakteristieke tijden 7 (wv) eenvoudig 
geschat worden met behulp van de BKW- 
oplossing, welke voor grote waarden van 
v en t de vorm heeft: R‚(v, t) Ze" x 

(1 +n) met n(t) = zu? exp(—)\t). De ka- 
rakteristieke tijd volgt nu uit de relatie 
n(T,) = 1, zodat: 


Tow) = NÀ log v?. 


Voor t >> To(w) geldt R‚(w, t) = 1 en 
voor ft < To (wv) geldt R‚(w, £) = 0. 

Een niet-uniforme nadering tot evenwicht 
met typische tijdschalen 7, (w) == log v* is 
ook te zien in het super-harde deeltjes- 
model, terwijl bij deeltjes met interacties 
zachter dan bij Maxwell-moleculen tijd- 
schalen voorkomen? als Tow) — v* of v°. 
Figuur 1 laat verder zien dat er voor een 
vaste willekeurig grote tijd steeds ener- 
gieën te vinden zijn, waarvoor R‚(v, t) 
niet in de buurt van de evenwichtswaarde 
ligt. We kunnen de oplossingen dus niet 
linearizeren om evenwicht. De niet-unifor- 
miteiten in andere modellen leiden tot 
dezelfde conclusies. 

De niet-uniforme nadering tot de Maxwell- 
verdeling laat zien dat de gelinearizeerde 
B. V. de relaxatieprocessen in de hoge- 
energiestaart van de verdelingsfunctie niet 
juist beschrijft. Tenslotte is het interessant 
om op te merken dat de hier bestudeerde 
niet-lineaire model-vergelijkingen ten 
nauwste samenhangen met de vergelijkin- 
gen ter beschrijving van coagulatie-pro- 
cessen!®, zoals kunnen optreden bij de 
vorming van melkwegstelsels, sterren, 


wolken en stofdeeltjes, alsook bij collo- 
idale deeltjes, polymeren en gels, tot zelfs 
bij het stollen van bloed. 
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Kinetische theorie van 
polymerenvloeistoffen 


R. Byron Bird en C.F. Curtiss 


De algemene kinetische theorie voor polymerenvloeistoffen van Curtiss, Bird en Hassager 
verschaft een gemeenschappelijk uitgangspunt voor het beschrijven van de rheologische 
eigenschappen van verdunde polymeeroplossingen en die van polymeersmelten. Deze 
theorie omvat de vroegere theorieën van onder meer Kramers, Kirkwood c.s. en Zimm, 
voor verdunde oplossingen als speciale gevallen, alsmede de onlangs gepubliceerde 
theorie van Doi en Edwards voor polymeersmelten. De theorie van Curtiss, Bird en 
Hassager geeft de mogelijkheid modellen te gebruiken voor polymeren, die uit wille- 
keurige verzamelingen van bolletjes, staafjes, en veertjes bestaan. Hierdoor is het moge- 
lijk het verband tussen rheologische eigenschappen en polymeerstructuur te bestuderen. 


Kinetische theorie voor de viscositeit 

van verdunde gassen 

Wij beginnen met een kort overzicht van 
de klassieke kinetische theorie van ver- 
dunde, éénatomige gassen. Hierin worden 
de moleculen als puntmassa’s (massa 77) 
beschouwd, waarvan de wisselwerking 
door bijvoorbeeld een Lennard-Jones 
6-12-potentiaal met parameters o en e 
beschreven kan worden. Het gas wordt 
dan voorgesteld zoals in figuur la, waarin 
de puntmassa’s (n puntmassa’s per vo- 
lume-eenheid) door de ruimte bewegen en 
af en toe met elkaar botsen. 


De kinetische theorie begint met de 
Liouville-vergelijking voor de ensemble 
distributiefunctie, en bestaat uit twee 
delen: í. De Liouville-vergelijking kan ge- 
integreerd worden om een reeks verge- 


Fig. 1. Modellen in de kinetische theorie: 
a. verdunde, éénatomige gassen; b. ver- 
dunde polymeeroplossingen-polymeer- 
ketens met N = 10 in een zee van struc- 
tuurloze oplosmiddelmoleculen; c. poly- 
meersmelt met polymeerketens (N = 10) 
die met elkaar in wisselwerking zijn. 








Nederlands Tijdschrift voor Natuurkunde A47 (4) 1981 


lijkingen te verkrijgen die de hiërarchie 
van BBGKY genoemd wordt. De eerste 
vergelijking van deze reeks is een verge- 
lijking in de faseruimte van één molecule, 
die echter de distributiefunctie bevat van 
twee moleculen. Om de Boltzmann-verge- 
lijking te verkrijgen wordt de zogenaamde 
moleculaire-chaos-benadering ingevoerd 
waarmee de reeks beëindigd kan worden; 
op die manier wordt een gesloten verge- 
lijking verkregen voor de distributie- 
functie f(r,p, t), die de kans beschrijft 
dat op het tijdstip f een molecule met 
impuls p zich op de plaats r bevindt. Deze 
vergelijking wordt vervolgens opgelost 
met een reeksoplossing, die voor een iso- 
therm gas luidt: — 
fr‚p.bd=fe- (l -B:Vu+ …), (1) 
waarin v de gassnelheid is; de tensor B bevat 
ingewikkelde integralen die afhankelijk 
zijn van de details van de dynamica van 
binaire botsingen. ii. De bewegingsverge- 
lijking van het gas kan ook van de Liouville- 
vergelijking afgeleid worden, waarbij een 
uitdrukking voor de spanningstensor ver- 
kregen wordt. Voor verdunde gassen luidt 


deze uitdrukking: 
— 2) ) ‚ ÂA) 


hes AE 


Hier betekent & >de gemiddelde waarde 
ten opzichte van de distributiefunctie 
f(r‚p, t). Volgens (2) vindt impulstrans- 
port uitsluitend plaats door de beweging 
van de moleculen (figuur 2a). 


Uit continuümmechanische overwegingen 
is het bekend dat, in lineaire benadering, 
de spanningstensor de volgende vorm 
moet aannemen: | 
AN ee 
n=Ppò —ny+(in—K)(V: v)ò, (3) 
waarin p de druk is,ö de eenheidstensor, 


y= Vv + (Vov)t een deformatiesnelheids- 
tensor, n de viscositeit, en Kk de volume- 


viscositeit. Als men de oplossing (1) in (2) 
substitueert en het resultaat met (3) ver- 
gelijkt, verkrijgt men een uitdrukking voor 
de viscositeit n waarin de parameters o en 
e voorkomen. Bovendien vindt men tege- 
lijkertijd dat k nul is voor verdunde, één- 
atomige gassen. 

De Boer en van Kranendonk! hebben, ge- 
bruik makend van deze theorie, de trans- 
porteigenschappen uitgerekend voor een 
gas waarin de moleculaire wisselwerking 
door een Lennard-Jones 6-12-potentiaal 
beschreven kan worden. De Boer en andere 
medewerkers hebben ook soortgelijke 
quantummechanische problemen onder- 
zocht voor He® 2, He3 3, ortho- en para- 
waterstof* , gassen bij hoge tempera- 
turen” en voor meeratomige gassen® . De 
resultaten van deze berekeningen zijn uit- 
voerig beschreven in MTGL', waaraan 
De Boer twee hoofdstukken heeft bijge- 
dragen. 

De rest van dit artikel is gewijd aan de 
kinetische theorie van polymerenvloei- 
stoffen, dat wil zeggen polymeeroplossin- 
gen en niet-verdunde polymeren (poly- 
meersmelten). Het zal blijken dat de volg- 
orde in de ontwikkeling van de theorie 
analoog is aan die in de overeenkomstige 
theorie voor gassen: de keuze van een 
model voor het molecule; het beschrijven 
van de vloeistof, het verkrijgen van een 
uitdrukking voor de spanningstensor en — 
uiteindelijk het uitrekenen van de vis- 
cositeit (en andere rheologische eigen- 
schappen) waarin de modelparameters 
voorkomen. 


Modellen voor polymeren en 
polymerenvloeistoffen 

Er zijn vele soorten polymeren, maar in — 
deze discussie beperken we ons tot keten- 
achtige macromoleculen zonder zijketens. 


Fig. 2. Bijdragen tot de spanningstensor: 
a. impulstransport door de moleculen in 
een verdund gas; b. vijf bijdragen ineen 
polymeeroplossing-ss en Typ komen in 
polymeersmelten niet voor. 
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Wij kiezen als een eenvoudig maar redelijk 
model het ‘parelsnoer’ van Kramers? ; het 
macromolecule wordt dan voorgesteld als 
een lineaire, uit N bolletjes (massa m) en 
Nt staafjes (lengte a) bestaande keten, 
waarbij de staafjes volkomen vrij draai- 
baar zijn. 

De oriëntatie van staafje k in de Kramers- 
keten kan door de twee bolcoördinaten 
0, en by worden gegeven, òf door een 
eenheidsvector uy. Bovendien kunnen 
twee andere eenheidsvectoren worden 
gedefinieerd: vx in de 0,-richting, en w‚ 
in de py-richting. Voor ieder staafje in de 
keten kan een operator: 


Ò Ì d 
men 
i sin 0 Ò dx 


Ò an 
Tee @ 
worden gedefinieerd. 

Een verdunde polymeeroplossing wordt 
vervolgens voorgesteld (figuur lb) als een 
mengsel van parelsnoerketens (de opge- 
loste polymeermoleculen)en een verzame- 
ling van puntmassa’s (de oplosmiddel- 
moleculen). De concentratie van de op- 
geloste polymeermoleculen wordt veron- 
dersteld zó klein te zijn dat deze polymeer- 
moleculen praktisch nooit met elkaar bot- 
sen en het is derhalve voldoende alleen 
maar de beweging van één polymeerketen, 
die zich in een zee van oplosmiddelmole- 
culen verwringt, te beschouwen; dit mole- 
134 





cule kan zich in alle mogelijke confor- 
maties bevinden. Een polymeersmelt 
wordt zoals in figuur 1c afgebeeld als een 
verzameling parelsnoeren, die voortdurend 
met elkaar in wisselwerking zijn. De bu- 
ren van een gegeven molecule zijn altijd 
dichtbij. Dit heeft een strenge beperking 
in zijdelingse bewegingen tot gevolg, maar 
een slangachtige beweging (Frans: “repta- 
tion”) in de richting van de ruggegraat 
van de keten is wel toegestaan. Hierbij 
wordt verondersteld dat de wisselwer- 
kingsenergie van twee ketens gelijk is aan 
de paarsgewijze som van de wisselwerkin- 
gen tussen de bolletjes van de twee ketens. 


Algemene theorie voor 
polymeeroplossingen en polymeersmelten 
De kinetische theorie voor verdunde op- 
lossingen en die voor polymeersmelten 
kunnen allebei opgesteld worden met één 
enkele formele theorie als uitgangspunt. 
Deze algemene theorie is door Curtiss, 
Bird en Hassager® ontwikkeld uit een ge- 
neralisatie van de vroegere werken van 
Kramers® en Irving en Kirkwood!®; een 
vereenvoudigde versie van de CBH-theorie 
kan gevonden worden in hoofdstuk 14 
van een leerboek over de kinetische theorie 
van polymeren!! . Een overzicht van de 
CBH-theorie is in figuur 3a gegeven. 

Wij beginnen met de bewegingsvergelijkin- 
gen van Hamilton voor alle puntmassa's in 
het systeem (dat wil zeggen het systeem 
van figuur 1b òf dat in figuur lc, gebruik 
makend van de coördinaten en impulsen 
van het massamiddelpunt van ieder poly- 
meermolecule, de gegeneraliseerde coör- 
dinaten en impulsen voor de inwendige 
vrijheidsgraden van het polymeermolecule 
en de coördinaten en impulsen van de op- 
losmiddelmoleculen (in zoverre deze aan- 
wezig zijn). Vervolgens wordt de Liouville- 
vergelijking opgesteld voor een distributie- 
functie in de faseruimte voor een ensemble 
van systemen; de Liouville-operator £ is 
uiterst ingewikkeld met afgeleiden naar 
alle faseruimtevariabelen. Als de Liouville- 
vergelijking vermenigvuldigd wordt met 
B,een willekeurige, tijdsonafhankelijke 
functie van coördinaten en impulsen en 
daarna geïntegreerd over alle faseruimte- 
variabelen, verkrijgt men de algemene 
transportvergelijking. 

We zijn nu in staat speciale keuzen voor 
B te maken. Wij kunnen bijvoorbeeld B 
zó kiezen dat we een vergelijking krijgen 
voor de tijdsontwikkeling van de distri- 
butiefunctie in de faseruimte van één 
molecule; deze vergelijking is een genera- 
lisering van de eerste vergelijking van de 
hiërarchie van BBGKY, en bevat de paar- 
distributiefunctie. Uit de vergelijking kun- 
nen twee belangrijke vergelijkingen wor- 
den afgeleid: i. door integratie over alle 
impulsvariabelen, kan een vergelijking 
worden afgeleid voor de tijdsontwikkeling 
van de configuratie distributiefunctie W 
voor één molecule (dat wil zeggen, de con- 
tinuiteitsvergelijking voor Y); üú. door te 
vermenigvuldigen met één van de gegene- 





Prof. dr. C.F. Curtiss promoveerde in 
1948 in scheikunde aan de University 
of Wisconsin in Madison, en is sinds- 
dien lid van het Theoretical Chemistry 
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try aan de University of Wisconsin. 

In 1960 werd hij tot hoogleraar be- 
noemd. Hij is een van de mede-auteurs 
van Molecular Theory of Gases and 
Liquids (uitgegeven door Wiley in 
1954). In de laatste jaren heeft hij met 
Bird, Hassager, en Armstrong meege- 
werkt aan Deel 2 (‘Kinetic Theory’) 

van het twee-delige boek Dynamics of 
Polymeric Liquids uitgegeven door Wiley 
in 1977). Zijn voornaamste research-in- 
teresse ligt in de verscheidene aspecten 
van niet-evenwicht statistische mecha- 
nica: transportverschijnselen in gassen; 
kinetische theorie van polymeeroplos- 
singen en polymeersmelten; de quantum 
en klassieke theorie van inelastische bot- 
singen tussen eenvoudige moleculen. 





raliseerde impulsen en daarna te integreren 
over alle impulsvariabelen, verkrijgt men 
een bewegingsvergelijking voor iedere in- 
wendige vrijheidsgraad van het polymeer- 
molecule. Uit deze twee vergelijkingen 
verkrijgen we uiteindelijk een tweede-orde 
partiële differentiaalvergelijking (de ‘dif- 
fusievergelijking’) voor de configuratie- 
distributiefunctie WP, maar voordat zoiets 
kan worden uitgevoerd is het nodig enkele 
veronderstellingen te maken, die straks 
besproken zullen worden. 

Een andere keuze voor B (figuur 3a) leidt 
tot de gewone hydrodynamische continu- 
iteitsvergelijking. Nog een andere keuze 
voor B levert de hydrodynamische bewe- 
gingsvergelijking op; in de bewegingsver- 
gelijking komt de spanningstensor voor, en 
deze wordt gegeven als de som van enkele 
nogal ingewikkelde integralen. De kine- 
tische bijdrage mx beschrijft het transport 
van impuls van de puntmassa’s (oplosmid- 
del zowel als opgeloste stof) en lijkt veel 
op vergelijking (2). De oplosmiddelbijdrage 
Tss Ontstaat uit de wisselwerking tussen 
twee oplosmiddelmoleculen. De bijdrage 
Tps beschrijft de intermoleculaire krach- 
ten tussen de oplosmiddelmoleculen en 
de opgeloste polymeermoleculen. De 
term 7, volgt uit de kracht die uitgeoe- 
fend wordt vanwege de spanningen in de 
schakels van de polymeerketen. En ten- 
slotte heeft spp te maken met de wissel- 
werking tussen twee polymeerketens. 
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Deze verscheidene bijdragen tot sr zijn 
geillustreerd in figuur 2b. Elke bijdrage 
wordt als een integraal gegeven waarin 

de bijbehorende intermoleculaire òf intra- 
moleculaire potentiële energie voorkomt; 
om deze uitdrukkingen te gebruiken is 
het echter nodig enkele benaderingen te 
maken. 


De diffusievergelijking voor de 
configuratie-distributiefunctie 

Integratie over alle impulsen om de bewe- 
gingsvergelijking te verkrijgen (zie figuur 
3a onderaan) geeft een term waarin de 
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paar-distributiefunctie voorkomt (Ys 
voor verdunde oplossingen òf W,, voor 
polymeersmelten). Een expliciete bereke- 
ning van WV òf W‚, kan vermeden wor- 
den door Ben bepaalde veronderstelling te 
maken. Wij merken op dat een gedeelte 
van de term die Ps, of V,, bevat, ge ïn- 
terpreteerd kan worden als een gemiddel- 
de kracht K, uitgeoefend op bolletje v 
van één polymeerketen door alle oplos- 
middelmoleculen (in een oplossing) òf 
door alle andere polymeermoleculen (in 
een smelt). Voor een verdunde oplossing 
nemen we aan dat K,, gegeven wordt door 


Fig. 3. Structuur van de kinetische 
theorie: a. algemene theorie voor poly- 
merenvloeistoffen, met nummers van 
vergelijkingen uit referentie 1 1;b. voort- 
zetting van de linker zijde van figuur 3a, 
met een beschrijving van de afleiding van 
vergelijkingen voor de configuratie-distri- 
butiefunctie, met nummers van vergelij- 
kingen uit referentie 11 of referentie 12; 
c. voortzetting van de rechter zijde van 
figuur 3a, met de afleiding van de span- 
ningstensor. De vector R,, duidt de plaats 
van bolletje v aan ten opzichte van het 
massamiddelpunt van het polymeer- 
molecule. < > is de gemiddelde waarde 
ten opzichte van Y, behalve in de laatste 
uitdrukking voor smelten, waarvoor < > 
de gemiddelde waarde ten opzichte van 
flu, o, t) is. Het onderschrift ‚}‚ is een 
codeformationele tijdsafgeleide!!, 12: 
Aiy=2Mdted. res (Vv)t . Â— A 

. (Vv). De bijdrage van het oplosmiddel 
totT is T.. 


de wet van Stokes (een dergelijke veron- 
derstelling wordt impliciet gemaakt in de 
theorie van Kramersë en wordt sindsdien 
gebruikt door alle onderzoekers): 


Ky, = {(([7,] — ú(r)), (5) 


dat wil zeggen de gemiddelde op bolletje 
vp uitgeoefende kracht is evenredig met 
het verschil tussen de ensemblegemiddelde 
snelheid van het bolletje, [[f, |], en de snel- 
heid van het oplosmiddel op de plaats r, 
van bolletje v; de dubbele haken [[ || dui- 
den de gemiddelde waarde aan in de im- 
pulsruimte van één polymeerketen. De co- 
efficiënt { is de ‘wrijvingscoëfficiënt’. 
Voor polymeersmelten gebruiken we een 
anisotrope vorm van de wet van Stokes, 
aangezien de aan een ketenschakel gebo- 
den weerstand in de u‚-richting niet ge- 
lijk is aan die in de v‚,- en wy -richtingen: 


Kja1 _Kj= MEP Ta 


— Ulrr+i) — vr), (6) 


{LF — fr] 


B = [8 —(l — uzurl. DD 


waarin v(r;) de snelheid is van de poly- 
meersmelt op de plaats r;, van bolletje k. 
In deze gemodificeerde wet van Stokes 
treden drie parameters op: een ‘wrijvings- 
coëfficiënt’ €; de ‘schakelspanningscoëf- 
ficiënt’ e; en de ‘ketenbeperkingscoëf- 
ficiënt’ 6. Voor e = 1 wordt vergelijking 
(6) weer isotroop;en voore = Oen p= 0 
worden de eindresultaten voor de span- 
ningstensor bijna dezelfde als die van Doi 
en Edwards!®, wier theoretische ontwik- 
keling echter sterk van onze methode 
afwijkt. De constante @ is ongeveer 0,3 
— 0,4 en wordt empirisch bepaald uit de 
afhankelijkheid van de viscositeit van het 
moleculairgewicht. (N.B. voor verdunde 
oplossingen is K , gelijkwaardig met F, in 
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referentie 11, maar voor polymeersmelten 
is K, gelijk aan N maal de F, van referen- 
tie 12.) 
Nadat de empirische uitdrukkingen voor 
K, in de bewegingsvergelijking zijn inge- 
voerd, kunnen deze en de continuiteits- 
vergelijking voor W gecombineerd worden 
om een tweede-orde partiële differentiaal- 
vergelijking voor W te verkrijgen. Deze 
vergelijking bevat echter de gemiddelde 
waarden [[úxu;]] in de impulsruimte. Hier 
zijn we weer verplicht een veronderstelling 
te maken om de vergelijking om te zetten 
in de configuratieruimte van één polymeer- 
| keten. Voor verdunde oplossingen wordt 
verondersteld dat de impulsdistributie-. 
functie in een niet-evenwichtstoestand ge- 
lijk is aan die in een evenwichtstoestand. 
Als we verder aannemen dat uz +1 — Up 
erg klein is (geringe kromming), verkrijgen 
we: 


Li; = KT/ma? Aij © —uzuj), (8) 


waarin de A‚; de Rouse-matrixelementen 
zijn: Aj; = 2 voork=j,= —1 voork =| 

+ 1;= 0 in alle overige gevallen. Voor poly- 
meersmelten wordt aangenomen dat de 
polymeerketens zich gemiddeld alleen 
heen en weer kunnen bewegen in de rich- 
ting van de ruggegraat van de keten (slang- 
beweging): 


Lia; = (2kT/ma® )ôu,Ööuj, (9) 


waarin ug = (1/2) (uy+1 — Ux 1); voor 
de uiteinden van de keten moet vergelij- 
king (9) enigszins gewijzigd worden. Na 
het invoeren van (8) en (9) verkrijgt men 
de diffusievergelijking voor de configuratie- 
distributiefunctie Y. 

Deze diffusievergelijking kan voor poly- 
meersmelten over alle u;’s behalve u, ge- 
integreerd worden, waaruit een vergelij- 
king ontstaat voor de distributiefunctie 
fx, tf) voor de oriëntatie van één schakel; 
deze functie geeft de kans dat schakel k 
zich in de oriëntatie u bevindt op een tijd- 
stip £. Hierna wordt de discrete keten 
door een continue keten vervangen; we 
gebruiken nu de variabele o in plaats van 
k/N, de relatieve afstand langs de keten, 
en fx (u, t) wordt vervangen door f(lu,‚o,£), 
waarvoor de differentiaalvergelijking luidt: 


2 
EE 


òf _ 1 
dt _\ do? 
ò 


& 


waarin \ = N3*Ê Ca? /2kT een karakteris- 
tieke tijdconstante is voor de vloeistof. 
De oplossing van (10) met randvoorwaar- 
den f= 1/47 bij a = 0,1 (voor alle t) en 
bij f= —oo (voor alle o) is: 


[u + Vv — uuu: var) ‚ (10) 


1 
fu,o, D= el 
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hm 


Hierin is y[Ol (rt, f°) een deformatietensor 


voor niet-oneindig kleine deformaties!!, 12 


Een dergelijke ontwikkeling voor verdun- 
de oplossingen is nog niet uitgevoerd. 


De spanningstensor 

In figuur 3c wordt een overzicht gegeven 
van de afleiding van een uitdrukking voor 
de spanningstensor 7 =pò +17, waarin 7 
dat gedeelte van de spanningstensor is dat 
in de evenwichtstoestand verdwijnt. Het 
is merkwaardig dat 7 — 7, voor verdunde 
oplossingen en 7 voor polymeersmelten 
dezelfde vorm hebben. Substitutie van de 
uitdrukkingen voor K, in deze formules 
geeft de moleculaire uitdrukkingen voor 
de spanningstensor (N.B: in de uitdruk- 
king voor Tr voor polymeeroplossingen be- 
tekent < > de gemiddelde waarde in de 
keten-configuratieruimte, maar in de uit- 
drukking voor r voor smelten betekent 

< > de gemiddelde waarde ten opzichte 
van f(u, 0, t), aangezien wij de geringe- 
kromming-benadering hebben gemaakt.) 
Volgens de continuümmechanica kan de 
spanningstensor voor visco-elastische 
vloeistoffen ontwikkeld worden in een 
reeks van de Rivlin en Ericksen tensoren 
(zie referentie 11, blz. 455): 


T= —ô1Y — Bur? + 62 


B Ne edn Ee ae 
Cl DVi (VD-À+ ATD + Een 


òt 
(12) 


De B, -term is de lineaire term, en de 5, - 
en B: -termen zijn van de tweede orde 
(overeenkomstig met de Burnett-verge- 
lijkingen in de kinetische theorie van ver- 
dunde, éénatomige gassen). Als (12) met 
de overeenkomstige moleculaire uitdruk- 
kingen wordt vergeleken, verkrijgen we de 
uitdrukkingen voor 6, , 62, 011, enz: 


mmm 


1. Verdunde oplossingen van Kramers- 
ketens 

(ns = viscositeit van het oplosmiddel; 
coëfficiënten na @, zijn niet exact) 


Brent ze WP Date? 3 


(N? — 1) (1ON® — 1I2N? +35N — 12) 
E 32400N 
ni? a* 


kT GE 


* Deze vergelijkingen werden afgeleid zonder 
de geringe-kromming-benadering (zie referentie 
11, blz. 635). 


t. (11) 


Ri [1 + (y[O1: uu)]3/2 


2. Smelten van Kramers-ketens 


Ì € 
Bi = (zo "5) NnkTN, (15) 


60 “ 90 

B, = En NnkTN? 16 

ZT 600 Nn ) ( ) 

6 ( er |) NakTX2. (17 
mj na n 

E \420 © 1050 ) 


Voor polymeersmelten zijn de coëfficiën- 
ten t/m de derde orde uitgerekend; boven- 
dien is de volledige ‘constitutievergelij- 
king’ (r als een tijdsintegraal over een 
functie van de deformatietensor yl®Ì) ver- 
kregen. De theorie geeft een zeer goede 
beschrijving van de meeste rheologische 
eigenschappen: de viscositeit bij oneindig 
kleine afschuifsnelheid is evenredig met 
M3+B: de viscositeit is een monotoon 
dalende functie van de afschuifsnelheid; 
spanningsrelaxatie, normaalspanningen, 
en andere interessante rheologische ver- 
schijnselen worden eveneens op een be- 
vredigende wijze voorspeld!“ 15. 
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Overeenstemmende toestanden en 
universaliteit bij het kritieke punt 


A. Levelt Sengers 


Samenvatting 

De wet van overeenstemmende toestan- 
den is een vrijwel volmaakte illustratie 
van de fysische werkmethode die een gro- 
te verscheidenheid van verschijnselen in 
uiteenlopende systemen begrij pelijk 
maakt door ze te herleiden tot een enkel, 
en in dit geval uiterst eenvoudig, begin- 
sel: mits in geschikte eenheden gemeten 
is het thermodynamische gedrag van alle 
fluïda hetzelfde. Geformuleerd door 
Van der Waals in 1880, vormt deze wet 
de basis voor onze theoretische en prak- 
tische kennis van de thermofysische 
eigenschappen van vloeistoffen en gassen 
en ligt dus bijvoorbeeld ten grondslag 
aan de hedendaagse chemische industrie. 
Als belangrijke vroege resultaten bespre- 
ken wij hier de theorie van mengsels van 
Van der Waals en de vloeibaarmaking van 
helium door Kamerlingh Onnes. Reeds in 
1881 legde Kamerlingh Onnes verband 
tussen de aard van de intermoleculaire 
krachten en de geldigheid van de wet van 
overeenstemmende toestanden. Een sta- 
tistisch-mechanische rechtvaardiging 
kwam echter pas veel later, rond de jaren 
1940, dank zij het werk van De Boer en 
anderen. Een zo algemeen beginsel als de 
wet van overeenstemmende toestanden 
bleek in de praktijk slechts in benadering 
te gelden. De gevonden afwijkingen wer- 
den, zoals vaak gebeurt in de natuurkun- 
de, gebruikt om het principe te wijzigen 
en te verfijnen. Verscheidene van de mo- 
derne wijzigingen zijn van groot prak- 
tisch belang voor het voorspellen van 
eigenschappen van stoffen die uit gecom- 
pliceerde moleculen bestaan. Maar de 
oudste en theoretisch het meest interes- 
sante uitbreiding was het quanteuze be- 
ginsel van overeenstemmende toestan- 
den, waaraan de naam van De Boer on- 
verbrekelijk is verbonden. 


In de laatste vijftien jaar heeft zich het 
beginsel van universaliteit bij het kritie- 
ke punt ontwikkeld. Dit beginsel schrijft 
gelijkheid van kritiek gedrag toe aan een 
klasse van systemen die behalve de flui- 
da ook magnetische stoffen, binaire le- 
geringen met orde-war overgangen, en 
een aantal theoretische modellen omvat. 
Het wordt wel voorgesteld alsof het be- 
ginsel van universaliteit de moderne 
plaatsvervanger is van de wet van over- 
eenstemmende toestanden. We zullen 
echter duidelijk trachten te maken dat 
dit niet zo is. 
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| Het beginsel van overeenstemmende 


toestanden 

Dit beginsel werd een eeuw geleden door 
Van der Waals geformuleerd in de vol- 
gende vorm: wanneer de druk P, de abso- 
lute temperatuur 7 en het volume V van 
een gas worden gemeten in eenheden van 
de kritieke druk P,, de kritieke tempera- 
tuur 7, en het kritieke volume V,, dan is 
de resulterende toestandsvergelijking 
voor alle gassen dezelfde. Voor de aflei- 
ding maakte Van der Waals gebruik van 
zijn toestandsvergelijking: 
(P+afV?)(V —b)=RT, (1) 
waarin R de gasconstante is, bd een para- 
meter met de dimensie van een volume 
die aangeeft welk deel van het molaire 
volume V is uitgesloten tengevolge van 

de eindige uitgebreidheid van de mole- 
culen, a een parameter met de dimensie 
van energie maal volume die de vermin- 
dering van de druk tengevolge van de 
moleculaire attractie in rekening brengt. 
De geldigheid van de wet van overeen- 
stemmende toestanden hangt echter niet 
af van die van de wet van Van der Waals. 
Als de toestandsvergelijking voldoende is 
gekarakteriseerd door niet meer dan twee 
constanten die van stof tot stof verschil- 
len, en als de vorm van de vergelijking 
niet afhangt van de gebruikte eenheden 
van P, Ten V, dan kunnen deze groothe- 
den dimensieloos worden gemaakt door 
geschikte combinaties van R en de twee 
constanten. In deze ‘gereduceerde’ eenhe- 
den is de vergelijking dan vrij van specifie- 
ke constanten en dus zijn de karakteristie- 
ke eigenschappen zoals de kritieke para- 
meters, de dampspanningskromme en de 
tweede viriaalcoëfficiënt dezelfde voor 
alle stoffen waarop de vergelijking toe- 
pasbaar is. 

Het staat ons dus ook vrij de kritieke 
parameters als reductieparameters te 
gebruiken; doet men dit dan wordt de 
oorspronkelijke formulering van Van der 
Waals hervonden. Het was in de vorige 
eeuw al bekend dat de wet van Van der 
Waals slechts zeer benaderd geldt. Voor 
het beginsel van overeenstemmende toe- 
standen mag men echter een wijdere gel- 
digheid verwachten. Een noodzakelijke, 
maar zeker niet voldoende toets is de 
waarde van het dimensieloze quotiënt 

Z. =P, VIRT, die voor alle stoffen de- 
zelfde dient te zijn. In tabel I ziet men 
typische waarden van dit quotiënt voor 





Mw. dr. A. Levelt Sengers studeerde fysi- 
ca aan de Universiteit van Amsterdam en 
promoveerde (evenals De Boer) bij 
A.A.F.M. Michels; haar proefschrift han- 
delde over de experimentele bepaling van 
de toestandsvergelijking van argon (1958). 
Daarna werkte zij een jaar bij J.O. Hirsch- 
felder aan de Universiteit van Wisconsin, 
onder andere aan een toepassing van het 
door De Boer ingevoerde quanteuze be- 
ginsel van overeenstemmende toestanden 
op de gecondenseerde fasen van de edel- 
gassen. Tot 1963 was zij verbonden aan 
het Van der Waalslaboratorium te Am- 
sterdam. In 1963 werd zij staflid op het 
National Bureau of Standards, waar haar 
belangstelling zich richtte op kritieke ver- 
schijnselen van fluida, dank zij de inspi- 
ratie van M.S. Green. Haar werk omvatte 
de toetsing van schaalwetten en het ex- 
perimentele onderzoek van kritieke ano- 
malieën. Aan het Instituut voor Theore- 
tische Fysica in Amsterdam bestudeerde 
zij in 1974/75 de geschiedenis van de 
Nederlandse fysica in de periode 1870- 
1910. Sinds 1978 leidt zij op het NBS 
een groep die werkt op het gebied van 
theorie, experiment en correlatie van 
thermofysische eigenschappen van gassen 
en hun mengsels. 





een aantal gassen. De ontoereikendheid 
van de wet van Van der Waals blijkt uit 
het feit dat alle experimentele waarden 
aanzienlijk lager zijn dan de theoretische. 
De kritieke waarden van Z, zijn daaren- 
tegen binnen een of twee procent gelijk 
voor de zware edelgassen en de zware 
twee-atomige gassen. Binnen de groep 
van niet-vertakte alkanen zien we een 
duidelijk verloop van Z,. De Z,-waarden 
van stoffen met waterstofbruggen zijn 
veel lager dan die van niet-geassociëerde 
stoffen. Merk ook op dat de Z,-waarden 
van de lichte een- en twee-atomige gas- 
sen hoger zijn dan die van de zwaardere. 
De meeste van de hier genoemde feiten 
waren rond de eeuwwisseling reeds be- 
kend. 


Het beginsel van mechanische 
equivalentie 
Aan de oorspronkelijke formulering van 
de wet van overeenstemmende toestan- 
den lagen geen echt microscopische of 
statistische beschouwingen ten grondslag. 
Het eerste microscopische beeld werd ge- 
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geven door Kamerlingh Onnes in 1881 en 
ging dus aan de statistische mechanica 
van Gibbs vooraf. Het was gebaseerd op 
het principe van mechanische gelijkvor- 
migheid of equivalentie, dat in de kine- 
matica en hydrodynamica zijn kracht al 
had bewezen. Weer wordt een dimensie- 
argument gebruikt, dat nu echter de gede- 
tailleerde beweging van de moleculen van 
systeem A afbeeldt op die van systeem B 
door een geschikte keuze van de schalen 
van lengte, tijd en massa. Er moet wor- 
den ondersteld dat de moleculen in de 
twee systemen gelijkvormige lichamen 
met gelijkvormige massaverdeling zijn 
(niet noodzakelijk sferisch symmetrisch) 
en dat de krachten tussen de moleculen 
in de twee systemen, op een evenredig- 
heidsconstante na, dezelfde functies zijn 
van de geschaalde afstand. In overeen- 
stemmende toestanden verlopen dan de 
individuele bewegingen van de moleculen 
in de twee systemen ‘mechanisch gelijk- 
vormig’, zodat niet alleen de drukken, 
maar ook de andere eigenschappen, zo- 
als de oppervlaktespanningen en de visco- 
siteiten, overeenstemmen. 


Men kan tegen deze beschouwing opwer- 
pen dat slechts het duiveltje van Maxwell 
de moleculen in twee systemen gelijkvor- 
mig kan laten dansen, zodat de eis van 
mechanische gelijkvormigheid te streng 
moet zijn. De statistische mechanica zou 
hier uitkomst brengen, maar dat gebeurde 
pas veel later, namelijk rond 1940. De 
kracht van het argument van Kamerlingh 
Onnes lag in het verschuiven van de di- 
mensiebeschouwing naar het microsco- 
pische vlak; de gevonden afwijkingen 
kunnen nu in verband worden gebracht 
met bepaalde niet-gelijkvormigheden van 
de moleculaire wisselwerking. Deze mo- 
gelijkheid werd reeds gesuggereerd door 
de manier waarop wij de stoffen in tabel 
Il rangschikten. Wij komen op deze zaak 
terug. | 


Vroege successen: de theorie van 
mengsels 

Binaire mengsels hebben een extra ther- 
modynamische vrijheidsgraad vergeleken 
met een-componentige systemen. Het 
gevolg is een verbluffende verrijking van 
het fasengedrag: behalve de fasenscheiding 
gas-vloeistof kunnen ook scheidingen 
vloeistof-vloeistof en gas-gas optreden. De 
wet van overeenstemmende toestanden 
gaf Van der Waals de sleutel tot de klassi- 
ficatie en beschrijving van het fasenge- 
drag van mengsels. Hij postuleerde dat 
een binair mengsel van constante meng- 
verhouding x in overeenstemmende toe- 
stand was met de twee pure componen- 
ten, door ze alle drie dezelfde toestands- 
vergelijking met twee specifieke parame- 
ters te beschrijven. Hij onderstelde dat 

de twee parameters van het mengsel func- 
ties waren van die van de twee pure com- 
ponenten en van de mengverhouding. De- 
ze eenvoudige veronderstellingen, gecom- 
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Fig. 1. Fasenevenwicht in mengsels, vol- 
gens computerberekeningen met het mo- 
del van Van der Waals!. Onder ziet men 
een P-T-projectie met de dampspannings- 
krommen van de pure componenten (1, 
2) en een driefasenevenwicht vloeistof: 
vloeistof-gas (L, Ly G) dat heel dicht bij 
de dampspanning van 1 ligt (—). Eén 
kritieke lijn (—-) loopt van het kritieke 
punt van 1 naar het kritieke eindpunt op 
L, LL, G. De andere loopt van het kritie- 
ke punt van 2 weg, tenslotte naar hoge 
drukken. Boven zijn T-x-diagrammen 
weergegeven voor hetzelfde systeem bij 
acht waarden van de gereduceerde druk. 


bineerd met de theorie-van heterogene 
evenwichten van Gibbs, leidden tot vrij- 
wel alle typen van fasenscheiding die nu 
bekend zijn. De hoeveelheid algebra die 
Van der Waals en vooral Van Laar rond 
de eeuwwisseling hebben verzet om tot 
deze klassificatie te komen is onvoor- 





Pe 


stelbaar. 

Onlangs hebben Van Konynenburg en 
Scott het programma van Van der Waals 
uitgevoerd op de computer! . In figuur 1 
ziet men het resultaat voor een mengsel 
uit klasse drie: de kritieke punten van 


\_de pure componenten zijn niet meer door 


Tabel l. Waarden van het kritieke quotient Z. =P. V./RT 





Edelgasen Niet-vertakte alkanen 
licht 3He 0,307 CH4 0,288 
“He 0,303 C‚ He 0,285 

Ne 0,308 C3 Hs 0,281 

zwaar Ar 0,290 C4H‚o 0,274 
Kr 0,291 CsH‚» 0,262 
Xe 0,287 

Tweeatomig ‘Geassocieerd” 

licht H‚, 0,301 HO 0,229 

zwaar N, 0,289 CH3 OH 0,224 
O, 0,288 C‚Hs OH 0,248 
F, 0,288 | 


Van der Waalsvergelijking 
0,375 
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Fig. 2. Dit is een plaatje uit het werk van 
Van Laar in 1904/05. Het stelt een soort- 
gelijk mengsel voor als dat van figuur 1, 
maar nu in de ruimte P-T-x. 


een continue kritieke lijn verbonden en er 
treedt bovendien vloeistof-vloeistoffasen- 
scheiding op met een tweede kritieke lijn. 
In figuur 2 vindt men een soortgelijk sy- 
steem, maar nu door de ogen en van de 
hand van Van Laar (1904). Het vergt enig 
ruimtelijk inzicht om de computerplaatjes 
in Van Laars kunstwerk terug te vinden. 
Het idee van Van der Waals, overeenstem- 
mende toestanden en concentratie-afhan- 
kelijke parameters, vormt tot op de huidi- 
ge dag het werktuig dat de chemisch in- 
genieur hanteert als hij mengseleigenschap- 
pen moet voorspellen. Dit werktuig is nu 
bekend onder de naam “one-fluid model” 
of “conformal solution theory”. 


Vroege successen: het vloeibaar maken 
van helium 

Kamerlingh Onnes begon in de jaren 
tachtig van de vorige eeuw zijn experimen- 
tele werk in Leiden met het expliciete 
doel de ‘permanente’ gassen vloeibaar te 
maken. In die tijd was het cascadeproces 
in gebruik: het trapsgewijs afkoelen met 
behulp van een serie gassen van afdalende 
kritieke temperatuur. Het gebruik van 
methylchloride en ethyleen leidt tot een 
produktieproces van vloeibare zuurstof 
dat door Kamerlingh Onnes in 1892 in 
werking werd gesteld. Met vloeibare zuur- 
stof kan men echter niet beneden het 
kritieke punt van waterstof komen. De 
joule-thomsontemperatuur, beneden wel- 
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ke een gas door expansie afkoelt, ligt 
echter veel hoger dan de kritieke en kon 
voor waterstof wel met vloeibare zuurstof 
worden bereikt. De vloeibaarmaking van 
waterstof moest echter wachten op de 
uitvinding van het dewarvat of de ther- 
mosfles (1892). Kamerlingh Onnes ver- 
loor de race tegen Dewar, die in 1898 
waterstof vloeibaar maakte door het cas- 
cadeproces met de joule-thomsonexpansie 
te combineren. Dat Kamerlingh Onnes 
echter als winnaar tevoorschijn kwam in 
de vloeibaarmaking van helium (1908) 
dankte hij aan de wet van overeenstem- 
mende toestanden. Door de eigenschap- 
pen van waterstof ver boven zijn kritieke 
punt te meten bij temperaturen bereik- 
baar met vloeibare zuurstof, en later die 
van helium bij temperaturen bereikbaar 
met vloeibare waterstof, kon hij, door ge- 
bruik te maken van de wet van overeen- 
stemmende toestanden met gassen waar- 
van het gedrag reeds bekend was, een 
schatting maken van de kritieke tempera- 
tuur, de verdampingswarmte en andere 
parameters die hij nodig had voor het ont- 
werpen van zijn compressoren, pompen 
en warmtewisselaars. Zo schatte hij een 
jaar voor de vloeibaarmaking van helium 
de kritieke temperatuur op 5,3 K, slechts 
0,1 K boven de waarde die wij nu ken- 
nen. Bij gelegenheid van het afscheid van 
Van der Waals van zijn leerlingen (11 juli 
1908), twee dagen na de vloeibaarma- 


king, schreef Kamerlingh Onneseenen- 


thousiaste brief naar zijn vriend waarin 
hij het werk van die historische dag be- 
schreef. De ochtend van de 9° was zijn 
staf om kwart voor zes begonnen met 
vloeibare waterstof te bereiden. Om een 


uur ’s middags was er 22 liter klaar waar- 
mee het apparaat werd gekoeld. “ Toen de 
laatste fles vloeibare waterstof aangekop- 
peld werd was er nog geen helium te 

zien. De toestel was echter bij een van 

de expansies volgelopen. Een weinig la- 
ter werd dan ook het oppervlak gezien’ 

… “Bij wat ik u nu schreef is uw theorie 
voortdurend leiddraad geweest. De kriti- 
sche temperatuur van een permanent gas 
uit de isothermen te berekenen herinnert 
nog eens op nieuwe wijze aan uw disserta- 
tie. De berekeningen zijn geheel volgens 
de wet van overeenstemmende toestan- 
den geregeld. Daardoor geleid schatte ik — 
al heb ik dat nu niet zo op papier gezet 

— 20 liter (waterstof) nodig te hebben. 
Had ik een paar liter minder geschat dan 
had de proef geen resultaat opgeleverd, had 
ik veel meer geschat, dan had ik het niet 
raadzaam geacht ze met mijn hulpmidde- 
len te ondernemen. 


De vloeibaarmaking van helium was het 
hoogtepunt maar tevens het eindpunt van 
een periode van bloei van de klassieke 
thermodynamica. De belangstelling van 
fysici voor dit werkgebied zou pas in de 
jaren dertig herleven. De verbindingsdraad 
die loopt van de vloeibaarmaking van he- 
lium naar het werk van Uhlenbeck en De 
Boer is die van de quantumeffecten die 
zich manifesteerden in het nieuwe, cryo- 
gene werkgebied. 


Empirische schaalfactoren 

Het toetsen van het beginsel van overeen- 
stemmende toestanden kan geschieden 
door voor elke stof de betreffende fysi- 
sche variabelen te reduceren of wel met 
de kritieke parameters, of wel met R en. 
de twee parameters die de wisselwerking 
karakteriseren. Gewoonlijk zijn deze re- 
ductieparameters minder precies bekend 
dan de te reduceren grootheden zelf. Om 
de experimentele nauwkeurigheid zoveel 
mogelijk te handhaven stelde Kamerlingh 
Onnes voor om de eigenschappen van gas- 
sen op elkaar af te beelden door het ge- 
bruik van twee empirische schaalfactoren 
per paar. Deze methode werd door mij in 
1960 opnieuw voorgesteld bij mijn toet- 
sing van overeenstemmende toestanden 
in argon en xenon? ; de methode werd op 
elegante wijze gemodificeerd door Been- 
akker en Van Eynsbergen in 1968, en cul- 
mineerde in twee artikelen van Kestin en 
medewerkers®»*. Dit laatste: werk toont 
aan dat met twee empirische schaalfacto- 
ren per stof overeenstemmende toestan- 
den simultaan gelden voor de tweede vi- 
riaalcoëfficiënt, viscositeit en diffusieco- 
efficiënt van het verdunde gas, waarin 
interacties tussen meer dan twee deeltjes 
kunnen worden verwaarloosd. De bestu- 
deerde gassen zijn de edelgassen, eenvou- 
dige niet-polaire twee-atomige gassen en 
achttien van hun binaire mengsels over 
een zeer groot temperatuurgebied. In fi- 
guur 3 vindt men de geschaalde krommen 


voor de tweede viriaalcoëfficiënt en de 
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Fig. 3. Overeenstemmende toestanden 
met empirische schaalfactoren voor edel- 
gassen en hun mengsels in het verdunde 
gas®. Het beginsel geldt nauwkeurig en 
simultaan voor de tweede viriaalcoëffici- 
ent (boven) en de viscositeit (beneden) in 
een zeer groot temperatuurgebied. De ho- 
rizontale lijnstukken duiden de meetge- 
bieden aan voor de betrokken gassen. 
Grootheden worden na reductie weerge- 
geven door het symbool van een sterretje 
te voorzien. 


viscositeit. Zolang men zich beperkt tot 
de verdunde gasfase, geldt het principe 
van overeenstemmende toestanden met 
een nauwkeurigheid die die van de beste 
metingen benadert (0,1%). De zo verkre- 
gen empirische schaalfactoren wijken 
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echter af van degene die men uit de kri- 
tieke parameters berekent, zoals mij even- 
eens was opgevallen toen ik de wet van 
overeenstemmende toestanden onder- 
zocht voor argon en xenon?. Het werk 
van Kestin et al. documenteert afwijkin- 
gen van overeenstemmende toestanden 
voor kritieke parameters, tripelpuntswaar- 
den en andere eigenschappen van verdich- 
te systemen wanneer men de schaalfacto- 
ren van de verdunde gasfase gebruikt voor 
reductie. Wij komen hierop terug. 


Overeenstemmende toestanden en 
intermoleculaire krachten 

Pas in de jaren dertig werd het thema van 
Kamerlingh Onnes, overeenstemmende 
toestanden door gepostuleerde gelijkvor- 
migheid van intermoleculaire krachten, 
weer opgenomen. Door Lennard-Jones 


was in 1924 een vorm voor de molecu- 
laire wisselwerking ingevoerd die nu naar 
hem is genoemd: 

V(r)= 4e{(rlo) 2 —(r/o)P }. (2) 
De intermoleculaire afstand is aangegeven 
door r, e stelt voor de wisselwerkingsener- 
gie bij het potentiaalminimum, o de af- 
stand r waarvoor de potentiaal verdwijnt. 
Vergelijking (2) heeft een zekere theore- 
tische rechtvaardiging, beschrijft de wissel- 
werking van sferische moleculen redelijk 
wel, en voldoet aan het criterium voor 
mechanische gelijkvormigheid. In 1938 
toonden De Boer en Michels door een di- 
mensie-argument aan dat de viriaalcoëffi- 
ciënten in de dichtheidontwikkeling van 
de druk van een gas overeenstemmende 
toestanden vertonen voor een potentiaal 
van deze vorm5. Door Pitzer (1939), Gug- 
genheim (1945) en De Boeré (1948) wer- 
den de condities voor overeenstemmende 
toestanden gepreciseerd door uit te gaan 
van de statistisch-mechanische toestand- 
som. Mits de intermoleculaire krachten 
gelijkvormig zijn en paarsgewijs additief, 
de inwendige vrijheidsgraden onafhanke- 
lijk zijn van de intermoleculaire afstanden, 
en klassieke statistische mechanica kan 
worden gebruikt, geldt de wet van over- 
eenstemmende toestanden voor de druk 
en voor alle evenwichtseigenschappen die 
niet afhangen van de inwendige vrijheids- 
graden. De Boer en Van Kranendonk 
breidden het principe uit tot de transport- 
eigenschap viscositeit (1948). 


De afwijkingen van het principe die expe- 
rimenteel zijn gevonden moeten dus wor- 
den toegeschreven aan het falen van een 
of meer van de bovenstaande condities. 
Voor de door Kestin et al. bestudeerde 
systemen staat het principe van gelijkvor- 
migheid van krachten minder onder ver- 
denking dan dat van paarsgewijze additi- 
viteit van krachten. De geldigheid van 
klassieke statistische mechanica, tenslot- 
te, is aan de orde gesteld in het werk van 
De Boer dat we nu bespreken. 


Het principe van overeenstemmende 
toestanden in de quantummechanica 
Reeds in het begin van deze eeuw was be- 
kend dat de lichte gassen helium en water- 
stof opmerkelijke afwijkingen van het be- 
ginsel van overeenstemmende toestanden 
vertoonden vergeleken met zwaardere 
een- en twee-atomige gassen. In 1912 op- 
perden Kamerlingh Onnes en Keesom in 
voetnoot 517 van hun lijvige artikel 

“Die Zustandsgleichung” het idee dat de 
hypothese van Planck die leidde tot quan- 
tisering van de faseruimte van oscillatoren 
in gebiedjes van de maat 4, ook implica- 
ties moest hebben voor de toestandsverge- 
lijking. Dit idee werd weer opgevat door 
Byk in 1921; gebruikmakend van het be- 
ginsel van mechanische equivalentie voer- 
de hij het dimensieloze werkingsquantum 
hIKT HU? MU? (V/N)H? in als derde 
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Tabel II. 
Parameter van De Boer, A* =ho V/me 


3 He 2,94 3 He — “He 2 
4 He 2,55 4 He — Ne 1,38 
Ne 0,53 “He — Ar 0,88 
Ar 0,17 *He — Kr 0,76 
Kr 0,09 “He — Xe 0,67 
Xe 0,06 Ne — Ar 0,34 
H. 1,73 
- 1,22 
N, 0,23 


parameter in de dimensieloze vrije ener- 
gie (m is de moleculaire massa, N het ge- 
tal van Avogadro). Wanneer de waarde 
van deze derde parameter belangrijk 
wordt, kunnen afwijkingen worden ver- 
wacht van het beginsel van overeenstem- 
mende toestanden. Dit is het geval voor 
helium en waterstof met hun lage waar- 
den van 7. Men voelt zich een beetje on- 
behaaglijk over dit huwelijk tussen het 
werkingsquantum uit de golfmechanica 
en het klassieke principe van mechanisch- 
equivalente beweging en over het feit dat 
de reductieparameters zelf zijn beïnvloed 
door quantumeffecten. Niettemin was 
Byks inzicht dat / expliciet moest optre- 
den en dan alsnog dimensieloos moest 
worden gemaakt, juist. In 1936 deden 
Uhlenbeck en Beth de eerste quantum- 
statistische berekeningen van de tweede 
viriaalcoëfficiënt voor harde bollen als 
een ontwikkeling in A2. De Boer en 
Michels formuleerden dit probleem in 
1938 voor de lennard-jonespotentiaal 
(2). Zij gebruikten de potentiaalparame- 
ters e en o als reductieparameters en 
voerden als gereduceerd werkingsquan- 
tum in de grootheid: 


A* = Alo =hl(ome), 


die nu in de Engelse literatuur bekend 
staat als de “De Boer parameter’. De fy- 
sische interpretatie die De Boer aan deze 
parameter gaf is dat A de de-brogliegolf- 
lengte voorstelt van de gemiddelde rela- 
tieve beweging van twee moleculen als de 
temperatuur gelijk is aan e/k. In tabel II 
geven wij typische waarden voor A*, 
voornamelijk ontleend aan referentie 4. 
Men ziet dat de parameter van De Boer 
waarden bereikt die groter zijn dan 1 
voor “He, “He, H‚, Ds», en voor een aan- 
tal mengsels met deze gassen als compo- 
nent. 


3) 


In de jaren veertig onderzocht De Boer 
systematisch de invloed van quantumef- 
fecten op het thermodynamische gedrag. 
In 1948 formuleerde hij het quanteuze 
beginsel van overeenstemmende toestan- 
dené, hetgeen inhoudt dat de met mole- 
culaire parameters gereduceerde thermo- 
dynamische eigenschappen van gassen bij 
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Fig. 4. Met moleculaire parameters gere- 
duceerde dampspanningskrommen van 
lichte en zware gassen, volgens De Boer 
en Lunbeck®. De gestippelde lijnen ver- 
binden de tripelpunten en de kritieke 
punten. De afwijkingen van overeenstem- 
mende toestanden zijn bijzonder groot 
en dienen te worden toegeschreven aan 
quantumeffecten. Grootheden worden na 
reductie weergegeven door het symbool 
van een sterretje te voorzien. 


corresponderende waarden van de onaf- 
hankelijke variabelen nog functies moeten 
zijn van de parameter A*. De Boer en 
Blaisse’ toonden aan dat de gereduceerde 
eigenschappen van de vaste edelgassen 
kunnen gerangschikt worden naar op- 
klimmende waarden van A*. De Boer en 
Lunbeck?® ordenden een groot aantal 
eigenschappen van de lichtere gassen als 
functies van A*. Hoe groot de afwijkin- 
gen van overeenstemmende toestanden 
zijn, kan men voor de dampspannings- 
kromme zien in figuur 4, die is overgeno- 
men uit hun artikel. In het voetspoor 
van Kamerlingh Onnes voorspelden Lun- 
beck en De Boer de kritieke parameters 
en de dampspanning van * He voor deze 
waren gemeten. Zij maakten ook schat- 
tingen van de kritieke parameters en de 
tripelpuntsparameters van de waterstof- 


isotopen HD, HT, DT en T,. 


Voor de beschrijving van de vloeistoffase 
berekenden De Boer en Lunbeck?® de 

eerste quantumcorrectie in het celmodel 
van Lennard-Jones en Devonshire in een 


ontwikkeling naar machten van n°. Hurst 


en ik hebben later de energieniveaus in 
het celmodel rechtstreeks berekend 
(1960). Uit onze berekeningen volgde dat 
zelfs in vloeibaar argon de quantumeffec- 
ten niét geheel verwaarloosbaar zijn, en 


ee 
Pand 
md 
pr eed 
mad 
amd 
men, 
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dus mogelijk bijdragen aan de met de 
dichtheid groeiende afwijkingen van over- 
eenstemmende toestanden voor de zware 
edelgassen? >. 

Het quanteuze beginsel van overeenstem- 
mende toestanden heeft toepassing ge- 
vonden in het beschrijven en verklaren 
van verschillen in thermofysische eigen- 
schappen van de isotopen van lichte gas- 
sen, waaraan Beenakker en Knaap veel 
werk deden in de jaren zestig, en in het | 
berekenen van de nulpuntsenergie van 
moleculaire vaste stoffen. Een recente 
toepassing is de voorspelling van het fa- 
sengedrag van atomaire waterstof en 
deuterium door Nosanov, die hiertoe 

de parameter van De Boer invoerde als 

een extra thermodynamische vrijheids- 
graad. Spectaculaire metingen aan spin- 
gepolariseerde atomaire waterstof gebeu- 

ren door Silvera e.a. in het Amsterdams 
Natuurkundig laboratorium van Silve- — 

rase, 


Het is vooral de fysische aanpak van De 
Boer geweest die door velen is nagevolgd. 
Ofschoon aan een van de voorwaarden 
voor de geldigheid van overeenstemmende 
toestanden niet was voldaan, verwierp hij 
het beginsel niet, maar breidde het uit 
met behulp van een derde parameter. De 
schending van andere voorwaarden kan op 
soortgelijke wijze worden behandeld. Zo 
voerde Rowlinson (1951) het dipoolmo- 
ment in als een derde parameter, terwijl 
Pitzer voor niet-sferische moleculen de 
“acentric factor” invoerde, die grote po- 
polariteit heeft gevonden bij de chemi- 
sche ingenieurs. De strengheid van De 
Boers formulering wordt echter zelden 
gehandhaafd. In de harde praktijk van de 
ingenieurswereld verwateren de molecu- 
laire parameters snel tot empirische aan- 
passingsconstanten, en de verworven ver- 
beterde beschrijving leidt niet tot verdie- 
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ping van de kennis van de moleculaire 
wisselwerking. 


Universaliteit bij het kritieke punt 

Het kritieke punt van een gas is een mar- 
ginaal stabiel punt waarbij het gas zich 
voor het eerst in fasen kan verdelen. In 
het kritieke gebied is het gas uiterst sa- 
mendrukbaar en de fasensplitsing wordt 
aangekondigd door fluctuaties in de 
dichtheid die zich over zeer vele molecu- 
laire afstanden uitstrekken. Een maat 
voor de grootte van de fluctuaties is de 
correlatielengte, die bij het kritieke punt 
divergeert. In de klassieke beschrijving, 
waarvan de wet van Van der Waals het : 
prototype is, worden deze fluctuaties ver- 
waarloosd, en die beschrijving moet dus 
_falen bij het kritieke punt. Zo voorspelt 
de klassieke theorie dat de top van de co- 
existentiekromme parabolisch is, terwijl 
hij in werkelijkheid een meer kubische 
vorm heeft en dat de soortelijke warmte 
G, eindig blijft, terwijl hij in werkelijkheid 
divergeert. 


In de laatste vijftien jaar is de studie van 
kritieke fluctuaties enorm ontwikkeld. De 
klassieke beschrijving is vervangen door de 
schaalwetten, ingevoerd door Widom 
(1965) en plausibel gemaakt door Kada- 
noff (1967), en is vervolgens op een meer 
solide fundering gezet door Wilson (1971) 
die de methode van de renormalisatie- 
groep invoerde, die op Nederlandse bo- 
dem bloeit dank zij het werk van Van 
Leeuwen en vele medewerkers. 


Het anomale deel van het kritieke gedrag 
treedt naar voren wanneer de correlatie- 
lengte zeer groot wordt. Men kan zich nu, 
met Kadanoff, afvragen wat er gebeurt 
met de beschrijving van een systeem on- 
der vergroting van de schaal waarop het 
wordt bekeken. Beschouwt men het eerst 
op de schaal van de gemiddelde molecu- 
laire afstand a, en dan op een nieuwe 
schaal La, dan is de correlatielengte in de 
nieuwe schaal relatief L maal zo klein als 
op de oude; het nieuwe systeem is dus 
verder van zijn kritieke punt dan het 
oude. Ook zijn de krachten tussen de 
groepjes van L* moleculen in het nieuwe 
systeem (d = dimensie) niet dezelfde als 
die tussen de moleculen in het oude sy- 
steem. Men kan zich voorstellen dat men 
het nieuwe systeem uit het oude zou kun- 
nen verkrijgen niet alleen door de schaal- 
verandering, maar ook door het oude sy- 
steem van zijn kritieke punt te verwijde- 
ren. Dit laatste kan gebeuren door ieder 
van de twee onafhankelijke variabelen, ge- 
meten vanaf het kritieke punt, te verme- 
nigvuldigen met een macht van de schaal- 
factor L. Door het anomale deel van de 
vrije energie voor deze twee transforma- 
ties te identificeren, vindt men dat het 
aan een gegeneraliseerd beginsel van ho- 
mogeniteit moet voldoen. Zoals altijd het 
geval is bij homogene functies, kan men 
dan door de variabelen op geschikte wijze 
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te combineren één onafhankelijke varia- 
bele elimineren. Zo verkrijgt men bijvoor- 
beeld voor de toestandsvergelijking Wi- 
doms schaalwet: 


Au*/{D Ap*| Ap*| ò-! } = g(x), 
X=Xo AT*/| Ap*|UÊ. (3) 


Het met kritieke parameters gereduceerde 
chemische potentiaalverschil Au* van een 
gas wordt geschaald met een macht Ö van 
het gereduceerde dichtheidsverschil 

Ap* = (p — p‚)/p‚; het is dan een functie 
g(x) van één variabele x, die evenredig is 
met het gereduceerd témperatuurverschil 
AT* =(T — TT, geschaald met een 
macht 1/9 van Ap *. De schaaltheorie 
voorspelt dat de kritieke exponenten Ö 

en B, en de schaalfunctie t(Xx) dezelfde 
zijn voor alle gassen: waarden van Ò en 6, 
en de functie g(x) kunnen worden bere- 
kend met de methode van de renormalisa- 
tiegroep. Alleen de twee constanten D en 
Xo kunnen nog van stof tot stof verschil- 
len. 


Het succes van de schaalhypothese blijkt 
uit figuur 5. Hier ziet men de volgens (3) 
geschaalde chemische potentiaal als func- 
tie van de variabele x voor vijf gassen, 
waaronder een edelgas, Xe, een gas met 
niet-sferische moleculen, CO, , twee 
quantumgassen, * He en * He, en een geas- 
socieerde stof, H‚ O. Ofschoon deze vijf 
stoffen niet aan de wet van overeenstem- 
mende toestanden voldoen (tabel 1), val- 
len de volgens (3) geschaalde chemische 
potentialen in een gebied, 0,01 < 

(T — TT, < 0,1, rond het kritieke punt 
op een universele curve!! 


Het beginsel van universaliteit bij het kri- 
tieke punt reikt echter nog veel verder. 


Fig. 5. Het geschaalde chemische-poten- 
tiaalverschil Au*/(D Ap*)| Ap*|Ô-* van 
vijf gassen in het kritieke gebied als func- 
tie van het geschaalde temperatuursver- 
schil x. Voor elk gas zijn twee vrije para- 
meters, D en xo, gebruikt**. De wet van 
overeenstemmende toestanden geldt niet 
voor deze vijf gassen. 





Dat men het systeem op een grovere 
schaal dan de moleculaire mag beschou- 
wen is een gevolg van het feit dat de cor- 
relatielengte groot is. Dan mag ook de 
gedetailleerde vorm van de moleculaire 
wisselwerking niets uitmaken, noch het 
feit of de moleculen zich al dan niet op 
roosterplaatsen bevinden. Het isingmo- 
del, bijvoorbeeld, is een theoretisch mo- 
del voor de ferromagnetische fasenover- 
gang, waarbij magnetische spins, die om- 
hoog of omlaag kunnen wijzen, op een 
rooster zijn geplaatst en slechts in wissel- 
werking staan met hun naaste buren. 
Door ‘spin omhoog’ te vervangen door 
‘bezet roosterpunt’ en ‘spin omlaag’ door 
‘lege plaats’ verkrijgt men een kunstma- 
tig model voor de gas-vloeistofovergang, 
het roostergas; ook is een transcriptie mo- 
gelijk naar een model voor orde-war-over- 
gangen in binaire legeringen. De universa- 
liteitshypothese zegt dat alle systemen 
die dezelfde ruimtedimensie d (hier 3) en 
spindimensie n (hier 1) hebben, in het 
kritieke gebied kunnen worden omschre- 
ven met dezelfde schaalfunctie g( en 
dezelfde kritieke exponenten Ge: Kan 
het gedrag van ferromagneten n. 1, 
binaire legeringen, het isingmodel «… alle 
fluida worden beschreven met de kromme 
van figuur 5? Nu moeten we voorzichtig 
zijn. Het beginsel van universaliteit doet 
slechts uitspraken over het thermodyna- 
mische gedrag dicht bij het kritieke punt, 
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waar de correlatielengte groot is. Al kan 
de theorie nu ook hogere termen in een 
ontwikkeling rond het kritieke punt voor- 
spellen, elke nieuwe term brengt een nieu- 
we ongespecificeerde coëfficiënt mee, zo- 
dat het universele karakter snel verloren 
gaat. De gegevens gebruikt in figuur 5 
gaan, zoals we ons nu realiseren, het 
asymptotische gebied (| 7 — T.I/T, < 
109) reeds ver te buiten en de universele 
kromme is ten onrechte verkregen door 
voor de kritieke exponenten waarden te 
gebruiken die van de theoretische circa 
10% verschillen. Gezien het succes van de 
empirische schaling lopen de niet-univer- 
sele correcties op het schaalgedrag voor 
deze vijf gassen kennelijk niet te ver uit- 
een, maar er is geen enkele reden om te 
verwachten dat dit ook voor andere leden 
van deze universaliteitsklasse zo zou zijn. 


Het beginsel van overeenstemmende toe- 
standen geldt onder de voorwaarde dat 
krachten tussen de moleculen gelijkvor- 
mig en additief zijn, en quantumeffecten 
kunnen worden verwaarloosd; het strekt 
zich dan uit over de gehele fasenruimte 
van de betrokken fluïda. Het beginsel van 
universaliteit geldt ongeacht de aard van 
de krachten (mits kort van dracht) of de 
aanwezigheid van quantumeffecten; het 
geldt niet alleen voor fluïda maar ook 
voor alle andere systemen in dezelfde uni- 
versaliteitsklasse, zoals ferromagneten met 
n= 1, binaire legeringen met orde-war-’ 
overgangen en het isingmodel; maar het 
beperkt zich tot de beschrijving van het 
anomale thermodynamische gedrag in een 
klein gebied rond het kritieke punt en 
kan daarom het beginsel van overeen- 
stemmende toestanden niet vervangen. 


Prof. P.H.E. Meijer, dr. H.P.R. Frederikse 
en prof. P.J. Brussaard verdienen dank 
voor hun opbouwende kritiek. 
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Theorie van 


de grenslaag bij fasescheiding 


J.M.J. van Leeuwen 


Bij een fase-overgang kunnen eigenschappen sprongsgewijs van de ene naar de andere 
iase veranderen. Zo is bij de damp-vloeistof overgang een verandering van de dichtheid 
met een factor 1000 niet ongewoon. Wanneer de twee coëxisterende fasen door een 
uitwendige oorzaak (bijvoorbeeld de zwaartekracht) ruimtelijk gescheiden worden, ver- 
anderen die eigenschappen in het grensvlak macroscopisch gezien discontinu, maar op 
microscopische schaal is er een grenslaag waarin de eigenschappen continu veranderen. 
Het profiel van de grenslaag is van praktisch belang voor ontmenging van binaire lege- 
ringen en vloeistofmengsels. Theoretisch wordt de grenslaag erg interessant in de buurt 
van een kritiek punt waar het onderscheid tussen de fasen verdwijnt en de grenslaag 


een zeer grote dikte kan bereiken. 


Er zijn verschillende systemen zoals damp-vloeistof, binaire legeringen en magnetische 
systemen waar de grenslaag op analoge wijze tot stand komt. De meest representatieve 
parameter voor zo’n grenslaag is de orde-parameter; in bovengenoemde voorbeelden 
respectievelijk de dichtheid, de concentratie en de magnetisatie. 


In dit artikel wordt aandacht gegeven aan de theorie van het profiel van de orde- 
parameter in de grenslaag en in het bijzonder aan de vraag of de fasescheidende kracht 
(zwaartekracht) een belangrijke invloed heeft op dit profiel. 


Klassieke Theorie 

De klassieke grenslaagtheorie gaat via Cahn 
en Hilliard! terug naar Van der Waals?. 
We zullen deze theorie schetsen in de con- 
text van het Ising-model. Weliswaar is dit 
een vrij schematisch model voor bovenge- 
noemde systemen maar het bevat de essen- 
tie van het probleem. Het Ising-model is 
gedefinieerd als een rooster waarop spins 
s; zitten, die slechts twee standen s; = 

+1 kunnen innemen. De energie E van 
een spinconfiguratie wordt gegeven door 
de uitdrukking: 


(1) 


De eerste term bevat de wisselwerking tus- 
sen een nabuurpaar van roosterplaatsen í 
en j. Als de koppelingsconstante J > 0, is 
het energetisch voordelig om naburige 
spins hetzelfde teken te geven, hetgeen 
tot ferromagnetisme kan leiden. De 
tweede term geeft de invloed van een uit- 
wendig magneetveld H; op een spin s;. De 
kans op een spinconfiguratie is evenredig 
met exp —E/KT. 


E=-J X SiS; — xH;s;. 
<ij> Ï 


Als H; > O zal de magnetisatie m; = Ss; 
naar een positieve waarde streven en bij 
H; < O naar een negatieve waarde. Wan- 
neer de temperatuur 7 beneden een kri- 
tieke waarde 7, daalt blijft m een van nul 
verschillende waarde houden als 4 — 0: 
de spontane magnetisatie, die het teken 
heeft van H. De twee fasen die bij H = O en 
T<T, coëxisteren hebben een tegenge- 
stelde (spontane) magnetisatie Mp. 


De fasen kunnen gescheiden worden door 


het veld H; van de plaats i te laten afhan- 
gen, bijvoorbeeld als H; = gz; waarbij z; 
de hoogte van roosterplaats i is. Met op- 
zet is de letter g voor de evenredigheids- 
constante gekozen omdat hij dezelfde rol 
speelt als de zwaartekrachtsversnelling bij 
vloeistof-damp evenwicht. In zo’n situatie 
zal het gebied met z; > 0 naar een magne- 
tisatie Mp tenderen en z; $ O naar —Msp. 
Rond het vlak z; = 0 ontstaat een grens- 
laag waar een geleidelijke overgang van 
Msp naar —Msp plaats vindt. | 
De eenvoudigste theorie hiervoor is de ge- 
middelde veldbenadering (mean field: MF- 
benadering). Deze houdt in dat de magne- 
tisatie m; niet alleen door H; bewerkt 
wordt maar ook door de aanwezige 
magnetisatie m;‚‚ op de buurplaatsen 

i + s van Â. In formule: 


(2) 


Voor een homogeenveld H; = H leidt dit 


m;= tgh (H; a MJ) KT. 


Fig. 1. De magnetisatie m als functie van 
het veld H voor T<T, volgens de MF- 
benadering. 
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versiteit van Amsterdam en promoveerde 
in 1962 bij prof. dr. J. de Boer op dia- - 
gramtechnieken in de statistische mecha- 
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van de Katholieke Universiteit te Nijmegen 
en werkte aan kinetische theorie, elektrisch 
geleidingsvermogen, He-mengsels en kri- 
tische verschijnselen. Hij was in die peri- 
ode ook post-doctoraal research fellow 
aan Cornell University in Ithaca (N.Y.) 
en aan het National Bureau of Standards 
in Washington (D.C.), visiting associate 
professor aan de Rockefeller University 
in New York en associate professor aan 
de Temple University in Philadelphia. In 
1979 werd hij benoemd tot hoogleraar in 
Delft. In de loop der jaren is zijn onder- 
zoeksterrein verschoven van de theorie 
van correlaties in vloeistoffen naar die 
van kritieke verschijnselen en in het bij- 
zonder naar de renormalisatietheorie 
voor spinsystemen. 


tot een homogene magnetisatie m; = m(H) 
en deze functie is geschetst in figuur 1. De 
magnetische toestandsvergelijking m (H) 
vertoont veel overeenkomst met de toe- 
standsvergelijking van Van der Waals voor 
vloeistoffen. In het bijzonder de (gestip- 
pelde) slinger die voor T <7, optreedt is 
analoog. Deze slinger maakt dat er bij 

H = O twee van nul verschillende waarden 
voor m overblijven als H — 0 en verklaart 
dus de spontane magnetisatie. 


Wanneer H; van de positie van í afhangt is 
(2) moeilijker op te lossen. Voor H; =8z; 
en g klein is bij niet al te lage temperatu- 
ren te verwachten dat m; geleidelijk van 
z; afhangt en we benaderen: 


EZ miis=CM; ta? 9? m;/dz; +… (3) 
S 


Fig. 2. Het profiel van de magnetisatie m 
in de grenslaag. 








144 


waarbij c het aantal naasteburen is en a 
de roosterafstand. Substitutie van (3) in 


_(2) resulteert in een differentiaal-verge- 


lijking voor m; = m(z;). In figuur 2 is een 
oplossing geschetst voor g — 0. Men kan 
zeggen dat dit profiel een uitgesmeerde 
reflectie is van de toestandsvergelijking. 
Opmerkelijk is dat een profiel overblijft 
voor g — 0, terwijl g = O alleen een homo- 
gene magnetisatie kan geven. Blijkbaar is 
de enige functie van het magneetveld 
H;= 27; om de symmetrie te breken en is 
de grootte van g niet van belang. 


De dikte van Kk van de grenslaag is in het 
algemeen van de orde van de roosteraf- 
stand a maar wanneer 7' de kritieke 7 
benadert divergeert Kk als: 

K >= (TT) * met v=1/2. (4) 
Bezwaren tegen de klassieke theorie 

De klassieke theorie heeft verschillende 
zwakke punten. De MF-benadering is on- 
bevredigend vooral in de buurt van 7. 

De kritieke eigenschappen worden door 
de MF-benadering kwalitatief verkeerd 
gegeven. De schaaltheorie van Kadanoff® 
en Widom* geeft een beschrijving van het 
kritiek gedrag waarin de singulariteiten 
van de thermodynamische functies ge- 
geven worden door samenhangende expo- 
nenten die in het algemeen willekeurige 
getallen zijn en geen eenvoudige gebroken 
getallen zoals de MF- en andere conven- 
tionele benaderingen die geven. Daarom 
moet de exponent v = 1/2 in (4) gewan- 
trouwd worden. Experimenteel!? vindt 
men inderdaad dat de grenslaagdikte ster- 
ker divergeert met een exponent v = 0,67. 
Het intrigerende van deze kritieke expo- 
nenten is dat ze universeel zijn ondanks 

de willekeurig aandoende waarde. 


Fisk en Widom® hebben dit bezwaar tegen 
de klassieke theorie ondervangen. Hun 
observatie is dat het profiel tot stand komt 
door een balancering van het effect van de 
tweede term in (3) en van de eerste term 
in (3) die op zich zelf tot de lokale, dat 
wil zeggen door H; gedicteerde waarde 
van Mm; zou leiden. Dus de tweede afgeleide 
van het profiel wordt door de toestands- 
vergelijking bepaald. Als we nu de klas- 
sieke toestandsvergelijking door de mo- 
derne vervangen dan krijgen we een pro- 
fiel dat in overeenstemming is met de 
schaalwetten. De waarde van v zoals die 
uit de schaalwetten volgt, is inderdaad 
p= 0,67 als we daarvoor de experimentele 
gegevens van de homogene fase nemen. De 
aansluiting van de Fisk-Widom-theorie 
met het experiment is bevredigend voor 
zover als de details van de grenslaag meet- 
baar zijn. 


Opmerkelijk is dat de fase-scheidende 
kracht geen rol speelt bij de tot stand- 
koming van de vorm van het profiel. We 
zullen later zien dat dit bedenkelijk is. 
Verder kleeft er aan deze procedure een 


irriterend theoretisch bezwaar. Het pro- 
fiel ziet vooral het deel van de toestands- 
vergelijking dat op het coëxistentiegebied 
slaat (de slinger in de kromme). Dit stuk 
van de kromme is niet meetbaar en theore- 
tisch kan men bewijzen dat m (H) een con- 
tinu stijgende functie is met eventueel een 
sprong. De toestandsvergelijking die Fisk 
en Widom hanteren is dus een ongeoor- 
loofde extrapolatie. Aan zo’n extrapolatie, 
die niet op strikt evenwicht kan slaan, is 
wel behoefte onder andere voor het ver- 
schijnsel van spinodale decompositie. Om- 
dat het grenslaagprofiel meetbaar is (met 
behulp van lichtverstrooiing) zou men 
eerder de zaak moeten omdraaien en de 
extrapolatie van de toestandsvergelijking 
op het grenslaagprofiel baseren. 


Renormalisatietheorie 

Het ligt voor de hand om voor een bevre- 
digende theorie van de grenslaag naar de 
renormalisatietheorie te kijken. Deze door 
Wilson opgezette theorie geeft een moge- 
lijkheid de kritieke exponenten te bereke- 
nen en de universaliteit er van vast te stel- 
len. Universaliteit veroorlooft ons ook 

om bij zulke vereenvoudigde modellen 
als het Ising-model te blijven, omdat de 
kwalitatieve aspecten in het model en 

de werkelijkheid dezelfde zullen zijn. 


Het principe van de renormalisatietheorie 
berust op het feit dat er een verband ge- 
legd kan worden tussen twee systemen 
met een verschillende maar vergelijkbare 
hamiltoniaan # en #. De hamiltoniaan 
is de uitdrukking die in de exponent van 
de Boltzmann-factor verschijnt om het 
relatieve gewicht van een configuratie te 
geven. Impliciet bevat #” dus ook het aan- 
tal deeltjes N (vrijheidsgraden, spins) en de 
temperatuur 7’, naast de interactie para- 
meters zoals het veld H, de naaste-buren- 
koppeling J enz. 4” is een hamiltoniaan 
met minder vrijheidsgraden dan ”/. Een 
verband tussen # en#{” ontstaat door het 
middelen over ongevaarlijke vrijheids- 
graden. In de veldentheoretische aanpak 
zijn dat de snel variërende componenten 
van het orde-parameterveld. Men kan ook 
een methode volgen, die speciaal in het 
Ising-model werkt, waarbij blokken van 
spins tot een nieuwe effectieve spin ver- 
enigd worden en #” is dan de effectieve 
hamiltoniaan voor deze nieuwe spins. 


Herhaald toepassen van de renormalisatie 
levert een rij #4, /, … de renormali- 
satie trajectorie in de #/-ruimte. Langs een 
trajectorie ligt een rij systemen met ver- 
gelijkbare eigenschappen. Correspondeert 
de begin #/ met een situatie van coëxi- 
sterende fasen dan ligt de hele trajectorie 
in het coëxistentie gebied. Renormalisatie 
is op te vatten als een stroom langs de re- 
normalisatie trajectoriën. Een schets van 
de projectie van deze stroom in het 7-H- 
vlak voor het Ising-model is gegeven in 
figuur 3. Trajectoriën mogen faseschei- 
dingslijnen niet doorkruisen en monden 
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Fig. 3. Renormalisatie stroompatroon 
geprojecteerd op de T-H-ruimte van het 
Ising-model. 


uit in vaste punten. Het Ising-model heeft 
onder andere vaste punten voor H = O en 
T=0,T=T, en 7=oo. Punten met 

T > Te stromen naar T= coen met TS 7e 
naar 7'= 0. De eindpunten geven het 
systeem onder extreme omstandigheden. 
Een situatie 7 > T, gaat naar de volledige 
wanordelijke situatie bij 7'= oo en een 
systeem met 7'< 7, eindigt bij volledige 
orde op 7'= 0. 


Een renormalisatietransformatie is slechts 
van waarde wanneer de eigenschappen van 
de systemen langs de trajectorie op een- 
voudige wijze met elkaar in verband ge- 
bracht kunnen worden. In eerste instantie 
betreft dat de thermodynamische eigen- 
schappen, maar ook de correlatiefuncties 
en het grenslaagprofiel kan langs de tra- 
jectorie bestudeerd worden. De eigen- 
schappen van de coëxisterende fasen kun- 
nen gevonden worden door te analyseren 
hoe ze langs de trajectorie veranderen en 
hoe ze bij het eindpunt 7'= Ò zijn. De kri- 
tieke singulariteiten volgen uit de manier 
waarop de stroom stagneert bij 7= 7. 
Met name de veldentheoretische aanpak 
heeft via de e-expansie (e = 4 minus de 
dimensie van het systeem) met indruk- 
wekkende precisie de kritieke exponen- 
ten berekend.” Bijvoorbeeld vindt men 
voor v de waarde 0,6300 + 0,0008 het- 
geen binnen de experimentele meetnauw- 
keurigheid valt. 


Toepassing van de renormalisatietheorie 
op het grenslaagprofiel betekent dat men 
het profiel moet berekenen uit het pro- 
fiel op 7'= 0 en uit de manier waarop het 
transformeert langs de trajectorie. Het ligt 
voor de hand om bij 7'= 0 een scherp 
grensvlak te veronderstellen. Zo kan men 
uit de renormalisatievergelijking het pro- 
fiel vastleggen bij eindige 7 zonder te moe- 
ten leunen op de extrapolatie van de toe- 
standsvergelijking in het coëxistentiege- 
bied. In feite wordt een nauwer verband 
gelegd met de correlatiefuncties, waardoor 
de exponent v in (4) dezelfde wordt als 
waarmee de correlatielengte divergeert bij 
T,. De resultaten van de renormalisatie- 
theorie sluiten aan bij de uitkomsten van 
Fisk en Widom. j 
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Oppervlakte golven 

Tot nu toe hebben we ons weinig bekom- 
merd om de vraag welke configuratie sig- 
nificant bijdragen tot de grenslaagdikte. 
Gaan we uit van een ideaal grensvlak op 
z = 0 dan kunnen we twee typen onder- 
scheiden: golven in het grensvlak en het 
verschijnen van eilandjes van —spins voor 
z > 0 en + spins voorz <0. 


Buff, Lovett en Stillinger? hebben naar de 
invloed van golven gekeken. Zij gaan uit 
van een grensvlak dat een oppervlakte- 
spanning oo heeft. Een golf leidt tot ver- 
groting van het grensvlak en zal een extra 
energie bevatten gelijk aan oo maal de 
oppervlaktevergroting. Daarnaast verzet 
de zwaartekracht zich tegen het ontstaan 
van golf. Als A (x,y) de hoogte van de 
golf ten opzichte van het ideale grensvlak 
is zal de energie E van een golf met niet al 
te kleine golflengte gegeven worden door: 


E(h)= f dxdy {oo IVA? +gAph? }/2. 
S 
(5) 


De eerste term vertegenwoordigt de opper- 
vlaktevergroting en de tweede term is de 
potentiële energie van een kolom ter 
hoogte h en een verschil massadichtheid 
Ap. Een golf komt voor met een kans 
gegeven door exp—E/kT. De golven leiden 
tot een modificatie van og en ook tot een 
grenslaagdikte K. Voor g — 0 volgt in 
dominante orde het verontrustende resul- 
taat: 


(6) 


Dus als g — 0 divergeert K. Voor een grens- 
lijn tussen twee coëxisterende fasen in 
twee dimensies zou het nog erger Zijn om- 
dat dan kK > g” 4. Dus de zwaartekracht 
moet eindig blijven wil de grenslaag niet 
onbeperkt groeien. De lange golven zijn 
hieraan debet en voor deze is (5) correct. 
Dit werpt een schaduw over alle voor- 
gaande beschouwingen waarin de limiet 
g> 0 ongevaarlijk leek zoals in de MF- 
benadering en de daarop geënte theorieën. 
Omdat de logaritme maar weinig met 2 
verandert is het effect niet dramatisch. 
Voor water bij kamertemperatuur be- 
draagt de logaritme in (6) een factor 6, 
zodat de schaal van Kk nog steeds door 
(kKT/oo)? bepaald wordt die van de orde 
van 10-8 cm is. 


Kk =- {(KT/o0) log (05/gAp kT}. 


De vraag blijft wat de betekenis is van deze 
continuüm beschouwing voor een micro- 
scopisch model. 


Het SOS-model 

Een microscopisch model waarin de im- 
portantie van oppervlaktegolven bestu- 
deerd kan worden is het SOS-model 
(SOS = Solid On Solid), een model dat 
uit de kristalgroei stamt. We kunnen het 
beschouwen als een ontaard geval van het 
Ising-model waarin de koppelingscon- 





Fig. 4. SOS-configuratie in het Ising- 
model. De getrokken lijn is de hoogte h 
tussen de twee fasen. 


stante J in één richting (de z-richting) on- 
eindig wordt. Deze limiet beperkt het aan- 
tal configuraties die tot de grenslaag bij- 
dragen drastisch. Eilandjes van verkeerde 
spin zijn zonder meer uitgesloten omdat 
de creatie daarvan oneindig veel meer 
energie zou vragen. Wanneer er een grens- 
vlak in het systeem loodrecht op de z-as 
geforceerd wordt zijn alleen deviaties van 
het ideale grensvlak mogelijk waarbij geen 
extra Z-binding verbroken zijn. Dit bete- 
kent een effectieve beperking tot opper- 
vlaktegolven die ook niet mogen over- 
hangen (zie figuur 4). 


Het SOS-model is een substantiële reductie 
van het Ising-model, maar niettemin nog 


| steeds niet oplosbaar in drie-dimensies. De 


spinconfiguraties die meedoen kunnen be- 
schreven worden door een discrete func- 
tie h die aangeeft waar het scheidingsvlak 
tussen + en — spins ligt. De energie van de 
toestand is evenredig met de extra grootte 
van het scheidingsvlak (dat wil zeggen het 
aantal verkeerde verbindingen in de x- en 


_Y-richting) en met een magnetische veld- 


term: 
Eh)=J E Ih-hjl-igEh;, M 
<ij> i Ol 


ien j doorlopen nu een twee (d—1) dimen- 
sionale verzameling en /, alle waarden O0, 
+1,+2, … (7) is een discrete versie van 
(5)en dit stelt ons in staat om de conclusie 
van Buff, Lovett en Stillinger op micro- 
scopisch niveau na te gaan. De twee- 
dimensionale versie van (7) is exact oplos- 
baar’ en leidt inderdaad voor alle 7'tot 
een divergentie van Kk als g=% zoals de 
oppervlaktegolven-theorie voorspelt. 


Het 3-dimensionale SOS-model (7) heeft 
een nieuw verschijnsel: de verruwingstem- 
peratuur 7, heteerst ontdekt in de kristal- 
groeil®. Voor T >T, gedraagt het SOS- 
model zich als de oppervlaktegolven- 
theorie voorspelt, voor 7 <T, blijft het 
oppervlak glad en eindig van dikte als 
g> 0. De renormalisatietheorie van Kos- 
terlitz en Thouless!! kan dit verschijnsel 
verklaren en de aard van deze faseover- 
gang in de grenslaag vaststellen. Het blijkt 
dat de overgang in de oppervlakte vrije 
energie uiterst zwak is: van eo-hoge orde, 
dus alle afgeleiden continu. 


Hiermede zijn we aangeland op een terrein 
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waar veel interessante verbanden, onder 
andere met het ontstaan van vortices in 
twee-dimensionale systemen, gelegd kun- 
nen worden en wat nog volop in ontwik- 
keling is. De hoop is dat langs deze weg 
een volledig bevredigende theorie van de 
grenslaag zal volgen. 
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Katastrofe-theorie; kwalitatieve theorie 
van niet-lineaire dynamische systemen 


J. Hijmans 


De katastrofe-theorie houdt zich niet specifiek bezig met onheil, maar met onver- 
wacht of abrupt verlopende, kwalitatieve overgangsverschijnselen. De theorie dankt 
zijn naam aan de Franse wiskundige René Thom, die met zijn in 1972 gepubliceer- 
de boek ““Stabilité structurelle et morphogenêse”’! enige opwinding teweegbracht; 
aanvankelijk alleen onder zijn vakbroeders in de algebraïsche topologie, maar later 
ook in een veel bredere kring van biologen, sociologen, psychologen, fysici en che- 
mici. Aan deze snelle opmars van de katastrofe-theorie waren de populariserende 
lezingen van de Engelse wiskundige Christopher Zeeman? en andere propagandis- 
ten zeker niet vreemd. Zowel de opwinding onder bewonderaars als de kritiek van 
de zijde van de ongelovigen® vinden hun verklaring ongetwijfeld in het wat onge- 
bruikelijke karakter van de theorie. Katastrofe-theorie is namelijk geen kwantita- 
tieve, van welomschreven grondprincipes uitgaande, theorie, die via exact gefor- 
muleerde wetten uitmondt in kwantitatieve, experimenteel toetsbare voorspellin- 
gen, maar een op universaliteit gerichte, algemene methodologie voor het klassifi- 
ceren en herkenbaar maken van kwalitatieve overgangsverschijnselen en voor het 
identificeren van de relevante variabelen, die voor het verloop van die verschijnse- 


len bepalend zijn. 


Toestands- en controlevariabelen 

In de katastrofe-theorie* maakt men on- 
derscheid tussen twee soorten variabelen: 
de toestandsvariabelen, vergelijkbaar met 
de extensieve grootheden (““densities’’) in 
de thermodynamica, en meestal aange- 
duid met x, ‚Xx, …, X,,, Zijn de variabelen 
waarin het verloop van het bestudeerde 
proces wordt uitgedrukt. De waarden van 
deze grootheden kunnen sprongen of 
abrupte veranderingen ondergaan, die als 
‘katastrofes’ worden aangeduid. De con- 
trolevariabelen, vergelijkbaar met de in- 
tensieve- of reservoirparameters (“fields’”) 
in de thermodynamica, hier aangeduid 
met u, ‚U, … U, Zijn de variabelen, 
waarmee het systeem gestuurd kan wor- 
den en die bepalen of een overgang al of 
niet plaats kan vinden. 

Als voorbeeld kan men denken aan een 
systeem van elkaar beïnvloedende che- 
mische reacties. Hierbij zijn dan de con- 
centraties x,,X5, … X,, van de n stoffen 
die bij die reacties geproduceerd of ver- 
bruikt worden de toestandsvariabelen, 
terwijl de temperaturen, drukken en che- 
mische potentialen, die voor het verloop 
van de reacties bepalend zijn, de rol van 
controlevariabelen spelen. Een dergelijk 
proces wordt beschreven door een sys- 
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teem van n gekoppelde 1° orde differen- 
tiaalvergelijkingen: 


de Pdf A, rs Br dn 
Es rs Ui) 

de dt As rs Wor nas (1) 
rs Uns var U | 


de, Jdf == A, (45% oes HS 
Bin Mzan 


of afgekort: 


dx/dt =A(x;u), (2) 
waar A (Xx, u) een afkorting is voor n, in 
het algemeen niet-lineaire functies van de 
concentraties x, ‚Xx, …, X,, die parame- 
trisch van de mm controlevariabelen u, , 
Us, … U, afhangen. 


Het faseportret van een dynamisch 
systeem? 

De katastrofe-theorie is niet primair ge- 
interesseerd in het oplossen van vergelij- 
kingenstelsels zoals (1), maar houdt zich 
bezig met de klassificatie van de mogelij- 
ke singulariteiten van een dynamisch sys- 
teem, dat wil zeggen zij bestudeert de 
ligging en de aard van die punten in de 
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controleruimte, waar de overgangen van 
het ene régime naar het andere plaats vin- 
den. 

Zo kan men bijvoorbeeld de bekende dif- 
ferentiaalvergelijking van de mathemati- 
sche slinger: 


d20/dt? =— (g/O sin 0, (3) 


waar 0 de uitwijkingshoek, Q de lengte van 
de slinger en g de zwaartekrachtsversnel- 
ling is, opvatten als een eerste-orde dyna- 
misch systeem in twee variabelen, 0 en 

w, namelijk: 


d8/dt = w;dw/dt=—(g/Q) sind. (4a, b) 


Van dit systeem kan men — gebruikma- 
kend van het behoud van de energie £ van 
de slinger — de baanvergelijking expliciet 
bepalen, maar voor de katastrofe-theorie 
is het slechts van belang, dat de toestands- 
ruimte, dus het (0, w)-vlak, uiteenvalt in 
kwalitatief verschillende gebieden, die 
Nederlands Tijdschrift voor Natuurkunde A47 (4) 1981 


w 


ENE 
NZLSE 


SPR) } 


E > mgql 


Ön 


IKO 


Fig. 1. Het faseportret van de mathematische slinger. 


gescheiden worden door zogenaamde 
separatrices (zie figuur 1): 


w=t Vgl cos z0,E= met. (Sa, b) 


De banen binnen de separatrices (E < 
mg) stellen periodieke oplossingen, dus 
slingeringen, voor, terwijl de buiten de se- 
paratrices gelegen banen (E > mg) met 
zuivere rotaties ‘over de top’ van de slin- 
ger corresponderen. De punten w = 0, 
0O=?2nn,(E=0Oenw=0,0=(2n+1) Tm, 
(E = mgQ), zijn elliptische (stabiele) res- 
pectievelijk hyperbolische (metastabie- 
le) vaste punten van het systeem. 

De slinger (4) is een voorbeeld van een 
Hamilton-systeem, dat wil zeggen een sys- 
teem waarvoor de bewegingsvergelijkin- 
gen kunnen worden uitgedrukt in de Ha- 
milton-vorm: 


| de òH(p,0) dp _oH(p,6) 
dt òp er E 90 
(6a, b) 
waar: 
H(p,0) =p? /2mQ? — mg& cos0, (6c) 


de Hamilton-functie is, uitgedrukt in 0, 
en de draai-impuls p = mQ? vo 

De moeilijkheid met Hamilton-systemen 
is, dat zij niet structureel stabiel zijn, dat 
wil zeggen dat hun “faseportret’ (de ver- 
zameling mogelijke banen in de toe- 
standsruimte) kwalitatief verandert door 
toevoeging van een infinitesimale dissipa- 
tieterm. Zo worden bijvoorbeeld in het 
geval van een zwak gedempte slinger, de 
afzonderlijke gesloten banen van figuur 1 
verenigd tot spiralen, die op een gecom- 
pliceerde manier door elkaar lopen (zie 
figuur 2). 


Daar men in praktische situaties eigenlijk 
altijd te maken heeft met aan kleine ver- 
storingen en fluctuaties onderhevige sys- 
temen, denkt men zich in de katastrofe- 
theorie een structureel instabiel systeem 
altijd opgenomen in een structureel sta- 
biele familie van naburige, licht verstoor- 
de systemen, geparametriseerd door be- 
paalde controleparameters, zoals de 
dempingsconstante in het voorbeeld van 
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Fig. 2. Faseportret van de gedempte 
slinger. 


de slinger. Deze structureel stabiele fami- 
lie wordt de ontplooiing (unfolding) van 
het systeem genoemd. 

Het faseportret van een structureel stabiel 
systeem (of familie van systemen) kan 
worden ingedeeld in stroomgebieden van 
bepaalde attractoren. Een attractor is een 
punt in de toestandsruimte ( een vast 
punt) of een verzameling punten (bijvoor- 
beeld een limiet cycel of een invariante 
torus), waartoe op den duur alle banen 
uit het stroomgebied willekeurig dicht na- 
deren. Het in figuur 2 weergegeven fase- 
portret van de gedempte slinger heeft al- 
leen punt-attractoren, namelijk de pun- 
ten w = 0, 0 = 2nm. Andere systemen 
(zoals bijvoorbeeld de Van der Pol-oscilla- 
tor) hebben echter ook limiet cycli, of 
meer ingewikkelde (“strange”) attracto- 
ren. 


Gradiëntsystemen 

Het meest bestudeerde (en, tot dusver het 
enige min of meer uitputtend in kaart ge- 
brachte) gebied van de katastrofe-theo- 
rie, is de dynamica van gradiëntsyste- 
men, die bepaald wordt door de mogelijke 
singulariteitstypen van potentiaalfunc- 
ties. 

Een gradiëntsysteem is een dynamisch 
systeem van de algemene vorm (2), waar- 


bij echter de evolutievector A (x;u) kan 


worden geschreven als de negatieve gra- 
diënt van een potentiaalfunctie: 


A (xu)=—V,V(x, 4). (7) 


Een gradiëntsysteem heeft uitsluitend 
punt-attractoren, namelijk de (attrac- 
tieve) stationaire punten xo van het sys- 
teem waar geldt: 


6) 


In de omgeving van die punten xo, waar 
de Hessiaan van de potentiaal ongelijk nul 
is, dus waar: 


Nee V(x;u)=0. 


det | V;; #0, 


Vj= 5 (9° V‚‚lòx,òx,), (9a, b) 
wordt het kwalitatieve gedrag volledig be- 
paald door het gelineariseerde systeem: 


dx; — Xi) 
dt 


n 
EE Vw — 


j=t X,0) — Xjo) + vs 


(10) 


dat door middel van een gelijkvormig- 
heidstransformatie x —x, = Sx' altijd op 
de vorm: 


dx;/dt= — Nx; + O(x!2) 


N#0(i=1,2,……,n), (11) 
kan worden gebracht. Een punt x, ‚waar 
alle eigenwaarden ), positief zijn is een 
attractor, een punt met uitsluitend nega- 
tieve )\, een repellor, en een punt met zo- 
wel positieve als negatieve )\, heet een ge- 
generaliseerd zadelpunt. Deze typen sta- 


_tionaire punten waarvoor alle ), # 0 


zijn, heten in de mathematische litera- 
tuur Morse-singulariteiten. 

In die punten x,, waar één of meer 
eigenwaarden van het gelineariseerde sys- 
teem, zeg À\, =À, =..= Ag = nul wor- 
den, terwijl de overige, Ag , > Àg 42» 
\, #0 zijn, wordt het lokale gedrag 
rondom de singulariteit essentieel bepaald 
door de hogere-ordetermen uit vergelij- 
king (10). 

Volgens de zogenaamde splitsingsstelling 
van Morse kan men in die punten altijd 
de n-toestandsvariabelen splitsen in & 
relevante variabelen en n — & irrelevante 
variabelen en het gedrag in de buurt van 
de singulariteit beschrijven als een Morse- 
singulariteit ten aanzien van de irrelevan- 
te variabelen, gesuperponeerd op een zo- 
genaamde residuale singulariteit in de re- 
levante variabelen alleen. De klassificatie 
van deze residuale singulariteitstypen 
vormt het eigenlijke onderwerp van 
Thoms theorie der elementaire katastro- 
fes. 


Het klassificatie- en het 
ontplootingsprobleem 

De residuale singulariteitstypen van een 
potentiaalfunctie worden bepaald door 
twee karakteristieken, namelijk de co- 


rang en de codimensie. De corang van 
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een singulariteit van een potentiaalfunc- 
tie, V(x), is de dimensie van de nul- 
ruimte van zijn Hessiaan (99), dat wil 
zeggen het aantal relevante toestands- 
variabelen, dat de residuale singulariteit 
karakteriseert. De codimensie kan men 
opvatten als het minimale aantal sym- 
metrie brekende verstoringen, die nodig 
zijn om het systeem terug te voeren tot 
in ‘algemene positie’ in de toestands- 
ruimte, dat wil zeggen om het systeem 
te veralgemenen tot een systeem met 
uitsluitend Morse-singulariteiten. 
De twee hoofdproblemen van de kata- 
strofe-theorie zijn enerzijds het klassifi- 
catieprobleem, dat wil zeggen het vinden 
van de elementaire katastrofe waarmee 
een bepaalde singulariteit “diffeomorf” 
is, en anderzijds het ontplootingspro- 
bleem, dit is het inbouwen van de sin- 
gulariteit in een familie van potentiaal- 
functies, die de betreffende singuliere 
potentiaal als bijzonder geval bevat, maar 
die beschouwd als familie in ‘algemene 
positie’ is. Twee functies zijn diffeomorf 
als er een (in het algemeen niet-lineaire) 
één-eenduidige, continue, differentieer- 
bare, dus omkeerbare transformatie be- 
staat, die de ene functie in de andere 
overvoert (in de omgeving van een sin- 
gulier punt). 
De universele ontplooiing van een singu- 
lariteit is een zogenaamde minimale ver- 
sale familie van (potentiaal) functies, 
V(x,u), die behalve van de & relevante 
toestandsvariabelen, x, ‚Xx, …, Xg, af- 
hangt van evenveel controlevariabelen 
Ui, Ur, Up als de codimensie (m) be- 
draagt, en die voldoet aan de eis, dat 
iedere functie in een infinitesimale om- 
geving van het singuliere punt V(x,0) 
diffeomorf is met een functie uit de fa- 
milie V(x,u). Zo heeft bijvoorbeeld de 
functie: 
V(x)=x*, (12) 
in de oorsprong een residuale singulari- 
teit van corang 1. Een infinitesimale om- 
geving van deze functie kan altijd worden 
weergegeven door de familie van alle 4°- 
graads polynomen: 


V(x,a, D, C, d) en 


=X* +AXS +bx? +ex+d, (12a) 
met infinitesimale coëfficiënten a, b,c en 
d. Deze familie is weliswaar versaal (bij- 
dragen van 5°“- en hogere graad kunnen 
willekeurig klein gemaakt worden ten op- 
zichte van x*, door dichter tot de oor- 
sprong te naderen), maar niet minimaal 
versaal, want door verschuiving van de 
nulpunten in x en V is het altijd moge- 
lijk (12a) te transformeren op de vorm: 


V(x,uj,ug)=X* +UjX* +Uuyx, (12b) 
die diffeomorf met (12a) is maar slechts 


van twee controle-variabelen u, en uz af- 
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hangt. Dat deze familie inderdaad een 
universele ontplooiing van (12) is kan 
men zien aan het feit, dat het vastleggen 
van de posities van de drie extrema die 
een 4°-graadsfunctie kan hebben, de 
specificatie van minimaal twee additio- 
nele parameters vereist. 


De singulariteitstypen van corang 1 en 2 


De elementaire katastrofes van corang Ì, 


de zogenaamde cuspoiden, worden be- 
paald door potentiaalfuncties in één 
variabele x, van de algemene vorm: 
Vjs K=) (Ari) 

Deze singulariteitstypen zijn door Thom 
voorzien van aan meetkundige voorstel- 
lingen ontleende namen, zoals de “fold” 
(k =3), de “cusp” (k = 4), de “swallows- 
tail” (k = 5) en de “butterfly” (k = 6). 
De residuele singulariteitstypen van corang 
2 zijn potentiaalfuncties V(x,,X2 ), waar- 


voor: 
ov |) _ 
0 


(ae), ={ 
dx, * 0 E ÒX 1 ÒXs 


(5e) 

mn! 2 pam) 0, 
0X2“ /9 

dus functies die gedomineerd worden door 
de 3°- of hogere-graadstermen in hun Tay- 
lor ontwikkeling om het singuliere punt. 
Wij bekijken allereerst het 3°-graads- 


gedeelte van de potentiaal, dat algemeen 
geschreven kan worden in de vorm: 











(14) 


Vz (X1,%2)=alXi HH X2) Xi + 


+2 X2) (Xi +H3X2), (15) 
en waarbij wij de volgende gevallen kun- 
nen onderscheiden: — 


1. De potentiaal heeft 3 reële verschil- 
lende nullijnen: u; Hus #3 = reëel. 

In dit geval kunnen wij de potentiaal door 
middel van een geschikt gekozen (lineaire) 
transformatie brengen op de standaard- 
vorm: 


V(Xi,X2)= (Xi +X2) (A1 — X2)%2= 


(Da) (16a) 


die door Thom wordt aangeduid als de 
“elliptic umbilic”. 

2. De potentiaal heeft 1 reële en 2 toe- 
gevoegd complexe nullijnen: u; = U; 

Ma =H*‚H3 = reëel. 

In dit geval is de potentiaal diffeomorf 
met de standaardpotentiaal van de hyper- 
bolice umbilic: 


2 3 
=Xi X2 TX 5 


V(xiXa)= Xi X2 +X2°. (Da) (16b) 


3. Wanneer de potentiaal twee samenval- 


* De aanduidingen met hoofdletters zijn door | 
VL Arnol’dé ingevoerde, aan de groepentheorie 
ontleende, symbolen. 


lende nullijnen heeft: u; =H» en u3 reëel, 
is hij diffeomorf met de potentiaal: 


(17) 


maar deze functie is niet eindig bepaald 
(dat wil zeggen heeft codimensie onein- 
dig). Daarom is het in dit geval essentiëel, 
dat wij ook de hogere-graadstermen ex- 
pliciet kennen. De eenvoudigste manier 
om de codimensie van de functie (17) 
eindig te maken is door toevoeging van de 


V(xiX)= Xi 3, 


_4°-graadsterm + x5*. Wij krijgen dan de 


standaard vorm van de “parabolic um- 
bilic”: 

Vg sai Hate (17a) | 
De enige andere manier om de singulari- 
teit (17) te ‘stabiliseren’ blijkt toevoeging 
van een hogere-graadsterm + x3 £ te zijn. 
Op die manier vindt men de complete 
reeks van ““conic umbilics””: 
VX, Xa) Xi Xa tak (Dis). (17b) 
4. Het laatste geval dat wij zullen be- 
schouwen is het geval, dat het 3°-graads- 
gedeelte van de potentiaal drie samenval- 
lende nullijnen heeft: U, =H2 = H3. 


In dit geval is de potentiaal diffeomorf 
met de functie: 


(18) 


die natuurlijk evenmin eindig bepaald is 
als (17). Het ‘stabiliseren’ van de poten- 
tiaal (18) tot een eindige bepaalde functie 
blijkt nu slechts mogelijk door toevoegen 
van een 4°-graadsterm + x3* of xj Xx2° 

of een 5°-graadsterm x,°. Dit geeft aan- 
leiding tot nog drie singulariteitstypen, 
namelijk: 


Vi) 


V(xi,xa)= Xi 23°, (E6) (18a) 
V(xj,x2)= Xi +X1X2, (Ey) (18b) 
en: 

rs ki A: (Eg) (18c) 


Weliswaar is hiermee de lijst van singulari- 
teitstypen nog niet uitgeput, maar alle 
overige typen, zoals bijvoorbeeld de func- 
tie: 


Er Fk: 
Tio) (19) 


hangen parametrisch af van één of meer 
continu variërende grootheden (zoals 
hier a), die mede bepalend zijn voor hun 
type. Arnold noemt deze singulariteits- 
typen, die van één of meer type bepalende 
parameters afhangen de één- of in het al- 
gemeen n-modale singulariteiten. De door 
Thom onderzochte elementaire katastro- 
fen, corresponderen met Arnol’ds 0-mo- 

dale singulariteiten. | 
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V(Xi,x2)=X1 


De universele ontplooiing van de 
elementaire katastrofen | 
Zoals wij al aan het voorbeeld (12b) van 
de cusp hebben gezien, zijn de cuspoïden 
(13) algemeen van codimensie k — 2,en 
kunnen hun universele ontplooiingen wor- 
den weergegeven door de (k — 2)-parame- 
terige familie 

k—2 
V(x‚u)=xk + > u;xi. (20) 

j= 
Dit is een consequentie van het feit, dat 
een k®-graadspolynoom in één variabele 
altijd tot de vorm (20) kan worden gere- 
duceerd door verschuiving van het nulpunt 


van x en aftrekken van de constante term | 


in V. 


Om de universele ontplooiing te vinden 
van de elementaire katastrofen van 
corang 2, zoals bijvoorbeeld de conic 
umbilics (17), sporen wij die functies op 
die diffeomorf zijn met de (singuliere) 
potentiaalfunctie zelf en die dus dienen te 
ontbreken in de universele ontplooiing 
van die functie. Beschouw bijvoorbeeld de 
functie: 
Va,y)= xt +y°, (21) 
die diffeomorf is met de hyperbolic 
umbilic (16b). Toepassen op (21) van een 
diffeomorfisme: 


XX =Xtdlx,y) H(O0,0)=0 (22a) 
met 
yay =y+VX,y) _ V(O,0)=0,(22b) 


en ontwikkelen in een Taylor-reeks, laat 
zien dat alle functies: 


Ny EE òV 
Vp) Vr EVP) + es 


v dv 23 
deni. 

En 23) 
diffeomorf zijn, met de oorspronkelijke 
potentiaal, en dus dat in de ontplooiing 
van (21) alle termen evenredig met: 

WV ge 


= IX 


Òx (20) 


of aas =3p*, 
py 

of met willekeurige combinaties daarvan 

dienen te ontbreken. Hieruit kunnen wij 

concluderen, dat in de ontplooiing van 

(21) alleen diè monomen voor kunnen 

komen, die in het schema: 


(25) 





buiten de ‘schaduwen’ van x° en y? val- 
len, met andere woorden dat de ont- 
plooiing van (21) wordt gegeven door de 
familie 


V(x‚v,u,v,w)=X? +3 +UX +UY +WXY. 


(26) 
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Het vinden van de universele ontplooiing 
van andere katastrofen van corang 1 , zoals 
de elliptic umbilic die traditioneel wordt 
voorgesteld door de met (16a) diffeo- 
morf functie: 


(27) 


is minder triviaal, omdat een der compo- 
nenten van de gradiëntvector: 


ov òV 


1 
Va) = 3 x° —xy?, 


=p. 1 =-2xy, 28 
a xy (28) 
geen monoom is. De combinaties: 
à òV , òV „ 
zy 
SE De Dap 
òV òV 
PN Wa 4 (29) 
Ox òy 


zijn dit echter wel, en dus kunnen in de 
ontplooiing van (27) alleen die monomen 
voorkomen, die niet behoren tot de scha- 
duwen van x°, y® of xy in het schema 
(25). Dit zijn x,y, x? en y?, maar van de 
laatste twee mag het verschil niet voor- 
komen, en dus alleen de som, x? + y?. 
Hiermee vindt men voor de ontplooiing 
van de elliptic umbilic: 


Ì 
V(x,y,u, vw) = Ei Kx + 


+UX HUP +Ww (Xx? +y?). (30) 
Analoge redeneringen leiden tot de ont- 
plooiingen van andere singulariteiten van 
corang 2. 


Meetkundige voorstelling 

Als men eenmaal de universele ont- 
plooiing van een elementaire katastrofe 
V(x,u) heeft gevonden, kent men in 
principe: 


i. het evenwichtsoppervlak, di. de 
verzameling van stationaire punten in 
de (Xx, 4u)-ruimte waar geldt 

VV, u) =0, 31) 
ii. de bifurcatieverzameling, bestaande 
uit die punten van de controle (u-)-ruimte 
waar behalve (31) geldt dat: 

det |V, (x,u) |= (32) 
Een volledige meetkundige voorstelling 
van de bifurcatieverzameling is slechts 
mogelijk voor de katastrofes van codi- 
mensie 3 of lager, dus voor de fold, de 
cusp,‚ de swallowstail, de elliptic- en de 
hyperbolic umbilic. 

Voor de “fold”: 


V(«x,u)=Xx° +ux 
is het evenwichtsoppervlak de parabool: 


òV/òx=3x? +u=0, (34) 


G3) 


a) 


b) 
diein jn 


Fig. 3. De “fold” katastrofe. a) het even- 


wichtsoppervlak 3x? +u = 0, b)de bifur- 
catieverzameling, u = 0. 


en de bifurcatieverzameling bestaat uit 
één enkel punt namelijk u = 0 (zie figuur 
Jk. 

Voor de ““cusp”: 


V(x,ur,uz)=X* +UX? +UX, (35) 
is het evenwichtsoppervlak: 
AV/dx=4xS +2ux +u=0, (36a) 


een 3°-graadsoppervlak in de (x, u, v)- 
ruimte (zie figuur 4a), terwijl de bifur- 
catieverzameling wordt gegeven door de 
lijn: 

8u? +27u° =0, (36b) 


in het (u, v)-vlak (figuur 4b), waar de dis- 
criminant van de vergelijking (36a) nul is. 


Voor de “swallowstail’’: 


V(x,u,v,w)=X° +ux* +ux? +wx, (37) 


Fig. 4. De “cusp” katastrofe, a) even- 
wichtsoppervlak, b) bifurcatieverzame- 
ling. 


a) 


ge 


b) 
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Fig. 5. Het bifurcatievlak van de 
“swallowstail”-katastrofe in de controle- 
ruimte. 


kan het evenwichtsoppervlak, dV/òx = 0, 
niet meer als geheel weergegeven worden, 
omdat het 4-dimensionaal is, maar de 
bifurcatieverzameling, die ontstaat door 
x te elimineren uit de condities à V/òx = 0 
en 9? V/dx? =O, heeft de in figuur 5 ge- 
schetste vorm. 


Ook voor de beide andere elementaire 
katastrofes van codimensie 3, de “elliptic’” 
en de “hyperbolic” umbilic, kunnen wij 
alleen de bifurcatieverzameling weergeven. 
Voor de elliptic umbilic (30) wordt het 
bifurcatievlak bepaald door uit de verge- 
lijkingen àV/òx =àV/dy = Oen (9° V/dx®) 
(92 V/ay?) — (92 V/axdy) =O, x en y te 
elimineren. Dit levert het in figuur 6 ge- 
geven oppervlak bestaande uit twee 
‘doornen’. 


Evenzo bepaalt men voor de hyperbolic 
_umbilic (28): 


V(x,‚y,u,v,w)=X? +7? +WXY + 


+UX +UYy, 


Fig. 6. Het bifurcatievlak van de “ellip 
tic umbilic”, | 


al 
pe 
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(28) | 





Fig. 7. Het bifurcatievlak van de “hyper- 
bolic umbilic”’ 


het in figuur 7 geschetste bifurcatie- 
oppervlak. 


Wij bespreken nu een drietal fysische toe- 
passingsgebieden van de theorie der ele- 
mentaire katastrofen. 


Brandvlakken in de optica” 

Een brandvlak kan ontstaan wanneer een 
evenwijdige bundel lichtstralen door een 
gekromde spiegel gereflecteerd of door 
een lens gebroken wordt. Het is de omhul- 
lende van alle gereflecteerde of gebroken 
stralen. In de limiet van de geometrische 
optica (golflengte — 0) wordt de intensi- 
teit oneindig in de punten van een brand- 
vlak (in de golfoptica wordt de intensiteit 
groot maar blijft hij eindig). 

Een bekend voorbeeld van een brandlijn 
is de scheidingslijn tussen een lichte en 
een donkere zone, die men ziet op het 


hen 


Fig. 8. “Cusp”’katastrofe op het opper- 
vlak van een theekopje (a); stralengang 


(b). 


a) 





b) 


d 
del 
Sj 
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Fig. 9. Lichtweg tot aan een punt 
Pf/x,y), via een punt O(cos, sind) op de 
rand van het kopje. 


oppervlak van een door de zon beschenen 
theekopje (zie figuur 8). 

Het ontstaan van deze brandlijn kan men 
gemakkelijk verklaren aan de hand van 
het principe van Fermat (principe van de 
kortste weg). Stel het oppervlak van de 
thee in het kopje voor door een cirkel van 
straal één, en laat een (praktisch horizon- 
tale) evenwijdige lichtbundel invallen in 
de x-richting (zie figuur 9). 


Voor een lichtstraal, uitgaande van een 
punt op afstand d links van de oorsprong, 
die het punt P met coördinaten (x,y) 
bereikt, na tegen de rand van het kopje 

te zijn gereflecteerd in een punt Q met 
coördinaten (cos6, sin), is de weglengte: 


fC«‚y,0)=d + c0s0 + [(x — cosl)? + 
+ (y — sin0}?}?. (38) 


Volgens het principe van Fermat moet 
deze lichtweg minimaal zijn als functie 
van de hoek 9, die de plaats van het punt 
Q bepaalt. Dit levert de conditie: 
òf(x,v, 0/00 =0, (39a) 
die na enige herleiding kan worden uitge- 
drukt in de vorm: 


(y — sin) cos20 = (x — cosl) sin20, 
(39b) 


waarin men gemakkelijk de reflectiewet 
(hoek van inval = hoek van terugkaatsing) 
kan herkennen. 

Een brandvlak wordt nu gevormd door 
die punten (x,y) waar behalve aan (39) 
ook voldaan is aan de conditie: 
d°f/o0* =0, (40) 
die betekent, dat (39) een dubbele wortel 
heeft. 

De twee voorwaarden (39) en (40) bepa- 
len tezamen een kromme in het (x, y)- 
vlak, die nephroïde heet. 

In dit geval speelt de hoek 9 de rol van 
toestandsvariabele, terwijl de coördinaten 
x en y van het punt P de rol van controle- 
variabelen spelen. De bifurcatieverzame- 
ling is diffeomorf met de cusp-katastrofe. 
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Dit kan men gemakkelijk inzien door de 
weglengtefunctie, f(x,y, 0), die hier de 
rol van potentiaalfunctie speelt, te ont- 
wikkelen om de coördinaten van het 
brandpunt (0 = 0, x = Ey = 0). 

Men kan ook de katastrofen van codimen- 
sie 3, dus de swallowstail en de elliptische 
en hyperbolische umbilic, zichtbaar maken 
als brandvlak. Dit doet men?’ door een 
laserstraal te laten vallen door een druppel 
die op een glasplaat is aangebracht, en 

het achter de glasplaat gevormde brand- 
vlak met een verplaatsbaar scherm af te 
tasten. De drie ruimtecoördinaten x,y 


en z spelen dan de rol van controlevariabe- 


len. 


Stabiliteit van schepen? 

De stabiliteit van een schip wordt in de 
eerste instantie bepaald door de positie 
van zijn zwaartepunt en van het aangrij- 
pingspunt van de opwaartse kracht, dus 
door de plaats van het zwaartepunt van 
het verplaatste water. Wanneer het schip 
kantelt of wanneer het wateroppervlak 
helt tijdens de passage van een golf, ver- 
plaatst het aangrijpingspunt van de op- 
waartse kracht zich langs een kromme 
B(6) ten opzichte van het schip (zie 
figuur 10). 


Voor een rechthoekig schip kan men deze 
kromme gemakkelijk bepalen. Bij kantelen 
van het schip over een hoek 0, wordt een 
sector van doorsnede 5 a? tg0 (waar 2a de 
breedte van het schip is) ondergedompeld, 
terwijl een diametraal daartegenover lig- 
gende sector van gelijke doorsnede droog- 
valt. Daardoor wordt het zwaartepunt van 
het verplaatste water verplaatst van de 
oorsprong, het punt B, naar een punt 
B(6) met coördinaten (x,y) op de para- 
bool: 


(41) 


waar p de kromtestraal in het punt B, 
voorstelt. Nu blijft de opwaartse kracht 
onder alle omstandigheden loodrecht op 
de kromme B(6) gericht, en de stabiliteit 
van het schip wordt daarom bepaald door 
de vraag of het door zwaartekracht en op- 
waartse kracht gevormde koppel het schip 
terugdraait naar de verticale stand of juist 
verder van de verticaal wegdraait, met an- 


y=x"|[2p, 
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Fig. 10. Zwaartepunt G en aangrijpings- 
punt B van de opwaartse kracht voor een 
rechthoekig schip. 


dere woorden of het zwaartepunt links of 
rechts van de normaal op B(O) ligt. 

De vergelijking van de normaal op de 
kromme B(6), aangrijpend in een punt 
met coördinaten x =ptg0,y= Dp tg20, 
is: 


(Y—y)tg0 + (A —x)=0, (42a) 
of: 
X+(Y-p)tgd -Fptg°0=0. (42) 


Deze vergelijking kunnen wij weer opvat- 





ten als de definitievergelijking van het 
evenwichtsoppervlak van een cusp-kata- 
strofe, waar 0 (of tg) de rol van toe- 
standsvariabele speelt, terwijl de coördi- 
naten (X, Y) als controlevariabelen fun- 
geren. 


Zet men Ô verticaal- en X en Y horizon- 
taal uit, dan krijgt men een regelvlak (zie 
figuur 1 la) gevormd door de verschil- 
lende normalen op B(6), ieder uitgezet op 
een hoogte 0 boven het (X, Y)-vlak, dus 
boven de doorsnede van het schip. Die 
lijnen waarvoor geldt dat (42) twee samen- 
vallende oplossingen heeft, dus waarvoor 
naast (42) ook geldt: 

3 2 
(Yp) Pte 0 =0, (43) 
vormen de raaklijnen aan de zogenaamde 
evoluente van de kromme B(0). Deze 
lijnen liggen verticaal boven een cusp- 
kromme, de bifurcatieverzameling 
27pX? =8(Y-p), die de scheiding aan- 
geeft tussen de mogelijke posities (X, Y) 
van het zwaartepunt G waarvoor het 
schip stabiel respectievelijk instabiel is. 
Het zal duidelijk zijn dat de vorm van het 
schip een grote invloed op de vorm van 





Fig. 11. Het evenwichtsoppervlak en 

de bifurcatieverzameling voor een recht- 
hoekig schip (a), en voor een elliptisch 
schip (b). 
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de kromme B(0), en daarbij op de vorm 
van de grenskromme tussen stabiliteits- 
en instabiliteitsgebied zal hebben. Zo is 

de kromme B(6) voor een schip met een 
elliptische doorsnede (een kano) zelf ook 
een ellips, waarvan de evoluente een omge- 
keerde cusp-kromme is, en in dit geval is 
juist het buitengebied van de cusp het 
instabiliteitsgebied waar de kano zal om- 
slaan (zie figuur 11b). 


Niet-lineaire oscillaties? »* 

In een groot aantal toepassingen heeft men 
te maken met (in het algemeen niet- 
lineaire) oscillatoren. Een van de een- 
voudigste modellen voor een niet-lineaire 
oscillator wordt voorgesteld door de Duf- 
fing-vergelijking: 


d?x/dt? + k(dx/dt) +x +ax® = 


=F cos Qt. (44) 
Als de dempingsconstante k,en de coëffi- 

ciënt a van de niet-lineaire term klein zijn 

en als Q dicht bij de resonantiefrequentie 

ligt, dus als: 

OQ=l+w, met k‚a,w 1, | (45) 
moet men oplossingen verwachten van de 

vorm 


x(t) = A cos (Q1-6) + B cos 3 (Qt-ô) 


En (46) 
Invullen van deze uitdrukking en gelijk- 
stellen van de coëfficiënten van cos Qt 


en sin Q\f in beide leden van (44) levert 
dan de condities: 


(1-Q? + À aA?) A cosô + 


+KkQA sins =F (47a) 
_kQA cos8 + (1-2? + 
+3 aA?) Asinô =0, (47b) 


en analoge condities voor de coëfficiën- 
ten van de hogere-ordetermen in cos 3 St, 
sin 3 Qt, enz. 

Uit (47) vindt men na invullen van (45) 
voor de fase: 

8 =arctg 4k/(3aA®* —8w), (48a) 
terwijl de amplitude 4 bepaald wordt 
door de conditie: 


3 
A°(G aA? Zeo = FP KLAP. (485) 


Deze 3°-graads vergelijking in A ° kan men 
weer opvatten als de vergelijking voor het 
evenwichtsoppervlak van een cusp-kata- 
strofe, waarbij A de rol speelt van toe- 
standsvariabele en w en a optreden als 
controlevariabelen. 

De bifurcatieverzameling die bepaald 
wordt door de conditie dat (48b) twee 
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Fig. 12. Evenwichtsoppervlak en bifur- 
catie verzameling voor de Duffing-oscil- 
lator. 


samenvallende wortels heeft, bestaat nu 
uit twee afzonderlijke cusp-lijnen (zie 
figuur 12) in het (a, w)-vlak. 


Voor a = O0 is de oscillator lineair en de 
amplitude is maximaal als aan de reso- 
nantievoorwaarde w = 0 is voldaan. Voor 
a > 0 is de oscillator ‘hard’, dat wil zeg- 
gen de terugdrijvende kracht neemt ster- 
ker dan lineair toe. In dat geval zal de 
amplitude A eerst geleidelijk toenemen 
met w, maar dan bij een bepaalde w 
sprongsgewijs afnemen tot op het onder- 
ste blad van het evenwichtsoppervlak. 
Tegelijk zal er een fasesprong optreden. 
Bij afnemende frequentie zal de amplitude 
weer terug omhoog springen, maar nu bij 
een lagere waarde van w, zodat hysterese 
optreedt. 

Een analoog effect, maar dan bij negatieve 
waarden van w, zal zich voordoen in het 
geval van een ‘zachte’ oscillator (a <0), 
waarvoor de terugdrijvende kracht zwak- 


ker dan lineair toeneemt met de uitwijking. 


Als wij in de Duffing-vergelijking, naast 
de 3°-graadsterm ax®, ook hogere-graads- 
termen opnemen, bijvoorbeeld een term 
bx’, krijgen wij te maken met hogere- 
ordecuspoïden, zoals bijvoorbeeld de 
swallowstail, waarvoor het evenwichts- 
oppervlak een 5°-graads‘oppervlak’ in de 
variabelen A °,w,a en b wordt, en de 
bifurcatieverzameling een oppervlak in de 
(w,a, b)-ruimte. 
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De potentiaaldistributietheorie en de 
statistische mechanica van vloeistoffen 


B. Widom 


Het principe van de potentiaaldistributie in de statistische mechanica zegt dat het kano- 
nische gemiddelde van de Boltzman-factor van de wisselwerkingspotentiaal, gepeild 
door een proefdeeltje, gelijk is aan het quotiënt van de dichtheid en de thermodyna- 
mische activiteit. Dit principe heeft vele uitbreidingen en toepassingen, in ruimtelijk 
homogene zowel als inhomogene vloeistoffen. In homogene fluïda leidt het principe 
tot betrekkingen tussen de paarcorrelatiefunctie en de thermodynamische eigenschap- 
pen; in ruimtelijk inhomogene fluïda verschaft het de meest directe weg tot de beschrij- 
ving in de gemiddelde-veldbenadering van het dichtheidsprofiel, in het grensvlak tussen 
twee fasen, en langs de curve waarbij drie fasen met elkaar in evenwicht zijn. 


Het is mijn bedoeling om hier een schets 
te geven van de ideeën en voornaamste 
toepassingen van wat bekend staat onder 
de naam van potentiaaldistributietheorie 
(of methode, respectievelijk theorema) in 
de statistische mechanica. Deze theorie is 
erg nuttig gebleken en is toegepast op bijna 
alle terreinen van de statistische mechanica 
van evenwicht. Het algemene principe is 
al bijna twintig jaar bekend !:? en speciale 
gevallen nog langer. 


De configurationele partitiefunctie van 
een systeem van N deeltjes in een volume 
V kan geschreven worden als: 


Ì 
On = ar en 


+WI/KT dr, … dins, (1) 


De dr; zijn volume-elementen, Uy —j is 
de wisselwerkingsenergie van de deeltjes 
1, …, N-I als een functie van hun con- 
figuratie, W is de extra wisselwerkings- 
energie die het gevolg is van de aanwezig- 
heid van het N* deeltje (het proefdeeltje) 
op een willekeurig vast punt in de vloei- 
stof. Vergelijking (1) kan dan herschreven 
worden als: 


On = zr VW-I!Oy1 (er WAT), 
(2) 


Hier duidt < * ) een kanonisch gemid- 
delde aan in de ruimte van N—1 deeltjes, 
waarbij de configuraties van de Nl 
deeltjes niet beïnvloed worden door de 
aanwezigheid van het proefdeeltje. Daar 
On 1/@n gelijk is aan À, de thermo- 
dynamische activiteit, en N/V aan p, de 
dichtheid, vindt men dat: 
pikt TTT, (3) 
dit is de grondvorm van ons principe. 

De reden voor de keuze van de naam 
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| ‘potentiaaldistributie theorie’ is als volgt: 


wanneer P(Y) de distributie is van de po- 
tentiaal, dat wil zeggen wanneer P(W)d WY 
de kans aangeeft dat de potentiaal die on- 
dervonden wordt door een proefdeeltje 
ligt tussen W en W +dW, dan kan het ka- 
nonisch gemiddelde < …- ) van elke func- 
tie van W berekend worden als f … P(Y) 
dW. Kiest men een benadering voor P(Y) 
dan verkrijgt men nuttige benaderings- 
methoden. 

Vergelijking (3) heeft een interessante in- 
verse die pas onlangs ontdekt is”. Laat ons 
Oy; schrijven als: 


Ì 
mann 


esp (-Uy-1/kT) dr, … drin 1 


On-1 = 


en: Pe 
WID! V ij 


fexp [—(Uy-W)/KT] dr, … dry 
Vv | 


ED oper 


(4) 


In deze vergelijking is het gemiddelde 
K=- ) verschillend van hetgeen we 
aangaven met { … ) :  …- ) iseen 
kanonisch gemiddelde in een systeem van 
N deeltjes, waarbij het deeltje dat de po- 
tentiaal W peilt er één is van de N,en niet 
een proefdeeltje, zodat de configuraties 
van de overige N—1 deeltjes nu beinvloed 
worden door de aanwezigheid van het N° 
deeltje. Omdat Oy _ 1/Qy weer gelijk is 
aan À en N/V aan p verkrijgen we uit (4): 


Npe Ke WE Po (5) 


wat de inverse is van (3). 
In een dichte vloeistof opgebouwd uit 


moleculen die elkaar krachtig afstoten op 
korte afstand is exp (— W/k7’) bijna altijd 





Dr. B. Widom (54) behaalde de doctors- 
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University in 1953, en was wetenschappe- 
lijk medewerker aan de Universiteit van 
North Carolina van 1952 tot 1954. Hij 
keerde naar Cornell terug in 1954 en is 
daar nu hoogleraar in de scheikunde. In 
het voorjaar van 1972 was hij als Van der 
Waals hoogleraar verbonden aan het Insti- 
tuut voor Theoretische Fysica en het Van 
der Waals Laboratorium aan de Universi- 
teit van Amsterdam. Hij werkt aan de 
thermodynamica en statistische mechanica 
van fasenevenwichten. 


erg klein, en dus is het niet praktisch om 
vergelijking (3) rechtstreeks toe te passen 
bij de berekening van de activiteit (ofwel 
de configurationele chemische potentiaal 
Hu = kT'In \) uitgaande van moleculaire 
configuraties die men bijvoorbeeld gene- 
reert in een computersimulatie; de con- 
figuraties die grote waarden van p/À op- 
leveren komen namelijk maar zelden voor. 
De inverse vergelijking (5) lijdt aan het 
tegenovergestelde probleem: exp(W/k7’) 
is bijna altijd enorm groot, maar ter com- 
pensatie is er een heel kleine gewichtsfac- 
tor in het gemiddelde « …-) ;en de 
waarde van \/p moet dan verkregen wor- 
den als het resultaat van een samenspel 
van heel grote en heel kleine factoren, 
hetgeen moeilijk nauwkeurig uit te voeren 
is. Gelukkig zijn er technieken ontwikkeld 
om deze moeilijkheden het hoofd te bie- 
den, en de hier geschetste ideeën, vooral 
vergelijking (3), hebben recentelijk 
waardevolle toepassing gevonden bij de 
berekening van chemische potentialen, 
zelfs in dichte vloeistoffen, uitgaande van 
met de computer gegenereerde configura- 
ties. We verwijzen naar het werk van 
Adams® voor een vloeistof van harde bol- 
len, en dat van Romano en Singer® voor 
modellen van vloeistoffen bestaande uit 
moleculen met realistische wisselwerkings- 
krachten; Powles® paste deze methode 
toe ter berekening van de dampspannings- 
kromme, en Shing en Gubbins’ boekten 
belangrijke vooruitgang door de ideeën 
van de potentiaaldistributietheorie te 
combineren met de techniek van ““um- 
brella sampling” in Monte Carlo-bereke- 
ningen. 


Er zijn ook andere nuttige resultaten!» 2 

bekend van dezelfde soort als (3). Bij- 

voorbeeld, in het geval dat de vloeistof 
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een mengsel van twee of meer componen- 
ten a, b, …, is, kan (3) toegepast worden 
op elk van de afzonderlijke componen- 
ten, zodat: 


Pala = (exp (-Wa/kT)), 


enzovoort. 


(6) 


WY, isde potentiaal die gevoeld wordt door 
een proefdeeltje dat zelf een molecule van 
de soort a is, terwijl het gemiddelde { …- Ò 
berekend wordt over de configuraties van 
alle moleculen in het mengsel. Een ander 
voorbeeld is dat van een eencomponentig 
systeem van moleculen waarin de energie 
additief is (een som van paarpotentialen); 
de gemiddelde potentiële energiedicht- 
heid, u, wordt dan gegeven door! : 
u=iN{ We PIT), (7) 
Een andere soortgelijke relatie levert de 
radiële distributiefunctie g(r), of de paar- 
distributiefunctie en de paarcorrelatie- 


functie die er nauw verband mee houden. 


De eigenschappen van deze functies zijn 
uitvoerig behandeld in het historisch over- 
zicht van De Boer? , die ook liet zien hoe 
wezenlijk de rol is die deze functies spe- 
len in de theorie van vloeistoffen. Laten 
we uitgaan van de definitie van g(r), 
namelijk dat p g(r) de gemiddelde lokale 
dichtheid is op afstand r van een mole- 
cule. Gebruiken we dan een soortgelijk 
argument als waarmee (3) werd afgeleid, 
dan vinden wel,? dat g(r) gegeven wordt 
door: 


g(r) PONT = D/p)? 


( el r)+ wr rk T ). ® 


In deze vergelijking zijn r en r’ +r de 
plaatscoördinaten van een tweetal proef- 
deeltjes die zich op een vaste afstand 

r= |r | van elkaar bevinden. De totale 
wisselwerkingspotentiaal is weer een som 
van paarpotentialen ®(7;;). 
Het produkt: 
y(r)=g(r)eP MIET, (9) 
dat optreedt in (8), is van fundamenteel 
belang, vooral in vloeistoffen van harde 
bollen. Deze grootheid is namelijk continu 
bij r = o (o is de doorsnee van de bol), 
ofschoon g(r) en exp [é(r)/KT] geen van 
beide continu zijn. Uit (8) en (9) verkrijgen 
wij dan: 

y(O) = (Mp)? Le PAT). (10) 
Voor harde bollen kan exp (—W/KT) 
alleen maar de waarde 0 of 1 hebben, dus 
zijn exp (—-W/KT)een exp (-2W/KT) 
identiek, en verkrijgen we uit (10) en (3): 
y(O)=Np (harde bollen), (11) 


welke uitdrukking werd afgeleid door 
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Hoover en Poirier”. 


Voor meer realistische paarpotentialen 
leidt vergelijking (10), met (3) en (7), tot 
een interessante ongelijkheid. Laten we 


als middelwaarde A van een willekeurige 


grootheid A definiëren: 
Ast Ae VE Ge TE Pe {I2) 


waarin < …” ) hetzelfde kanonisch ge- 
middelde is dat we reeds eerder gebruikt 
hebben. (Wanneer we voor A de grootheid 
A =exp (W/KT) kiezen, dan is de middel- 
waarde die we met een horizontale streep 
aanduiden dezelfde die werd aangegeven 
door « … ) in (4) en (5).) Dan ver- 
krijgen we uit (3), (7), (10) en de defini- 
tie (12): 


y(0) = Np) exp -2u/pkT) 


exp (-W/KT) exp (U/KT). (13) 
Er geldt echter voor elke middelwaarde 
van elke willekeurige grootheid dat de 
middelwaarde van de e-macht van die 
grootheid tenminste even groot is als de 
e-macht van de middelwaarde. Dus vinden 
we: 

y(0) > (Mp) eFW/PET, (14) 
welke ongelijkheid de grootheid y (0) in 
verband brengt met een zuiver thermo- 
dynamische eigenschap van een vloeistof, 
namelijk met u/p, de gemiddelde poten- 
tiële energie per molecule. Henderson'® 
heeft de waarde van het rechterlid van — 
(14) geschat en vond dat deze van de orde 
van 10° was voor een dichte vloeistof; 
volgens vergelijking (14) is deze grote 
waarde nog maar een ondergrens voor 

y (0). 


De zeer grote waarde van y (0) is vooral 
het gevolg van de korte afstandsafstotings- 
krachten tussen de moleculen in de dichte 
vloeistof. In de buurt van een kr:tiek 
punt, waar de vloeistof lang niet zo’n 
hoge dichtheid heeft, is y (0) veel kleiner. 
Een schatting kan als volgt verkregen wor- 
den. Denk aan de definitie van het roos- 
tergas: het volume is in cellen verdeeld, er 
is niet meer dan één molecule per cel, en 
tussen elk paar moleculen in nabuur- 
cellen bestaat een constante negatieve 
wisselwerkingsenergie. Als elke cel c buren 
heeft (zodat c het coördinatiegetal is in 
het model), dan vinden we in de benade- 
ring van Bethe en Guggenheim (quasi- 
chemische benadering): 


6 


ld 
nn (5) | 
kr.pt. c—-2 


De waarde hiervan is slechts 35 voor c =6, 
terwijl voor hogere c-waarden deze uit- 
drukking naar de waarde 2 gaat. Dus, 
zoals we al vermoedden, is er geen reden 


om aan te nemen dat y (0) erg groot is bij 
het kritieke punt. 


Onder de meest interessante toepassingen 
van de potentiaaldistributietheorie vallen 
diegenen waarbij benaderingen gevonden 
worden voor de eigenschappen van vloei- 
stoffen met behulp van de gemiddelde- 
veldbenadering (““mean-field approxima- 
tion”), welke reeds ten grondslag ligt aan 
de toestandsvergelijking van Van der 
Waals. De potentiaaldistributietheorie 


| blijkt de meest directe weg te zijn die leidt 


tot de toestandsvergelijking van de homo- 
gene vloeistof en ook tot de structuur van 
het grensvlak tussen coëxisterende fasen 
in heterogeen evenwicht in beide gevallen 
in de gemiddelde-veldbenadering, waarbij 
Van der Waals zelf reeds de wezenlijke 
uitkomsten verkreeg. De Boer!! heeft 


_een helder en gezaghebbend verslag over 


de gemiddelde-veldbenadering geschreven, 
waarin alle essentiële resultaten voor ho- 
mogene systemen gevonden kunnen wor- 
den; in zijn boeiende historische schets 
van het werk van Van der Waals, en van 
de recente herleving van de belangstelling 
ervoor, heeft de Boer aangetoond!? hoe 
de theorie van inhomogene vloeistoffen 
(de theorie van de capillariteit) voor Van 
der Waals een natuurlijke uitbreiding was 
van de inzichten, die bij hem leidden tot 
zijn toestandsvergelijking. Dit verband 
blijkt ook duidelijk in de formulering van 
van Kampen!* en komt eveneens naar vo- 
ren in de potentiaaldistributietheorie. 


Laten we ons de moleculen van een vloei- 
stof voorstellen als harde bollen die elkaar 
aantrekken. We zullen de afstoting exact 
behandelen maar de aantrekking in de 
gemiddelde-veldbenadering. We beginnen 
met de homogene vloeistof. Stel dat w 

de waarschijnlijkheid is dat de harde kern 
van een proefdeeltje geplaatst kan worden 
op een willekeurige positie in de vloeistof 
zonder dat het de harde kernen van een 
van de andere reeds aanwezige moleculen 
overlapt, en als W,, het attractieve deel 
van de potentiaal is die door het proef- 
deeltje ondervonden wordt, gaat (3) over 
in: 

pl\=w {exp (-W‚,./KT) ). (15) 
In de gemiddelde-veldbenadering’* ver- 
vangen we nu het gemiddelde van de e- 
macht in (15) door exp (—-W_#/KT), 
waarbij W_ # een constante attractieve 
wisselwerkingsenergie voorstelt die het 
proefdeeltje overal ondervindt waar het in 
de vloeistof past: 


Worse f Dory, W)dr=-2ap. (16) 
v>o0 


In deze vergelijking is &,,‚‚ (”) het attrac- 
tieve gedeelte van de intermoleculaire po- 
tentiaal, ov is weer de diameter van de 
harde kern, en: 
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(17) 


is de a van Van der Waals. De waarschijn- 
lijkheid w kan uit (15) herkend worden 
als de waarde van p/À voor een vloeistof 
van harde bollen zonder attractie; noem 
deze waarde p/\ p, waarbij Anp = Any (£) 
berekend moet worden bij de dichtheid p 
van ons model. Dan wordt de toestands- 
vergelijking van de vloeistof in deze bena- 
dering gegeven door: 


a= —5 S Ü ttr (7) dr, 
r>o 


Mo, T) = App (p)eZLPIKT, (18) 
Deze uitdrukking is equivalent met de be- 
kende vergelijking p = php — ap° voor de 
druk p. 

We hebben deze behandeling van de ho- 
mogene vloeistof gegeven om te laten zien 
hoe de potentiaaldistributietheorie kan 
worden toegepast; het is echter van meer 
belang hem toe te passen op de inhomo- 
gene vloeistof, in het bijzonder op de be- 
paling van het dichtheidsprofiel p (z) als 
een functie van de hoogte z, bijvoorbeeld 
in een grensvlak vloeistof-damp.* Verge- 
lijking (3) moet dan gegeneraliseerd wor- 
den tot!®: 

plrijk= Le VIRE |, (19) 
waarbij het gemiddelde nu niet meer ge- 
nomen wordt op een willekeurig punt in 
de vloeistof, maar op de plaats r, zoals 
aangegeven; maar À stelt nog steeds de 
ruimtelijk-homogene thermodynamische 
activiteit voor. Vergelijking (19) geldt al- 
gemeen, voor elke vorm van inhomogeni- 
teit. In het geval van een horizontaal grens- 
vlak mag r vervangen worden door z. 

In de gemiddelde-veldbenadering is het 
rechterlid van (19) een eenvoudige func- 
tionaal van de dichtheidsverdeling p(r), 
zoals we direct zullen zien. Daar À dan uni- 
form is, gaat (19) over in een functionaal- 
vergelijking voor p(r), of meer in het bij- 
zonder voor p(z), in het grensvlak. 

Het gemiddelde veld wordt nu r-afhanke- 
lijk, zodat (16) vervangen moet worden 
door: 


Par) = ks Pater CI o(r+r) dr, 
(20) 


dus een functionaal van de dichtheidsver- 
deling. Als p(r) weinig verandert over af- 
standen van de orde o geldt nog steeds: 
A= App lor)lexp [Yr P)/KT), (25) 
in de gemiddelde-veldbenadering, net als 
in (16)-(18). Maar nu is A » lp (r)] gelijk 
aan Ànp (p), berekend voor de lokale 
dichtheid p = p(r), terwijl de A zonder 
aangegeven argument nog steeds de ruimte- 
lijk-homogene thermodynamische activi- 
teit aanduidt. Als p (7) nu ook weinig ver- 
andert over de werkingssfeer van ®,,,, (0), 
dan kunnen we (20) vereenvoudigen. We 
zullen hiertoe spoedig overgaan, maar 
voorlopig handhaven we (20) nog in zijn 
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meer algemene, niet-lokale gedaante. 


Als we nu (17), (18) en (20) in rekening 
brengen, en de configurationele chemische 
potentiaal u = kT In A invoeren, dan kun- 
nen we (21) herschrijven als: 


—S Pater) lor +7) POldr, = 
=ulo(r), T|— u. (22) 


Hierin is ulp(r), T] = KT In Ip (r), 7] 
terwijl \(p, 7) gegeven wordt door (18). 
Dus kan u[o(r), 7] beschouwd worden 
als de chemische potentiaal van de hypo- 
thetische vloeistof die homogeen gehou- 
den wordt bij de dichtheid welke gelijk is 
aan de lokale dichtheid van de inhomo- 
gene vloeistof; de u zonder argumenten is 
weer de ruimtelijk-homogene chemische 
potentiaal van de inhomogene vloeistof. 


Vergelijking (22) is een functionaalverge- 
lijking voor de dichtheidsverdeling p(r) 
die wij zochten. In een grensvlak moet 
p(r) opgevat worden als p(z), en dan moet 
(22) opgelost worden onder de condities 
p(ee) = pg en p(—ee) = pj, waarbij pg en 
py de dichtheden zijn van de gas- en de 
vloeistoffase bij de temperatuur 7 (of, 
wat daarmee gelijkstaat, bij de chemische 
potentiaal u, die zelf een functie van T'is). 


Als p(r) of p(z) inderdaad slechts weinig 
verandert over de werkingssfeer van 
Pater 7), dan mag p (r_ +7’) in (22) ont- 
wikkeld worden rond r,‚en de ontwikke- 
ling afgebroken worden na de tweede 
orde; wanneer de inhomogeniteit alleen 
in de z-richting optreedt, verkrijgen we 
dan: 


md*p(z)/dz? =ulp(2),T]-u, (23) 

waarin m het tweede moment van 

—Ôs ttr 7): 

mei FP Oster Dr, (24) 
r>o 


is, dus een positieve constante. Dit zijn de 
beroemde formules van Van der Waals 
voor de bepaling van het dichtheidspro- 
fiel van een grensvlak vloeistof-damp. 

In de buurt van het kritieke punt mogen 
we u(Pp, 7) — uKOontwikkelen naar machten 
van p—pe en 7-7,, de verschillen van p 
en 7 van hun respectievelijke waarden bij 
het kritieke punt. Als we nu weer afbre- 
ken na de tweede term, vinden we: 


(Pp, T) —u=-— A(T.-T) (pp) + 


+B(o-p). (25) 


Uit (25) en (23) verkrijgen we dan het _ 
klassieke hyperbolische tangensprofiel!®: 


p(@) — pe =VA(T-TDIIB 
tanh [VA(T,-T)/2m zl. (26) 


Op dezelfde wijze kunnen we het rooster- 
gasmodel analyseren, dat we reeds eerder 
hebben ingevoerd. Noem —e (e > 0) de 
wisselwerkingsenergie tussen moleculen in 
nabuurcellen. Neem aan dat op het niveau 
z elke cel c° buren boven zich heeft op het 
niveau z + 1,en evenzo c' buren beneden 
zich op het niveau z — 1 „en dan nog c—2c' 
buren in zijn eigen laag, op het niveau z, 
zodat het totale coördinatiegetal weer c is. 
(Merk op dat z in dit voorbeeld dimensie- 
loos is; het is de hoogte gemeten in een- 
heden van de hoogte van een cel.) Het 
model laat geen meervoudige celbezetting 
toe, dus de potentiaal W die ondervon- 
den wordt door een proefdeeltje is + oo 
als er al een molecule van de vloeistof in 
dezelfde cel zit. Dan gaat (19) over in: 


DON = [IO] exp War KD), 
(27) 


=[1-p(@lexp [WOI KT], 
(28) 


waarbij de factor 1—p(z) de waarschijn- 
lijkheid is dat het proefdeeltje een wille- 
keurige cel op het niveau z leeg aantreft. 
(Deze dichtheid is weer dimensieloos; het 
is de deeltjesdichtheid in eenheden van 
het reciproke volume van een cel.) Verge- 
lijking (27) is exact, en (28) is de gemid- 
delde-veldbenadering. Opnieuw wordt 
Wor gegeven door (20), zodat we voor 
het roostergas vinden dat: 


Wine (2) = —elc'o(Z+1) + (e-2c)p(2)+ 


+c'p(e-1) (29) | 


=—elep(z)+e' Ap(z)l, 


waarbij A° de tweede-differentie operator 
is. 

Als de vloeistof homogeen is, met een uni- 
forme dichtheid p(z) =p,en met A°p =0, 
dan verkrijgt men uit (28) en (29) de toe- 
standsvergelijking: | 
Np, T)=[p/(l—pylerctPIEE, (30) 
voor het roostergas in de gemiddelde- 
veldbenadering. Uitgedrukt in de ermee 
corresponderende configurationele che- 
mische potentiaal u(p, 7) geven de ver- 
gelijkingen (28)-(30) dan voor p (z) van 
de heterogene vloeistof de differentie- 
vergelijking: 

ce A°p(z)=ulp(z), Tu, (31) 
waarin ulo (z), 7] — u dezelfde betekenis 
heeft als tevoren (zie de discussie na ver- 
gelijking (22)). Dit is het analogon van 

de niet-lokale integraalvergelijking (22) 

in het model van harde bollen met attrac- 
tie. Als p (z) langzaam variëert wordt de 
tweede differentie A° p(z) gelijk aan de 
tweede afgeleide d° p(z)/dz®, zodat (31) 
het analogon wordt van (23), met m = ce 


(wat nu de dimensie van een energie heeft 
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omdat de huidige z en p dimensieloos 
zijn). 

Deze techniek, waarbij grensvlakprofielen 
afgeleid worden-met behulp van het po- 
tentiaaldistributieprincipe door 

(exp (—-W/KT) ), in gemiddelde-veld- 
benadering te schrijven als een functionaal 


van p(r) heeft onlangs een nieuwe toepas- 


sing gevonden; namelijk op de bepaling 
van het concentratieprofiel bij een drie- 
fasenlijn. In sommige gevallen zijn drie 
vloeistoffasen met elkaar in evenwicht in 
drie horizontale lagen, paarsgewijs van el- 
kaar gescheiden door twee gewone twee-. 
fasengrensvlakken, een tussen bovenste en 
middelste laag, de ander tussen middelste 
en onderste laag. Even dikwijls echter 
nemen de fasen in drie-fasenevenwicht 
niet de vorm aan van horizontale lagen; zij 
ontmoeten elkaar daarentegen met eindige 
contacthoeken langs een lijn — de drie- 
fasenlijn. In dit geval zijn zij paarsgewijs 
van elkaar gescheiden door drie grens- 
vlakken, elk tussen twee fasen, waarbij 
die drie grensvlakken ook op de drie- 
fasenlijn samenkomen. Ver van de drie- 
fasenlijn vandaan hebben deze vlakke, of 
bijna vlakke scheidingsoppervlakken 
dichtheids- of concentratieprofielen net 
als bij een gewoon grensvlak; die pro- 
fielen kunnen beschreven worden in de 
gemiddelde-veldbenadering, met behulp 
van de theorie die we reeds geschetst heb- 
ben. Nu willen we echter ook de ruimte- 
lijke variatie van de dichtheden bepalen 
dicht bij de drie-fasenlijn, net zoals we de 
variatie van de dichtheid door elk van de 
grensvlakken tussen twee fasen reeds be- 
paald hebben. 

We beschouwen weer een eenvoudig mo- 
del, namelijk een mengsel van drie. compo- 
nenten a, ben c zodanig dat de moleculen 


van dezelfde soort geen wisselwerking met _ 


elkaar hebben terwijl paren van ongelijke 
moleculen elkaar afstoten als harde bol- 
len met doorsnede 0. Dit is het model van 
Guerrero et al}? We behandelen eerst een 
homogene fase van deze modelvloeistof 
en vinden À, (a, Pp» Pc) enz, uit (6) in 
de gemiddelde-veldbenadering. De inter- 
acties zijn zodanig gekozen dat exp(—Y,/ 
kT) alleen maar O of 1 is; de waarde 0 
wordt verkregen als een proefdeeltje van 
soort a binnen de afstand o komt van de 
oorsprong van een deeltje van type b of c 
in de vloeistof; zoniet, dan heeft de 
e-macht de waarde 1. Dus is p„ /)\g de 
waarschijnlijkheid dat een willekeurig 
punt in de vloeistof niet binnen de afstand 
o ligt van het centrum van enig molecule 
van het type b of c. We benaderen deze 
waarschijnlijkheid door wat hij zou zijn als 
de b- en c-moleculen lukraak geplaatst 
zouden zijn (hetgeen, gezien het feit dat 
ze WASSER niet Seat waar kan Rin: 
Dan vinden we!’ 


Pal Xa >exP Eik +), enz. GD 
waarbij pa … en \ … alle gemeten worden 
in eenheden (4/3) (20)®. Dit vervangt 


nu het eerder verkregen resultaat van de 
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gemiddelde-veldbenadering, omdat er in 
dit geval geen aantrekkingskrachten tus- 
sen de moleculen zijn, alleen maar afsto- 
tingen tussen harde bollen. De benadering 
(32) is thermodynamisch consistent; hij 
voldoet aan de strenge eisen ò In A,/dpp = 
ò In \,/òp,, enz, waaraan andere denk- 
bare benaderingen mogelijkerwijs niet zou- 
den voldoen. Dus vinden we voor de toe- 
standsvergelijking van de vloeistof, in - 
deze benadering: 


Àa Pas Abs Pc) = Pa EXP (Pv +Pc), enz.(33) 


Met behulp van (33) kan men toestanden 
vinden!” waarin het mengsel zich in drie 
fasen scheidt; in een van deze toestanden 
iS Pa > Pp + Pe; de andere twee worden 


uit de eerste verkregen door a, b en,c te 


permuteren. Vanwege de symmetrie moet 
in zulk een drie-fasenevenwicht gelden 
dat As = Ap = Àc. En ook alweer vanwege 
de symmetrie moeten de fasen elkaar langs 
een drie-fasenlijn ontmoeten, met drie 
gelijke contacthoeken (de hoeken tussen 
de twee-fasenvlakken) van 120°. 7 
Nu willen we het dichtheidsprofiel p‚ (x,y), 
Pv (X,y) en Pp (x,y) bepalen, waarbij x, 
de twee ruimtecoördinaten zijn die lood- 
recht staan op de richting van de drie- 
fasenlijn (er verandert natuurlijk niets in 
de richting langs de lijn, die de rich- 

ting is, evenwijdig aan alle drie de twee- 
fasenvlakken tegelijk). Daartoe beschou- 
wen we (6) in de vorm: 


Pa (Xx, Na = (exp (-Va/kT) bea 


die van toepassing is voor een inhomo- 
gene vloeistof, en we benaderen het 
rechterlid weer door te veronderstellen 
dat een proefdeeltje van type a deeltjes 
van type b of c alleen maar ziet alsof ze 
willekeurig geplaatst waren, maar met lo- 
kale dichtheden py (x;y)en po (x,y). 
Het resultaat is een benadering voor 
(exp(— Wa/kT')) „ in de vorm van de 
e-macht met in de exponent een integraal 
functionaal van py (x,y) en p. (x,y) die 
sterk analoog is met de functionaal (20), 
maar die een zuiver geometrische factor 


| bevat in plaats van de intermoleculaire 


potentiaal!®. 

Ik zal dit resultaat hier niet weergeven, 
want het zal binnenkort beschikbaar zijn 
voor de belangstellende lezer '*, maar ik 
zal de vereenvoudigde uitdrukking geven 
welke daaruit volgt in de veronderstelling 
dat de dichtheden langzaam veranderen. 


We mogen dan de integraal functionaal ont- | 


wikkelen en de ontwikkeling afbreken na 
de tweede term, net zoals we (23) uit (22) 
verkregen. Dan vinden we dat de dicht- 
heden pa (x,y), enz., door de volgende 
vergelijking bepaald worden: 


5 Vv. (pp + Pc) = In [Aa (Da, Pps Pe) AL, 


enz., (34) 
waar V? de twee-dimensionale Laplaciaan 
02 /ox? +9? /dyp* is, waar x en y gemeten 
worden in eenheden o en dus dimensieloos 
zijn; en waar À, (Pa, Pp, Pc) de door (33) 


ongelijksoortige moleculen, en van de 


bepaalde functie is, terwijl À de gemeen- 
schappelijke constante waarde van de acti- 
viteiten À,, Ap en A, in het drie-fasen- 
evenwicht is. Vergelijking (34) is duidelijk 
analoog aan (23) of (31); nu echter is de 
activiteit Aa (Pax, Pp Pc) van het hypothe- 
tische systeem met uniforme dichtheden 
Pa, -… Enz. gelijk aan de lokale dichtheden 
Pa (x,y), enz. via (34) afhankelijk van de 
ruimtelijke variaties van de dichtheden 
van de andere componenten dan a, enz. 
Dit is een gevolg van de afstoting tussen 


afwezigheid van een wisselwerking tussen | 
gelijksoortige moleculen in dit model. De 
fysisch-relevante oplossingen van (34), 

Pa (x,y) enz, zijn di welke overgaan in 
de passende constante dichtheden p„, 
enz. binnen elke fase. Ver van de drie- 
fasenlijn vandaan, en in de juiste richtingen 
bekeken, gaan deze oplossingen ook op 
correcte wijze over in de dichtheidspro- 
fielen aan de grensvlakken tussen twee 
fasen. 

Deze nieuwe toepassing maakt het nut 
van de potentiaaldistributietheorie duide- 
lijk. Ofschoon de theorie, zoals we reeds 
gezien hebben, de grondslag vormt voor 
belangrijke praktische methoden in de 
theorie van homogene vloeistoffen, is het 
waarschijnlijk juist in de toepassing op 
niet-homogene vloeistoffen dat zijn ware 
kracht en veelzijdigheid tot hun recht 
komen. 
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